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CHAPTER  1 


INTRODUCTION 

— J  The  importance  of  designing  control  systems  which  are  robust  or 
insensitive  to  variations  in  the  plant  parameters  has  long  been  appreciated. 
However,  the  rapid  advances  in  design  techniques  for  multivariable  systems  has 
heightened  interest  in  the  study  and  design  of  robust  systems.  The  purpose 
of  this  report  is  to  provide  an  up  to  date  survey  of  the  work  in  this  field 
and  summarize  the  results  of  research  in  this  area  conducted  at  the  Coordinated 
Science  Laboratory. 

The  report  begins  in  Chapter  2  with  a  description  and  examples  of 
the  robust  control  problem.  Chapter  3  provides  a  survey  of  research  in  the 
field  of  robust  control.  It  is  apparent  from  this  survey  that  the  work  can 
be  divided  into  two  areas.  The  first  assumes  unstructured  perturbations  and 
analyzes  worst  case  effects.  The  second  considers  large,  structured  perturba¬ 
tions.  The  parameter  space  design  method  presented  in  Chapter  4  is  directed 
at  the  second  area.  The  tools  of  Chapter  4  are  applied  to  a  fighter  aircraft 
example  in  Chapter  5.  Chapter  6  presents  an  optimization  approach  to  the  same 
problem.  Finally,  Chapter  7  summarizes  the  report  and  presents  several 
directions  for  future  research. 


The  objective  of  a  basic  control  system  design  problem  is  to  satis 
a  set  of  performance  specifications  for  a  given  dynamical  system.  The  robus 
control  system  design  problem  adds  to  the  basic  control  problem  by  requiring 
the  performance  criteria  to  be  satisfied  under  a  specified  class  of  perturba 
tions  to  the  dynamical  system.  Typical  examples  of  performance  criteria  for 
which  the  overall  closed  loop  system  must  be  robust  are: 

1)  Stability  or  nice  stability  (e.g.  defined  by  constraints  on 
eigenvalue  locations). 

2)  Limited  deterioration  of  a  performance  index. 

3)  Limited  deviation  from  an  ideal  behavior,  e.g.  constraints 
on  step  responses  or  frequency  responses  or  on  the  return 
difference. 

4)  Limited  deviation  from  a  reference  behavior,  e.g.  deviation 
from  a  nominal  trajectory  or  a  reference  model  response. 

5)  Tracking,  i.e.  zero  asymptotic  error  for  a  class  of 
reference  and  disturbance  inputs. 

6)  Limited  demand  on  control  tul  and  control  rate  lul. 

The  classes  of  perturbations  which  are  considered  can  be  grouped 
in  two  categories:  structural  perturbations  and  system  parameter  perturba¬ 
tions.  Some  examples  of  structural  perturbations  are: 
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1)  Sensor  failures. 

2)  Actuator  failures. 

3)  Switching  from  automatic  to  manual  control. 

4)  Change  in  system  order  due  to  a  failure.  Example:  An 
aggregate  description  for  several  power  generators  or  a 
traffic  flow  or  economic  variables  must  be  dissolved  into  a 
more  detailed  description  of  transients  between  individual 
components  in  failure  situations. 


Parametric  perturbations  are  due  to  uncertainties  in  the  plant  model  and  in 
the  controller  implementation.  Examples  are: 

5)  Analytically  known  dependence  of  a  plant  model  on  uncertain 
physical  parameters.  Example:  The  linearized  equations  of 
a  crane  with  physical  parameters  mc=crab  mass,  m^  =  load  mass, 
.4=  rope  length,  g  =  gravitational  constant,  and  state  variables 
x^  =  crab  position,  X£ = crab  velocity,  x^  =  rope  angle  and  x^  = 
rope  angular  velocity  are 


p 

i 

0 

0“ 

”0  ~ 

j  0 

0 

mLg/mc 

0 

1 

1 

X  =  ’ 

;  0 

0 

0 

1 

xH - 

—  m 

c 

0 

:  o 

0 

2 

-OJ 

0 

(1) 


with  m  = (mc + m^) g/mc &.  Input  u  is  the  force  accelerating 
the  crab.  The  crane  may  operate  with  an  unknow  load  mass 
m^  between  the  empty  hook  and  the  maximum  mass,  for  which 
the  crane  is  designed.  It  may  also  operate  with  an  unknow 
constant  rope  length  between  zero  and  the  height  of  the  crane. 
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6)  Numerically  known  dependence  of  a  plant  model  on  an  unknown 
physical  parameter  vector  Example:  linearized  equations 
of  longitudinal  motion  of  an  aircraft  depending  on  altitude 
and  speed 

x  =  A  (9)  +  B(Q)u  (2) 

A  and  B  may  be  given  for  J  typical  flight  conditions  in  the 
flight  envelope.  Aj=A(0j),  Bj3SB.(£j)»  j=l,2,...J. 

7)  Known  dynamics,  which  have  disappeared  in  a  simplified  design 
model  by  linearization,  truncation  of  structural  modes,  model 
reduction,  neglecting  of  actuator  and  sensor  dynamics.  In 
some  cases  it  may  be  possible  to  pull  out  all  uncertainties 
as  illustrated  by  Fig.  2 . 1,  where  for  P  =  0  the  nominal  plant 
N  is  obtained. 


Figure  2.1.  Nominal  plant  N  with  perturbations  P. 

In  simple  cases  the  perturbations  P  can  be  expressed  as  a 
diagonal  matrix  of  linear  or  nonlinear  operators. 

8)  Unknown  dynamics,  which  cannot  be  modeled.  In  this  case  only 
vague  assumptions  about  perturbations  6A,  6B,  5C  of  the 
system  matrices  A,  B,  C(x=Ax  +  Bu,  y  =  Cx)  or  perturbations 
6G(s)  of  the  transfer  function  matrix  G(s)  can  be  made. 
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9)  Quantization  effects  and  time  delays  in  controller 
implementation. 

10)  Variance  of  components  in  mass  produced  control  systems 
and  circuits. 

These  examples  of  system  properties  and  perturbations  show  that  many  special 
combinations  can  be  specified.  Therefore  many  different  definitions  of 
"robust  control"  can  be  found  in  the  literature. 

The  design  problem  for  a  robust  control  system  may  be  formulated  in 
one  of  the  following  three  forms; 

1.  Given  a  system  property,  determine  the  class  of  perturbations 
with  respect  to  which  the  system  property  is  robust.  Design 
the  controller  such  that  the  class  of  admissible  perturbations 
is  extended  in  the  direction  of  the  really  expected  perturba¬ 
tions  . 

2.  Given  a  class  of  perturbations,  determine  the  maximum  deviation 
from  a  desired  system  behavior  which  occurs  under  the  worst 
perturbation  in  the  given  class.  Design  the  controller  such 
that  the  maximum  deviation  is  minimized. 

3.  Given  a  system  property  and  a  cLass  of  perturbations,  determine 
if  there  exists  a  set  of  controllers  for  which  the  system 
property  is  robust  under  the  class  of  perturbations.  If  yes, 
select  one  on  other  criteria  than  robustness.  If  not,  relax 


specifications . 
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2.2.  Controller  Structure  and  Other  Design  Considerations 

Design  problems  for  control  systems  are  usually  parameterized  by  the 
assumption  of  a  controller  structure  which  defines  a  vector  of  design  para¬ 
meters.  Two  typical  assumptions  for  the  controller  structure  are  adaptive 
controllers  or  fixed  gain  controllers.  One  extreme  is  the  attempt  to  obtain 
as  much  information  about  the  perturbations  as  possible  by  on-line  identifica¬ 
tion  and  failure  detection.  Then  ideally  the  structure  and  parameters  of  the 
controller  are  adapted  in  order  to  achieve  the  best  possible  performance  of 
the  control  system  given  the  momentarily  available  information.  An  intrinsic 
difficulty  of  this  approach  is  that  plant  inputs,  which  admit  a  fast  and 
accurate  identification,  are  not  good  to  achieve  the  best  performance  and  vice 
versa.  Also  a  tradeoff  between  a  fast  failure  detection,  identification  and 
adaptation  and  a  reliab le  one,  which  avoids  false  alarms  and  noise  sensitivity 
of  the  adaptation,  must  be  made. 

The  other  extreme  is  the  attempt  to  find  a  fixed  gain  controller 
which  accomodates  a  specified  class  of  perturbations.  In  this  approach  it 
may  be  necessary  to  sacrifice  some  performance  in  the  nominal  case  in  order 
to  achieve  robustness  for  the  perturbed  situation.  The  assumed  controller 
structure  may  be  state  feedback,  or  static  or  dynamic  output  feedback.  Note 
that  full  state  feedback  is  not  the  most  general  controller.  Information 
about  the  unknown  parameters  is  contained  in  past  states,  their  processing  in 
a  dynamic  system  can  therefore  improve  the  performance.  In  some  cases  the 
unknown  parameters  can  be  introduced  as  additional  states,  which  may  be 
estimated  and  fed  back.  Most  design  techniques  are  restricted  to  linear 
systems,  thus  a  linear  controller  is  usually  assumed. 

Only  the  case  of  a  fixed  gain  controller  is  usually  called  robust 
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control.  However  it  should  be  apparent  from  the  previous  discussion  that 
robustness  is  also  a  desirable  feature  for  an  adaptive  control  system.  The 
fixed  gain  solution  indicates  whether  a  more  complex  adaptive  system  is  needed 
at  all,  or  how  far  one  has  to  go  adaptive.  Practical  solutions  to  the  robust¬ 
ness  problem  frequently  are  in  between  the  two  extremes  and  combine  features 
of  both  cases.  Examples  are: 

1.  Gain  scheduling  with  switching  or  continuous  variation  between 
fixed  linear  feedbacks  dependent  on  a  measurement  of  an 
environment  condition,  which  has  an  influence  on  the  plant 
parameter  values. 

2.  Variable  structures  with  state  dependent  switching  between 
fixed  linear  feedbacks. 

3.  A  fixed  gain  robust  controller  may  be  used  for  stabilization, 
with  an  adaptive  controller  then  used  for  improvement  of 
performance.  (Some  techniques  for  design  of  adaptive 
controllers  assume  an  open  loop  stable  system). 

4.  A  fixed  gain  robust  controller  may  be  used  as  backup  for  the 
case  of  a  failure  in  the  adaptive  system  or  in  a  gain 
scheduling  system.  Air  data  measurements  (e.g.  dynamic 
pressure)  are  not  very  reliable. 

5.  Under  external  noise  an  adaptive  system  may  not  adjust  fast 
enough  to  a  fast  change  in  plant  parameters.  In  such  cases, 
it  may  be  possible  to  switch  to  a  fixed  gain  robust  system 
until  the  identification  has  followed  and  adaptation  can 
improve  the  performance. 

6.  Adaptive  control  theory  usually  does  not  deal  with  problems 
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of  structural  Identification  (e.g.  failure  detection)  and 
structural  adaptation  after  a  failure  has  been  detected. 

However  problems  are  related:  Fast  structural  identification 
may  lead  to  false  alarms,  in  particular  under  noisy  conditions. 
Slow  and  reliable  structural  identification  may  leave  the 
system  in  a  failed  unstable  configuration  for  a  while.  The 
control  should  be  designed  to  provide  robustness  of  stability 
with  respect  to  the  failure  to  ensure  that  nothing  very  bad 
happens  until  the  failure  is  detected  reliably. 

7.  Robust  fixed  gain  control  may  be  combined  with  some  redundancy 
concepts.  Various  levels  are  possible: 

a)  Passive  redundancy  be  paralleled  components.  For  example 
the  50%  gain  reduction  margin  of  LQ  designs  offers  the 
possibility  of  using  two  paralleled  sensors  or  actuators 
such  that  in  case  of  a  failure  the  gain  is  reduced  only 
by  50%. 

b)  Removal  of  failed  components.  Even  if  a  component  failure 
can  be  tolerated,  as  far  as  stability  is  concerned,  it 
may  be  necessary  in  the  long  run  to  remove  a  failed 
component,  e.g.  to  close  a  leaking  gas  jet  valve  or  to 
remove  a  bias  term  entering  into  a  control  system  from  a 
sensor  failed  at  a  nonzero  constant  value. 

c)  Analytic  redundancy  may  help,  if  an  adaptive  observer 
provides  an  estimate  for  a  missing  signal. 

d)  Hardware  redundancy,  e.g.  majority  voting  in  a  multiplexed 
system  can  bring  the  system  back  to  its  original  performance. 
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However  this  part  of  the  system  ideally  should  not  be 
vital  for  stability  (see  6). 

With  the  availability  of  cheap  computers  there  are  few  constraints  on  the 
complexity  of  the  controller  structure.  However  there  are  several  practical 
limitations  and  additional  aspects  for  the  design.  Some  of  them  are  given  here. 

1.  A  main  constraint  on  the  controller  complexity  is  given  by  the 
presently  available  design  methods.  Also  it  is  a  question  with 
which  methods  the  designer  has  experience  and  for  which  he 

has  design  software  available.  Ideally,  control  theory  should 
provide  the  designer  with  convenient  tools,  e.g.  for  the 
computer-aided  design  of  control  systems,  instead  of  demanding 
that  the  designer  has  to  put  all  thinkable  tradeoff  situations 
into  one  scalar  performance  index  or  set  of  inequalities. 

2.  In  many  control  problems  structural  limitations  are  mainly 
dictated  by  the  cost,  availability  and  reliability  of  sensors 
and  actuators.  Thus  output  feedback  and  saturation  of  control 
or  control  rate  are  important  design  considerations.  For 
robust  control  systems  with  control  constraints  it  is  a 
particularly  important  rule  of  thumb  to  make  only  physically 
reasonable  requirements.  For  example,  one  should  not  try  to 
make  a  slow  system  fast  or  a  fast  system  slow  (i.e.,  do  not 
use  one  reference  model,  but  fast  and  slow  reference  models 
for  different  operating  conditions,  or  demand  only  invariance 
of  damping  or  maximum  overshoot,  not  of  natural  frequency  or 
time  of  maximum  overshoot).  Of  course  such  considerations 


depend  on  which  property  of  the  application  must  be  robust. 

In  the  design  of  an  oscillator  the  frequency  must  be  constant 
whereas  in  the  design  of  a  crane  the  frequency  is  less 
important.  Also  it  would  be  not  reasonable  to  require  an 
unchanged  performance  in  cases  of  sensor  or  actuator  failures. 
Thus  emergency  specifications  should  be  given  for  these  cases. 

3.  Some  special  consideration  must  be  given  to  the  case  for 

which  there  is  a  man  (operator,  pilot)  in  the  outer  loop.  A 
pilot  does  not  want  to  be  a  passenger.  He  may  want  to 
identify  the  controlled  aircraft  by  "playing"  with  the  input 
signals.  Control  schemes  which  give  him  the  same  feeling  for 
a  wide  range  of  parameter  variations  may  be  dangerous  if  the 
dynamics  suddenly  become  bad  beyond  an  assumed  range  of 
parameter  variation.  The  pilot  needs  a  warning  before  the 
"cliff".  This  is  another  reason  why  the  dynamics  should 
change  with  changing  parameters. 

A  plant  operator  may  want  to  switch  one  or  more  loops  from  automatic  to  manual 
control.  It  is  desirable  that  he  always  sees  a  stable  system  or  better  one, 
which  he  can  easily  control.  A  man  can  control  an  unstable  plant  provided 
the  eigenvalues  in  the  right  half  plane  are  close  to  the  origin.  He  has  more 
problems  if  he  has  to  control  fast  modes,  even  if  they  are  slightly  damped. 

In  other  words,  the  imaginary  axis  is  not  necessarily  the  best  emergency 
boundary  for  sensor  failures. 

The  problems  of  actuator  and  sensor  failures  look  similar  if  we 
interpret  them  as  a  row  or  a  column  of  the  feedback  matrix  being  switched  to 


zero.  For  the  human  operator  these  cases  are  quite  different,  since  for  many 
variables  he  has  sensors,  which  he  can  use  as  backup.  This  is  rarely  the 
case  for  actuators. 
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CHAPTER  3 

A  SURVEY  OF  ROBUST  CONTROL  RESULTS 


3.1.  Introduction 

The  discussion  of  the  preceding  chapter  has  demonstrated  that 
robustness  is  an  Important  issue  in  areas  of  control  system  design  ranging 
from  fixed  gain  controllers  to  completely  adaptive  control  systems.  The 
intent  of  this  report  is  to  concentrate  on  fixed  gain,  linear  time  invariant 
control  systems.  Hereafter  in  this  report,  references  to  robust  control 
system  design  will  assume  this  structure. 

There  are  two  basic  philosophies  to  the  analysis  and  design  of 
robust  control  systems.  The  first  assumes  that  the  perturbations  are  largely 
unstructured.  The  objective  is  to  design  the  controller  for  the  worst  perturba¬ 
tion  and  evaluate  or  bound  the  size  of  the  permitted  perturbations. 

The  second  approach  assumes  that  the  structure  of  the  disturbances 
and  their  size  are  known  a  priori.  The  design  objectives  are  to  minimize  the 
sensitivity  of  the  closed  loop  systems  in  the  known  perturbations  directions. 

Each  of  the  two  approaches  have  advantages  and  disadvantages.  The 
first  is  more  likely  to  provide  robustness  with  respect  to  unmodeled  errors 
and  hence  requires  less  accurate  models.  The  results  are  often  very  conserva¬ 
tive.  The  second  requires  accurate  models  of  the  perturbations  which  one 
likely  to  occur.  However,  the  controL  effort  is  directed  where  it  is  needed 
most. 

The  purpose  of  this  chapter  is  to  review  the  work  which  has  been 
done  on  the  robust  controL  problem.  The  survey  is  divided  into  two  sections 
corresponding  to  the  approaches  outlined  above. 


Sensitivity  and  Unstructured  Perturbations 


■j  •  d-  • 

3.2.1.  Frequency  Domain  Methods 

The  main  reasons  for  the  use  of  feedback  are  stabilization  and  the 
preservation  of  desirable  system  properties  in  spite  of  noise  inputs  and 
perturbations  of  system  parameters.  The  reduction  of  nonlinear  distortions 
was  an  essential  reason  for  the  use  of  feedback  amplifiers,  (see  Black  £ 1J ) - 
The  reduction  of  nonlinearity  by  high  gain  feedback  has  been  further  investi¬ 
gated  by  Cruz  [2]  and  Desoer  and  Wang  [3], 

In  frequency  design  methods  the  concept  to  compensate  the  loop, 
such  that  high  gains  are  possible  without  instability,  is  the  classic  rule  of 
thumb  for  the  reduction  of  noise  and  uncertainty.  Bode  [4]  expressed  it  in 
terms  of  gain  and  phase  margins  and  a  sensitivity  function,  which  was 
generalized  to  the  multivariable  case  by  Cruz  and  Perkins  [5j.  A  sensitivity 
matrix  S(s)  relates  the  output  errors  Ec(s)  due  to  perturbations  in  a  feedback 
system  to  the  output  errors  Eq(s)  due  to  the  same  perturbations  in  a  correspond¬ 
ing  open  loop  system  by  Ec(s)  * S(s)Eq(s)  .  The  sensitivity  matrix  S(s)  is  the 
inverse  of  the  return  difference  matrix,  for  the  loop  of  Fig.  3.1. 

S(s)  *  [I  +  G(s)K(s)H(s)]'1  (3) 


Figure  3.1.  Feedback  system,  return  difference  for  loops 
broken  at  a. 
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Note  that  G(s)  is  the  actual  plant,  which  may  be  expressed  by  the  nominal 
design  model  G^s)  and  a  perturbation  6G(s),  i.e.  G(s)  =GN(s) +6G(s).  If 
the  known  G^(s)  is  used  in  eq.  (3)  instead  of  the  unknown  G(s),  then  all 
results  are  local,  i.e.  restricted  to  small  6G(s).  For  a  reduction  of 
sensitivity  it  is  sufficient  that 

S^(-ju))S(jiu)  -  I  ^  0  (neg.  semidefinite)  (4) 


over  the  frequency  band  of  interest,  or  in  terms  of  the  return  difference 
F  (s)  =  I  +G(s)K(s)H(s) 

FT(-jt»>)F(ju>)  -  I  iO.  (5) 

Hsu  and  Chen  [6]  proved  the  relationship 


det  F(s) 


closed  loop  characteristic  polynomial 
open  loop  characteristic  polynomial 


(6) 


Thus,  if  no  cancellations  occur,  closed  loop  stability  can  be  analyzed  using 
det  F(s).  MacFarlane  [7]  studied  the  eigenvalues  Hj(s),  j=l,2,...,m  of 
F(s)  and  showed  that  the  closed  loop  is  stable,  if  all  characteristic  fre¬ 
quency  loci  Pj(jw),  j  *1,2 . in  satisfy  the  Nyquist  criterion.  He  also 

proved  a  necessary  condition  for  the  system  to  be  optimal  in  the  sense  of  a 

CO 

T  T 

quadratic  criterion  j  (y  Qy  +  u  Ru)dt: 

0 

I  Pj  (j^)  I  a  1  for  0  £  - 50  j  =  1,2, . . .  ,m  (7) 


or 


det  F ( j«0  I  2  1 


for  all  ao. 


(8) 


These  results  have  the  graphical  interpretation  that  the  complex  plane  plots 
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of  Idet  F(jaj)l  or  I P  ^  ( j^)  I  must  not  penetrate  the  interior  of  the  unit  disc. 

It  follows  from  this  that  the  characteristic  frequency  loci  of  an  optimal 
proportional  feedback  controller  have  infinite  gain  margin  and  at  least  60° 
phase  margin. 

Robustness  of  stability  with  respect  to  gain  and  phase  changes  may 
also  be  achieved  in  design  by  Rosenbrock’s  inverse  Nyquist  array  [8],  Here 
I-rGo^(ju)>  with  Gq(s)  =  G(s)K(s)H (s)  (see  Fig.  2)  is  analyzed  graphically  and 
modified  in  the  design.  A  standard  technique  in  multivariable  control  system 
design  is  to  use  compensation  or  feedback  to  decouple  or  approximately  decouple 
a  multivariable  system  into  several  single  input  systems,  which  may  be  designed 
by  single-loop  techniques.  Rosenbrock  [ 8 J  uses  the  criterion  of  diagonal 
dominance  for  approximate  decoupling. 

Doyle  showed  by  counterexamples  [9]  that  these  methods  can  lead  to 
highly  optimistic  margins  for  individual  loop  gains,  even  if  only  very  small 
margins  exist  for  simultaneous  change  of  several  loop  gains.  Already  in  the 
single-input  case,  gain  and  phase  margins  are  insufficient  to  characterize 
what  happens  for  simultaneous  gain  and  phase  perturbations.  Another  difficulty 
is  that  by  compensation  or  feedback  for  diagonal  dominance  the  actual  loca¬ 
tion  of  the  uncertainty  is  obscured. 

Doyle  [9]  examines  the  properties  of  the  return  difference  using  the 
concepts  of  singular  values,  singular  vectors  and  the  spectral  norm  of  a 
matrix.  The  singular  value  of  a  matrix  A  are  the  non-negative  square  roots 
of  the  eigenvalues  of  A*A,  where  A*  is  the  conjugate  transpose  of  A.  Since 
A*A  is  Hermitian,  its  eigenvalues  are  real.  The  singular  values  give  a  measure 
of  how  close  A  is  to  being  singular.  The  ratio  of  the  smallest  singular 
value  £  and  the  largest  one,  c,  is  the  condition  number  £/a.  One  may  also 
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interpret  the  singular  values  as  generalizing  to  matrices  the  notion  of  gain. 
This  characterization  is  of  great  practical  value,  since  good  software  to 
compute  a  singular  values  is  widely  accessible  [10].  Using  this  singular 
value  concept  Doyle  proved  the  following  robustness  theorem. 


r 

I  +  L(s) 

A  (s) 

y 

i 

;  i 
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— ► 

FP-7 041 

Figure  3.2.  Perturbation  by  L(s). 


In  the  system  of  Fig.  3.2,  let  G(s)  be  rational,  square,  invertible  and  such 
that  the  nominal  closed  loop  with  L(s)  *0  is  stable,  i.e.  G(I  +  G)  ^  =  I  +  G  ^  is 
stable.  If  the  system  is  perturbed  by  L(s),  which  by  itself  is  stable,  then  the 
perturbed  system  is  stable  if 

£(I+G  >c(L(juu))  for  all  a).  (9) 

For  this  theorem  Sandell  [11]  gave  a  different  proof,  in  which  G(s)  need  not 
be  rational.  £(I+G  ^(j^))  is  a  frequency  dependent  measure  of  robustness  in 
terms  of  gain  margins.  For  the  eigenvalues  X  of  A  (here  =  I+G  ^(j^))  generally 
the  relation 

£(A)sk(A)U:(A)  (10) 

holds.  It  is  possible  that  the  smallest  eigenvalue  is  much  larger  than  £(A). 
Thus  the  minimum  singular  value  £  gives  a  more  reliable  measure  of  robustness 
than  the  smallest  eigenvalue.  In  fact  Doyle  constructed  an  example,  where 
the  diagonal  dominance  approach  as  well  as  the  characteristic  loci  approac'.i 
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generates  a  Nyquist  or  Inverse  Nyquist  plot,  which  shows  +  «  db  gain  margin 
and  90°  phase  margin.  However  the  system  is  only  marginally  stable. 

The  problem  of  uncertainties  due  to  a  reduced  order  design  model  is 
interrelated  with  the  question  of  which  modes  of  the  system  must  be  influenced 
by  the  control  and  which  others  should  ideally  not  be  influenced  at  all.  In 
vehicle  control  it  may  for  example  be  desirable  to  control  the  rigid  body 
dynamics  fast  and  accurately,  i.e.  with  a  reasonably  high  bandwidth,  without 
interf erring  with  structural  vibrations.  In  frequency  domain  design  techni¬ 
ques,  this  is  achieved  by  a  40  db/decade  roll  off  beyond  the  design  band¬ 
width.  This  aspect  is  frequently  ignored  in  state  space  design  techniques. 

In  all  design  techniques  it  is  important  to  study  carefully  the  behavior  in  a 
frequency  range  above  the  bandwidth,  where  modes  are  still  sufficiently 
controllable  and  observable,  such  that  the  control  may  move  them  into  the 
right  half  s  plane. 

Stein  and  Doyle  [12]  give  a  design  example  for  a  CH-47  helicopter 
with  two  control  inputs.  They  apply  singular  value  analysis  and  the  robustness 
condition  (9).  Rotor  dynamics  and  rate  limits  are  translated  into  a(L(juu)) 
using  a  result  of  Safonov  [13].  The  two  singular  values  were  made  approxi¬ 
mately  equal  and  the  bandwidth  in  both  loops  was  increased  as  much  as  G(L(j'x)) 
admitted.  A  lev  pass  helped  to  meet  the  "roll-off"  requirement.  The  example 
also  showed  that  these  methods  nay  lead  to  very  conservative  results  in  cases 
of  large  variations  of  parameters  in  specific  directions,  here  the  flight 
condition  variation. 

3.2.2.  State  Space  Methods 

Single-input  linear  quadratic  state  feedback  regulators  have  a 
return  difference  greater  than  unity  at  all  frequencies,  a?  was  shown  by  Kalman 
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[14].  Anderson  and  Moore  [15]  showed  that  this  fact  implies  a +60°  phase 
margin,  infinite  gain  margin  and  50  percent  gain  reduction  tolerance.  Safonov 
and  Athans  [16]  generalized  this  result  to  the  multiinput  case: 

x  =  Ax  +  Bu 

u  =  -Kx  (11) 

with  m  inputs  u^. 

The  feedback  matrix  K  is  determined  by  solving  a  Riccati  equation  minimizing 

00 

J  =  j  (x^Qx  +  u^Ru)dt  (12) 

0 

with  Q  positive  definite  and  R  =  diag[r^, . . . ,rj ,  r^  >0. 

The  individual  inputs  u^  are  perturbed  to  without  interaction 

between  them,  i.e. 


x  =  Ax  +  B?cu  with  TJu  =  I  . 


Let  each  perturbation  71.  be  linear  time  invariant  with  proper  rational  stable 

1  Jy^O") 

transfer  function  P^s).  Its  frequency  response  is  P^(juj)  =  a^uu)  »e 

Then  the  closed  loop  remains  stable  under  a  phase  perturbation  yL  (u>),  with 

£  60°  for  all  iw.  It  also  remains  stable  under  a  gain  perturbation 

a^Cuu)  2  0.5  for  all  w. 

Note  that  this  emphasizes  the  importance  of  the  bandwidth  of  the 
control  system.  The  60°  phase  margin  without  bandwidth  limitations  is  not 
sufficient  to  accomodate  neglected  error  dynamics  since  physical  actuators 
have  at  least  90°  phase  laR  at  high  frequencies.  For  this  reason  Otto  Smith 
[17]  used  the  "complex  gain  margin",  i.e.  the  minimum  distance  of  G(j^)  to 
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the  critical  point,  scaled  by  the  local  frequency  increment  along  G(jtu) . 

This  approximates  the  negative  real  part  of  a  dominant  pair  of  eigenvalues. 

A  multivariable  measure  for  the  distance  of  G(jtu)  from  the  critical  point  has 
been  discussed  already  in  form  of  the  singular  values  of  the  return  difference. 

Doyle  [18]  showed  by  counterexample  that  the  margins  may  be 
arbitrary  small  if  the  state  is  replaced  by  a  state  estimate  from  an  observer 
or  a  Kalman  filter.  In  his  example,  the  gain  margins  were  arbitrarily  small 
in  both  the  positive  and  negative  db  direction.  To  improve  the  margin  in  this 
situation,  Doyle  and  Stein  [19]  developed  a  "design  adjustment  procedure", 
which  introduces  fictituous  noise  at  the  control  input  to  the  plant.  In  this 
procedure  the  observer  eigenvalues  tend  to  the  finite  transmission  zeros  and 
to  infinity.  Thus  the  procedure  works  only  for  minimum-phase  plants.  The 
procedure  is  essentially  the  dual  of  Kwakemaak's  sensitivity  recovery  method 
[20].  This  however  drives  the  plant  poles  instead  of  the  observer  poles  to 
the  transmission  zeros,  which  may  lead  to  large  control  inputs  u. 

Gain  and  phase  margins  may  be  much  smaller  in  discrete  time  linear 
quadratic  state  feedback  systems.  Jacques  Willems  and  van  de  Voorde  [21] 
give  bounds  for  the  single-input  case,  which  show  that  the  system  may  be  very 
sensitive  to  feedback  gain  variations.  This  is  not  surprising,  since  the 
hold  element  may  be  approximated  by  a  phase  shift  of  one  half  sampling 
interval. 

Safonov  and  Athans  [16]  also  generalize  a  single-input  result  by 
Anderson  and  Moore  [15],  which  is  useful  for  actuator  nonlinearities.  If 
the  perturbation  operator  in  eq.  (13)  describes  a  time  varying,  memoryless 
nonlinearity  Tt^u^  =  f  ^  (u,  t)  ,  then  it  is  a  sufficient  condition  for  the  closed 


loop  stability,  that 
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j  ^  f(u,t)  £  M  for  some  M  <  “  and  for  all  t.  (14) 

For  example, for  an  actuator  saturation,  stability  is  guaranteed  if  the  inputs 
do  not  exceed  twice  the  saturation  level. 

Comparisons  of  numerous  optimization  techniques  for  insensitive 
control  systems  were  made  by  Harvey  and  Pope  [22,23]  for  wing  load  alleviation 
for  the  C-5A  aircraft  and  by  Vinkler  and  Wood  [24]  for  a  lateral  autopilot 
for  a  rudderless  remotely  piloted  vehicle.  A  minimax  technique  by  Salmon  [25] 
and  an  uncertainty  weighting  technique  by  Porter  [22]  were  judged  superior  to 
six  other  techniques  in  the  first  report;  both  however  failed  in  the  comparison 

[27] .  Here  an  expected  cost  technique  by  Ly  and  Cannon  [26]  and  a  multistep 
guaranteed  cost  technique  by  Vinkler  and  Wood  [27]  came  out  better  than  four 
other  techniques.  In  [23]  an  information  matrix  approach  by  Kleimann  and  Rao 

[28]  compared  favorably  with  other  techniques. 

In  problems  with  insignificant  constraints  on  the  control  inputs, 
the  weighting  matrix  R  in  a  quadratic  criterion  may  be  small.  This  leads  to 
high  gain  solutions  as  they  were  discussed  in  the  previous  section.  A 
comparison  of  various  high  gain  feedback  systems  is  made  by  Young,  Kokotovic 
and  Utkin  [29].  This  comparison  also  includes  variable  structure  systems, 
which  in  their  sliding  mode  are  insensitive  to  parameter  variations  and  distur¬ 
bances,  similar  to  the  high-gain  system  [30],  Young  [31]  applied  this  concept 
to  the  design  of  an  adaptive  model  following  control  system  and  compared  the 
results  for  the  longitudinal  motion  of  a  Convair  C-131B  aircraft  with  other 
model  following  techniques. 

A  special  case  of  a  high  gain  control  system  is  useful,  if  the 
reference  or  disturbance  input  signals  can  be  exactly  modeled  and  asymptotic 
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tracking  or  disturbance  rejection  is  required.  The  use  of  integrators  in 
the  loop  for  zero  stationary  errors  in  step  and  ramp  responses  is  a  classical 
recipe.  Also  for  other  inputs  an  internal  model  of  the  input  can  be  used, 
e.g.  a  tuned  oscillator  (notch  filter)  for  disturbance  rejection  of  helicopter 
rotor  vibrations,  whose  frequency  is  regulated.  Such  a  high  gain  at  particular 
frequencies  makes  asymptotic  tracking  robust  to  plant  parameter  variations 
as  long  as  the  loop  remains  stable.  This  robustness  problem  was  studied  by 
Davison  [32]  and  others.  In  sampled-data  systems  the  internal  model  is  to  be 
implemented  in  continuous  time,  if  the  tracking  property  is  required  also 
between  the  sampling  instants  [33]. 

Some  common  problems  in  all  high  gain  concepts  are 

•  Measurement  noise  goes  highly  amplified  to  the  actuator  inputs. 

•  High  values  for  1 u I  and  1*1  may  occur. 

•  Non-cooperative  efforts  of  the  actuators  may  occur. 

The  LQG  design  method  offers  a  systematic  way  to  avoid  these  difficulties  by 
increase  in  the  R  matrix  and  by  the  use  of  a  Kalman  filter. 

3,3.  Robustness  with  Respect  to  Large  Perturbations  in  Known  Directions 
3.3.1.  Parameter  Methods 

In  the  methods  of  Section  3  relatively  little  knowledge  about  the 
parametric  perturbation  is  assumed.  The  results  are  therefore  primarily  valid 
for  small  perturbations.  In  some  cases  information  is  obtained  about  how  big 
the  perturbation  is  allowed  to  be  in  order  to  maintain  stability. 

In  situations  where  large  perturbations  in  known  directions  occur, 
the  previous  methods  generally  lead  to  very  conservative  results.  In  this 
section  some  tools  are  discussed  by  which  such  perturbations  can  be  accomodated 
in  the  design.  The  next  chapter  and  [34]  also  describe  a  parameter  space 
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meChod  which  is  applicable  to  this  problem. 

In  typical  design  examples  not  only  the  mathematical  model  of  the 
plant  is  uncertain,  but  also  the  formulation  and  relative  weight  of  many  design 
criteria.  Some  of  these  criteria  are  in  form  of  inequality  constraints;  others 
are  to  be  minimized.  It  is  artificial  to  put  all  of  them  together  into  one 
scalar  performance  index,  which  is  then  minimized  over  the  parameters  in  an 
assumed  controller  structure.  For  the  designer  an  interactive  computer-aided 
design  procedure  is  more  useful,  where  he  can  make  higher  level  decisions  of 
how  to  change  requirements  after  each  computer  solution  or  failure  to  find 
a  solution.  Hie  computer  may  have  to  solve  a  nonlinear  programming  problem  in 
each  design  step.  Various  aerospace  problems  have  been  formulated  and  solved 
this  way.  Schy  [37,38]  deals  with  a  lateral  stability  augmentation  system  for 
a  fighter  airplane,  Hauser  [39]  with  an  autopilot  for  a  flexible  space  vehicle. 
Further  design  examples  are  given  by  Karmarkar  [40]  and  Kanarachos  [41],  It 
is  convenient  to  formulate  all  design  criteria  for  each  operating  point  as 
components  of  a  performance  vector  £.  It  may,  for  example,  contain 

•  bounds  on  the  individual  feedback  gains  |k„|. 
and  for  each  flight  condition  specifications  on 

.  eigenvalue  location. 

•  deviation  from  nominal  response  for  typical  reference  and 
disturbance  inputs. 

.  bounds  on  the  control  rate  I  ill  for  typical  reference  and 
disturbance  inputs. 

Kreisse lmeier  and  Steinhauser  [42]  use  in  an  example  with  five  flight  conditions 
of  a  F4-C  aircraft  a  40  dimensional  vector  £.  A  vector  constraints  <c  (i.e. 
componentwise  £ c is  given  and  the  feedback  gains  K  are  the  solution  of 
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the  problem 

Min (Max  g. (K)/c  }.  (17) 

K  i 

Using  an  algorithm  described  in  [43]  Kreisselmeier  and  Steinhauser  obtain  a 
Pa re to -optimal  solution.  Figure  3.3  shows  some  reference  step  responses  of 
this  design  for  an  F4-C.  It  is  stable  in  the  five  flight  conditions.  The 
open  loop  responses  on  the  left  side  show  that  the  aircraft  is  slow  in  flight 
condition  1  (landing  approach).  Here  a  slower  reference  response  was  given 
than  for  the  high  speed  condition  2  and  4.  The  desired  reference  response 
was  specified  as  gi(t)  =  gM(«it)  where  for  each  flight  condition  i»l,2,...,5 
an  appropriate  time  scale  was  chosen.  This  resulted  in  the  insensitive 
closed  loop  responses  on  the  right  side  of  Fig.  3.3,  which  required  only  a 
relatively  small  control  rate  Idl.  The  same  feedback  resulted  in  similarly 
good  disturbance  responses. 

Also  the  results  of  Shy  [38]  showed  that  an  amazingly  large  variation 
of  parameters  can  be  accomodated  by  a  fixed  gain  controller,  if  the  requirements 
were  in  good  agreement  with  the  physical  limitations.  These  designs  result 
in  low  gain  solutions,  and  the  dynamics  change  in  an  acceptable  or  desirable 
way  as  the  physical  parameters  vary. 

3.3.2,  Integrity:  Robustness  with  Respect  to  Sensor  and  Actuator  Failures 
If  an  actuator  or  sensor  is  connected  to  a  high  gain,  then  its 
failure  is  a  larger  perturbation  than  in  a  low  gain  situation.  Thus  requirements 
for  robustness  with  respect  to  actuator  and  sensor  failures  tend  to  result  in 
low  gain  solutions.  Even  more  important  Is  the  aspect  of  avoiding  non- 
cooperative  efforts  of  actuators.  If,  for  example,  one  input  alone  places  some 
eigenvalues  in  the  right  half  plane  and  another  is  needed  to  bring  them  back 


uncontrolled 


insensitive  control 


Figure  3.3.  Response  Co  a  step  command:  F-4  (Phantom)  aircraft  at  5  extremal 
flight  conditions  (altitude  0  ...  40  000  ft,  Mach  number  .2  ... 
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into  the  left  half  plane,  then  apparently  no  robustness  of  stability  with 
respect  to  actuator  failures  can  be  achieved. 

One  approach  to  achieve  robustness  of  stability  with  respect  to 
certain  failures  is  to  try  to  extend  gain  reduction  margins  to  include  gain 
zero.  Belletrutti  and  MacFarlane  [44]  use  the  term  "high  integrity"  for 
robustness  with  respect  to  certain  failures.  They  check  the  stability  condi¬ 
tions  for  gains  reduced  to  a  small  e  using  Nyquist  stability  criteria  for 
characteristic  loci  of  principal  submatrices  of  the  return  ratio.  In  this 
analysis  the  loop  must  be  broken  at  the  point  where  the  actual  failure  may 
occur  and  thus  the  gain  reduction  margin  is  needed.  Owens  [45]  derived 
necessary  and  sufficient  conditions  for  integrity  of  systems  with  multivariable 
proportional-integral  controllers . 

Solheim  [46]  formulated  the  integrity  problem  in  the  context  of 
quadratic  optimal  control.  In  examples  an  increased  integrity  is  obtained  with 
an  increased  weight  R  on  the  control  in  the  quadratic  criterion,  another 
indication  that  the  solution  will  tend  to  a  low  gain  solution.  Wong,  Stein 
and  Athans  [47]  show  the  following  gain  reduction  result  for  LQ  regulators: 

The  matrix  A  (A)  =  A  +  BAK  with  A  =  diag[a, . .  .at  ],  where  K  minimizes 
c  1  1  m  ’ 

CO 

x'Qx  +  u'Rxdt  for  A  =  I,  is  stable  for  all 

A  >  j[I  -  (R1/2K'Q'1KRL/2)"L] .  (18) 

This  generalizes  the  bound  or^>0.5  from  [16].  The  recommendation  is,  from  a 
purely  robustness  standpoint,  to  choose  Q  and  R  such  as  to  maximize 

\n  =  X  .  l(R1/2K,Q_1KIl1/2)-1}.  (19) 

0  nun  7 

Kreisselmeier  [48]  proposes  to  modify  the  quadratic  criterion,  where  for  each 
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considered  failure  situation,  a  quadratic  criterion  is  formulated  and  the 
overall  criterion  is  a  weighted  sum  of  these  terms. 

In  failure  situations  it  may  be  desirable  to  specify  other  emergency 
boundaries  in  the  eigenvalue  plane  than  only  the  imaginary  axis.  This  problem 
is  treated  by  parameter  space  methods  in  [34J  and  in  Chapter  4.  The  concept 
is  illustrated  for  the  case  of  sensor  failures  in  Fig.  3.4.  A  nominal  region 
for  the  eigenvalue  location  and  a  larger  emergency  region  are  mapped  into 
the  space  of 


Figure  3.4.  Illustration  of  failure  robustness  and 
emergency  boundaries. 

feedback  gains.  It  is  assumed  that  the  syscera  is  represented  in  "sensor 
coordinates",  then  a  failure  of  a  sensor  for  state  variable  corresponds  to 
switching  k_^  to  zero.  The  projection  of  point  1  on  the  k^  axis  is  outside 
the  emergency  boundary,  i.e.  the  emergency  specification  is  not  robust  with 
respect  to  a  sensor  failure  k?  =0.  It  is,  however,  robust  with  respect  to 
k^  =  0.  For  all  points  in  the  shaded  area  the  emergency  specifications  are 
robust  with  respect  to  either  sensor  failure.  An  alternative  to  this  robust 


27 


solution  would  be  in  this  example  to  omit  sensor  1  and  to  use  multiplexed 
sensors  for  x^  and  failure  detection. 

In  the  multiinput  case  a  sensor  failure  is  equivalent  to  changing  a 
column  of  the  K  matrix  to  zero  and  an  actuator  failure  is  equivalent  to 
changing  a  row  of  K  to  zero.  In  [34]  an  actuator  failure  example  is  studied, 
where  the  problem  is  formulated  such  that  the  eigenvalues  are  placed  in  a 
nominal  position  with  two  actuators  and  move  as  little  as  possible  towards 
the  stability  boundary  for  failures  of  either  one  of  two  actuators. 

Apparently  a  necessary  condition  for  robustness  with  respect  to 
failures  is  that  the  insufficiently  damped  eigenvalues  (outside  the  specified 
region)  remain  controllable  and  observable  after  the  failure.  In  the  crane 
example,  the  sensor  for  the  crab  position  x^  is  essential,  because  x^  is  not 
observable  by  other  states.  In  such  situations  it  is  apparently  misleading  to 
use  high  gain  feedback  and  to  show  gain  reduction  to  only  a  few  percent  of  the 
high  gain.  For  failures  of  essential  actuators  and  sensors  only  redundant 
components  can  help. 
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CHAPTER  4 

PARAMETER  SPACE  DESIGN  OF  ROBUST  OONTROL  SYSTEMS 

4.1.  Introduction 

In  this  chapter  a  new  tool  for  the  design  of  robust  control  systems 
is  proposed.  First  the  type  of  robustness  problem  for  which  the  tool  can  be 
applied  is  described. 

Robustness  of  a  control  system  is  defined  in  terms  of  system  properties 
which  are  invariant  under  a  specified  class  of  perturbations.  The  system 
property  considered  in  this  paper  is  "nice  stability"  as  specified  by  a  region 
F  in  the  eigenvalue  plane,  in  which  all  eigenvalues  must  remain  in  spite  of 
perturbations.  The  perturbations  may  be  large  changes  of  physical  parameters 
of  the  plant  or  failures  of  actuators  and  sensors  or  inaccurate  implementation 
of  the  control  law. 

The  following  assumptions  are  made. 

1.  Only  linear  plants 

x(k+l)  =  A  x(k)  +  B  u(k)  or 

x(t)  =A  x(t)  +  B  u(t)  (20) 

X*  **  [Xj_  . . .  xnJ ,  U*  =  [uL  . . .  u  ] 

are  considered.  It  is  assumed  that  eq.  (20)  is  written  in 
"sensor  coordinates",  i.e.  all  measured  variables  are  state 
variables  x^.  It  may  be  part  of  the  design  task  to  decide 
which  states  are  to  be  measured. 

A  and  B  may  depend  on  a  physical  parameter  vector  Only 
some  typical  values 


2. 
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A.  =  4(6..),  Bj  =6(0^  j=l,2...J 


may  be  given.  Also  the  required  system  property  I\  may 
depend  on  the  operating  point  _0 ^ . 

3.  The  simplest  assumed  controller  structure  is  state  feedback 


u  =  -  K  x .  (22  ) 

It  may  not  be  possible  to  make  all  plant  models  A  ^  ,  B^  nicely 

stable  with  the  same  fixed  K,  i.e.  to  have  all  roots  of 

J 

rr  det(-\l -A.  +  B.K)  *  0  (23) 

j=l  3  3 

in  the  specified  region  T  in  the  ?i-plane.  In  this  situation 

the  designer  may  decide 

a)  To  relax  the  eigenvalue  region  specification  such  that  a 
solution  exists. 

b)  To  use  several  feedback  gains,  each  of  which  nicely  stabilizes 
a  group  of  pairs  A_.,B^.  The  gains  can  then  be  scheduled 
depending  on  a  measurement,  which  admits  a  distinction 
between  the  groups. 

c)  To  use  linear  dynamic  feedback  of  order  m  with  controller 
state  vector  x  ,  i.e,  try  to  find  a  state  feedback 


u 

*ji 

i 

=.  - 

1 

u 

— p 

M  N  1 

_  _ i 

which  nicely  stabilizes  the  augmented  system 


The  controller  M,N,L  with  n  inputs  and  p  outputs  may  be 

written  in  a  canonical  basis,  such  that  the  feedback  matrix 

in  eq.  (24)  contains  pn-Hnn-Hnp  design  parameters.  ((p-Kn)  x  (n4m) 

2 

coefficients  in  eq.  (24)  of  which  m  are  normalized  by  the 

choice  of  an  mxm  transformation  matrix). 

d)  To  use  nonlinear  feedback,  e.g.  an  adaptive  system 

estimating  the  physical  parameter  vector  _6  in  eq.  (21). 

The  tool  proposed  in  this  paper  is  useful  for  problems  of  the  types  a),  b), 
and  c ) . 

If  some  states  are  not  available  for  feedback,  then  the  corresponding 

columns  of  K  in  eq.  (22)  or  of  K  and  M  in  eq.  (24)  are  zero. 

4.  A  sensor  failure  is  equivalent  to  switching  all  elements  of  the 

corresponding  column  of  K  (or  K  and  M)  to  zero.  An  actuator 

failure  is  equivalent  to  switching  all  elements  of  the 

corresponding  row  of  K  (or  K  and  L)  to  zero.  Assume  that  M 

such  failures  or  failure  combinations  have  to  be  considered, 

which  lead  to  M  crippled  feedback  matrices  K^,  m = 1,2 , . . . ,M. 

The  design  goal  is  to  find  K  such  that  all  roots  of 
M  J 

rT1  ,rr.  (?>l-A.+B  .K  )  are  inside  an  "emergency  region"  in  \-plane. 
m=t  j=i  ~  j  “j- ra 

Apparently  it  is  a  necessary  condition  for  nice  stability  to 
be  robust  under  such  failures,  that  the  plant  modes  outside 
the  specified  region  F  in  the  eigenvalue  plane  remain  control¬ 
lable  and  observable  after  the  failure.  This  fact  may  be  used 


in  the  decision  for  which  sensors  and  actuators  redundancy 


is  necessary.  Components  which  are  needed  to  make  the 
insufficiently  damped  modes  observable  and  controllable  are 
"essential1’. 


Example  4.1:  (Essential  sensor)  Consider  a  crane  with  the  physical  parameters 
mc  = crab  mass,  = load  mass,  i  = rope  length,  g  = gravitational  constant.  Its 
state  variables  are  x^ = crab  position,  = crab  velocity,  x^ = rope  angle,  and 
x^  =  rope  angular  velocity.  For  small  rope  angles  the  linearized  state 
equations  are 
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(2  6) 


with  (Up =  (mc+m^)g/mc2.  Input  u  is  the  force  accelerating  the  crab.  Eigen¬ 
values  are  [0,0,  -jx^] .  The  observability  analysis  shows  that  x^  is  not 

observable  by  x^t  x^>  or  x^'  Since  one  of  the  zero  eigenvalues  is  unobservable, 
the  crab  position  sensor  is  essential  for  stabilization. 


5.  It  is  assumed  that  desirable  features  of  the  dynamic  behavior 
of  the  control  system  can  be  specified  by  a  region  F  in  the 
eigenvalue  plane  .  Examples  are 

a)  the  stability  region,  i.e.  the  left  half  s  plane  or  the 
interior  of  the  unit  circle  in  the  z  plane, 

b)  military  specifications  for  damping  and  natural  frequency  of 
modes  of  an  aircraft, 

c)  in  some  problem  formulations  it  is  convenient  to  define  a 


family  of  regions  with  a  parameter  r.  In  the  z-plane 
this  may  be  a  family  of  nonintersecting  circles  (see  Fig. 
4.1).  For  r  =  0  it  is  a  deadbeat  solution,  and  r  in  the  range 
0.3  to  0.5  corresponds  to  well  damped  transients.  (The 
right  shift  of  this  circle  excludes  heavily  oscillatory 
solutions,  the  circles  approximate  the  usual  logarithmic 
spirals  for  constant  damping  augmented  by  a  constraint  on 
Izl).  For  r  =  l  the  stability  boundary  is  obtained. 

Similarly  in  the  left  half  s-plane,  a  family  of  hyperbolas 
(guaranteeing  a  minimum  damping  and  a  minimum  negative  real 
part  of  the  eigenvalues)  may  be  introduced,  augmented  by 
parallels  to  the  imaginary  axis  in  the  right  half  plane 
(see  Fig.  4.2).  The  equations  for  these  families  of 
boundaries  will  be  given  later. 

Three  types  of  design  problems  for  robust  control  system  may  be  distinguished. 

1.  Given  a  system  property  and  a  controller  initially  designed  for 
nominal  parameter  values,  under  what  perturbations  of  the  para¬ 
meters  is  the  property  robust?  Modify  the  controller  such  as 
to  extend  the  admissible  class  of  perturbations.  Example: 

Try  to  extend  the  gain  and  phase  variations  under  which 
stability  is  robust.  It  is  difficult  in  this  approach  to 
accomodate  large  perturbation  in  known  directions,  e.g. 
large  variations  of  the  physical  parameter  vector  2.  in  ea.  21. 
Typically  very  conservative  results  are  obtained. 

2.  Given  a  class  of  perturbations  and  a  system  property  with  a 


parameter  r,  e.g.  the  family  of  boundaries  f  ,  for  a  given 
controller  which  is  the  best  value  of  r  such  that  1’  is  robust 
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under  all  perturbations  in  the  given  class?  Modify  the 
controller  to  improve  r.  The  system  property  may  be  formulat¬ 
ed  such  as  to  include  several  design  aspects  and  an  optimization 
can  be  performed. 

3.  Given  a  system  property  and  a  class  of  perturbations,  does 
there  exist  a  state  feedback  solution  such  that  the  property 
is  robust?  If  not,  does  there  exist  a  linear  dynamic  feedback 
controller  of  order  m,  eqs.  (24)  and  (25),  such  that  the 
property  is  robust?  Find  the  set  of  admissible  controllers 
and  select  one  based  on  criteria,  other  than  robustness,  e.g. 
based  on  simulations  with  a  nonsimplif ied  nonlinear  plant  model. 

It  is  primarily  this  third  problem  formulation  for  which  the  proposed  design 
cool  can  be  applied.  The  concepts  however,  are  useful  also  in  the  second 
problem.  The  design  tool  basically  consists  in  mapping  boundaries  from  the 
\  plane  into  boundaries  in  the  parameter  space  X,  whose  coordinates  are  the 
elements  of  the  state  feedback  matrix.  Then  all  tradeoffs  between  dynamic 
requirements,  robustness  requirements,  and  bounds  on  the  feedback  gains  can 
be  made  in  X  space. 

The  details  and  most  of  the  examples  are  worked  out  for  single  input 
plants  with  state  feedback 

u  =  -k,i  =  -[k1  k2  ...  kn]x  (27) 

or  state  output  feedback,  i.e.  some  k^ = 0.  For  multi-input  plants  the  basic 
result  is  formulated. 

Parameter  space  methods  have  a  long  tradition,  mainly  in  Russia 
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and  Yugoslavia.  Siljak  [49]  gives  a  historical  review  of  the  work  by 
Vishnegradsky,  Neimark,  Mitrovic,  and  others.  Siljak  generalized  these 
parameter  mapping  methods  significantly.  A  typical  procedure  for  a  continuous 
time  system  is  to  assume  a  controller  structure  with  two  free  parameters  c* 
and  3.  Determine  the  closed-loop  characteristic  polynomial 

n  i 

P(s)  =  =  0.  (28) 

Substitute  s=a  +  ju)  and  separate  eq.  (28)  into  its  real  and  imaginary  parts: 
Re(a,(u,a,P)  =  0,  Im(a,<ju,a,  3)  =  0.  Assume  these  nonlinear  equations  have  a 
solution 

a  ~  a(a,u;) ,  @  =  3(ct,uj).  (29) 

Equation  (29)  allows  mapping  c,u>  pairs  on  the  boundary  into  the  or-3-plane. 

The  image  boundaries  divide  the  a-3-plane  into  regions  characterized  by  the 

number  of  eigenvalues  inside  and  outside  the  s-plane  region. 

In  the  present  paper  the  control  system  structure  is  restricted  to 

partial  state  feedback.  This  permits  simplifying  the  determination  of  eq.  (29) 

by  pole  placement  methods.  Consider  for  example  a  second  order  single-input 

system  with  k^  =  a;,  ^=3  in  eq.  (27).  In  classical  parameter  plane  methods 

2 

P(s)  *  det(s_l-A+b  k' )  =  pQ  (a,  8)  +  p^(a,  3)s  +  s  =  0  is  determined  and  with  s=a+juj 
solved  for  ce  and  3.  In  the  method  proposed  in  this  paper  the  p^  are  expressed 
in  terms  of  a  and  tu  by 

2  2  2 
P(s)  =  (s-a+joj)  (s-c- ju;)  =  s  -2?s+g  +'ju 

7  (30) 

=  PQ  (<?/■»)  +P1(a)s  +  s“  =  0. 


Then  by  pole  placement 
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kL  =a(po,pL)  =a(CT,tu) 

k2  “P<po»pl>  =P(a>U))- 


(31) 


Thus  the  mapping  equation  (29)  is  obtained  in  a  different  way. 

More  generally  for  an  nth  order  single  input  system  in  both 
approaches  an  n  dimensional  parameter  space  &  with  coordinates  p^  is  introduced 
as  an  intermediate  step  between  the  set  of  eigenvalues  A={X^  . . .  and  the 
^-space.  The  relation  between  A  and  X  can  be  expressed  in  both  directions: 

a)  From  X  to  9  by  the  characteristic  equation  P(X)  =  det  (Xl-A+b  k*  ) , 
from  &  to  A  by  numerical  factorization  of  P(\). 

b)  From  A  to  @  by  multiplication  of  elementary  factors  P(\)  = 

•••  (^“^n)»  from  &  to  X  by  pole  placement. 

Apparently  direction  b)  is  much  simpler  than  direction  a).  In  this  paper  only 
direction  b)  is  used. 

In  the  next  section  pole  placement  is  reformulated  as  a  linear 
mapping  from  #  to  X  space.  This  is  then  used  in  the  third  section  to  map  not 
only  a  trial  point  from  &  to  X  or  from  A  to  X>  but  to  map  boundaries. 


4.2.  Single  Input  Pole  Plac  e.  ient--A  Linear  Mapping 
4.2,1.  State  Feedback 

Theorem  1:  Givea  a  polynomial  F(\)  =  p  +p,\  +  ...+p  ,\n  ^+A.n,  a r.  nxn  matrix 

-  o  i  n-i 

A  and  an  n  x  1  vector  b  such  that  det  R#0,  R=[b,A  jj  . . .  A°  ^b] .  The  unique 
solution  to  det (\I-A+b  k’ )  =  P(\)  is 

k'  =  £*E  (32) 


where 


P*  -  [PQ  PL  •••  Pn-1  lJ»  E  = 


e'A 


I  .u 

e'A 


.-1 


and  e'  is  the  last  row  of  R  . 

Proof:  Existence  and  uniqueness  of  the  solution  were  shown  by  Rissanen  [50] 
by  transformation  to  control  canonical  form. 

Let  F=A-b  k1  and  expand  powers  of  F  as  follows: 

(1) 


F°  =  A°  =  I 


F=A-bk' 


F2  =A2  -  Abk1  -  b  k'F 


(2) 

(3) 


j,n  ^  lfa  ^  An  k'F  bk'F^*^ 


(n+1) 


Multiply  the  first  equation  by  pQ,  the  second  by  p^,  etc.,  the  (n+l)st  row 
by  one  and  add  the  equations 


n-L 


P(F)  =  P(A)  -  [b,A  b  ...  A  b] 


k' 


By  Cayley-Hamilton  P(F)  =0.  Then 


=  R_1P(A),  R  =  [b,A  b  ...  A^b] 


(33 


k'  must  satisfy  the  last  row,  i.e. 


k'  =  e'P(A),  e'  -  [0  ...  0  1]R 


-1 


(34 


Explicitly 
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k‘  =  e’[p  I  +  p.,A+...+p  ,An  ^+An] 
—  —  or-  1—  n-I—  - 


k'  =  [Pc  Pi-.-Pn.!  U 


e ' 
j'A 

X 


n 


=  £*E. 


Q.E.D. 


The  form  (34)  of  the  result  was  derived  in  [33].  The  (n+1)  xn  matrix  E  is  a 

convenient  representation  of  a  controllable  pair  A,b.  It  maps  a  vector 

=  [PQ  p^  •  •  •  Pn  in  &  space  into  a  vector  k'  =£*E=  [jd1  1]E  in  X  space. 

E  is  evaluated  only  once  for  a  given  pair  A,b.  The  mapping  of  each  trial 

2  2 

design  point  in  £  space  then  requires  only  n  multiplications  and  n  additions. 
This  compares  favorably  with  mapping  a  trial  design  point  from  the  parameter 
space  of  quadratic  criteria  via  the  Riccati  equation  into  X.  space.  This  is 
an  advantage  for  computer-aided  design  methods,  in  which  many  trial  design 
points  have  to  be  mapped  and  displayed  graphically. 

Example  4.2:  (Foie  placement,  output  feedback,  gain  scheduling).  For  the 
crane  of  Example  4.1 


E  = 


2 

im  / g  0  i  m  /g 

c  c 

I  0  ItaJ  g  0 


2 

j e  mc/g 


o 

o 


0  -ini  0 

c 

0  0  -imc 


0 


0  (mc+mL)g  0 


kl  =  Pq^c78 

k2  =  Pl-V8 

2 

k3  =  p0£  mc/s  ‘  p2J&nc  +  •IVfmL)g 
k4  =  ^mc^pl'i7g'p3)* 


(35) 
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This  result  admits  some  conclusions  about  state  output  feedback  and  gain 
scheduling 


i) 

for  stability  pQ  >  0, 

i.e.  k^  >  0 

ii) 

for  stability  p^  >  0, 

i.e.  k2  >  0 

ill) 

k^ = 0  implies  the  constraint 

p2  (mc-HnL)g/imc> 

iv)  =  0  implies  the  constraint  p^  =  p^i/g, 
v)  a  gain  scheduleing  for  different  loads  m^  can  be  implemented  as 

k3  =  k30  +  mL8  Wlth  k30  =  Po£2mc/g_P2imc+mc8* 

The  other  k^  remain  unchanged.  With  this  gain  scheduleing  the 
closed  loop  eigenvalues  do  not  depend  on  the  load. 

in  numerical  calculations  with  large  n  the  accuracy  of  the  vectors  jj'A*-, 
1=1,2, ... ,n,  must  be  checked.  One  test  is  to  let  pQ  = p^ = . . .  = Pn_^ = 0* 

Then  k'  =  e'An.  Evaluate  det  (XI-A  +  b  e '  An)  =  pQ  +  p^X  +  . . .  +  p^_^  +  \n.  The  p^ 
should  ideally  be  zero.  Their  magnitude  is  a  measure  for  the  error  in  £'A.n. 
Another  convenient  test  follows  from  the  definition  of  e' 

.  0  k=0, 1, . . . ,n-2 

e'A^b  -  (36) 

1  k=n- 1 . 

This  relation  also  implies  that  e',  e^'A  . ..  e'An  ^  are  invariant  under  state 
feedback  (A,b)  -(A-bk',  b).  If  the  inverse  of  eq.  (32)  is  needed,  the  last 
row  of  E  can  be  brought  to  the  left  side 
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k'-e'An  =  £,W~\  p'  =  [po  p^  ...  pn_^] »  W"1 


e ' 
e’A 


*,.n-l 
e  'A 


£'  =  W(k'-e'An). 


(37) 


It  was  shown  in  [33]  that  the  columns  of  W  can  be  evaluated  recursively  by 
Leverrier's  algorithm,  which  also  gives  the  coefficients  a^  of  the  charac¬ 
teristic  polynomial  of  A,  det (XI-A)  =  a^  +  a^X  +  . . .  +  ^  +  \n.  In  this 

case  it  is  more  convenient  to  express  the  last  row  of  E  by  the  Cayley-Hamilton 
theorem  in  terms  of  the  previous  ones.  Then 


is'  =  (£'-£' )W-1,  a1  =  [ap  a1***an..1J 

js'  *  Wk*  +  a ' 


(38) 


4.2.2.  Partial  Pole  and  Gain  Assignment,  Output  Feedback 

So  far  the  mapping  from  &  to  X  has  been  formulated.  This  is  easily 
extended  to  mapping  from  A  -  ...  Xn)  to  X .  P(X)  may  be  written  as 

P(\)  =  (X-\1)(\-X2)  ...  (\-,\n).  (39) 

Then  by  eq.  (34) 

k'  =  e'  (A-^I)  (A-\2I)  . . .  (A-XnI)  .  (40) 

This  form  of  the  result  admits  consecutive  assignment  of  one  eigenvalue  after 
the  other.  P(X)  may  be  factorized  in  any  form 


P(X)  =  Q (X)R(X) 

R(X)  =  r  +  r.X+...+r  ,X°  ^  +  X® 

o  1  g-1 

Q(^)  -  %  +  <i]}+  ...  +Vg-l>n"8"1  +  xn"S* 


(41) 
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Let  for  example  the  n-g  roots  of  Q(X)  be  fixed,  then 


k*  =  e'P(A)  =  e 'Q (A)R(A)  *  e'QR(A) 


*  tr0  •••  rg-i  ll 


s’, 


e  ’A® 


— *Eq  =  [£’  13Eq- 


(42) 


It  is  also  possible  to  fix  g  gains  k^=£*n^,  where  n^  is  the  ith  column  of  E. 
These  g  linear  equations  in  the  p  coefficients  may  be  used  to  eliminate  jr 1 
[55].  £*  is  written  as 


£*=[Pr 


P  i  •  •  •  r  -i  2-3 

rn-l  1  o  g-1 


qi  •••  Vg-1  1 


i  o 


qn-g-l 


*q  ,1 

n-g- 1  ! 


£*  =  [r*  1] 


where  S  is  a  g  x(n+l)  matrix  and  _t’  a  1  x(n+l) vector.  Let  k£  be  the  fixed  gains, 
which  for  convenience,  are  chosen  to  be  the  last  g  gains  in  k’.  Then 

"S' 

Je’  -  V£K]  =  &  1]l 


t' 


K  =  —  — b  b 


which  can  be  solved  for 


-i 


(43) 


if  the  gxg  matrix  S  is  invertible.  Note  that  this  condition  does  not 


depend  on  the  vaLues  of  k'  Thus  this  is  the  same  problem  as  in  output  feedback, 
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k£  = 0,  where  certain  pole  locations  cannot  be  achieved  [51],  The  singularities 
of  SE^will  require  further  study,  is  determined  by 

&  -  £*E  -  (£'S  +  t*)E  .  (44) 


Assigning  the  n-g  eigenvalues  of  Q (X.)  determines  j>.  The  remaining  eigenvalues 
can  be  determined  by  factoring  the  residual  polynomial  R(\)  with  coefficients 
given  by  eq.  (43). 

Example  4.3:  (partial  pole  placement)  For  the  crane  let 

P(s)  =  (qQ +  qLs +s2)(rQ +  rLs +s2)  (45) 

where  q  and  q  are  fixed  and  r  and  r,  remain  as  free  parameters, 
o  l  o  i  _  _ 

e* 

«Q  =  e'Q(A)  =  [qQ  1]  e'A 


T[,o  h 


k«  =  e’  .  R(A) 


q0  qL  ^qo-g  iqL 

im 

“  —  Iro  rl  1]  0  qo  _ql8  iqo'g 

_°  0  (^-qo)g  -qLS 

Let  for  example  the  natural  frequency  m  of  the  pendulum  be  unchanged,  but 

2  2 

introduce  a  damping  d,  i.e.  Q(s)  =  u.p  +  2du)pS  +  s  ,  then 


rnC+T\ 

c 

m  Z 

1j 

c  Z 

is’  •  [ro  rj  l] 

0 

”c™L 

-eg 

ra  Z 

L 

i  ° 

0 

0 

-eg 
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with  c  *  2d  V  (mc+iaL)mc/gi. 


Example  4.4:  (fixed  feedback  gains)  For  the  crane  of  Example  4.2  let  and 


k4  be  fixed,  i.e.  kj[  =  [k2  k3]  ,  1^  =  1^  k^] 


4mc/g 


5.  = 


J  g 


(mc-HnL)g 


4  = 


fmc/g 


f2mc/g 


Then  by  eq.  (43) 


r‘  -  tr0  r].]  =  tki  k4+q1J&ncl* 


q0Wg  qi^  mc/g 


^nc(q0^/g-i) 


The  inverse  exists  if  q  /g/i  and  q  £0.  Then 

o  o 


kls 

ro  q  in,  ’ 


qlmc'kiql/qo+k4^ 

mc  (qo^/ 8_1) 


and  with  eq.  (44) 


2  *  q*-g  (49> 

°  2 

qig  k,q,  k 

k,  =  i(k,  /q  -m  )  (q  +  ; —  )  +  g  ( —  - h  m  +  itl  ) 

3  i  o  c/v^o  q0^-g  q0f-g  qQ  c  V 

k^  will  be  fixed  by  the  following  consideration:  Assume  a  force  limitation 
lu(t)  I  sU  for  all  t  for  a  typical  operation  of  the  crane,  i.e.  a  displacement 
of  a  load  at  rest,  x(0)  = [L  0  0  0] ' ,  L  > 0  (e.g.  length  of  a  loading  bridge) 
to  a  final  position  x(t  )  “  [0  0  0  0] * .  Typical  responses  of  sufficiently 
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stabilized  cranes  show  an  initial  peak  u(0)  of  the  force  as  the  maximum  value 
of  |u(t)l.  A  simple  approach  to  avoid  saturation  is  therefore  to  meet  a 
necessary  condition  by  fixing  |u(0)l  =U  and  checking  the  conditions  for 
u(0)/u(0)  SO.  Then  I u ( t )  i  for  t>0  may  be  checked  in  a  simulation.  Here 


u (0)  =  -k'x(0)  =  -kjL 

u(0)  =  -k'  (A-b  k1 )  x(0)  =  L  k^k,,  -  k^/jO/m^  =  I,  k^p^ 


(50) 


Thus  u(0)/u(0)  = -1/p^ £ 0  for  all  stabilizing  feedbacks  and  lu(0)l  =U  results 
in  k^  = U/L. 

It  is  desirable  to  avoid  the  difficult  measurement  of  the  rope 
angular  velocity  x^=9.  Thus  k^  =  0  is  chosen.  Then  by  eq.  (48) 


r  =  Ha _ 

o  Lim  q  * 
cno 


rl  = 


q1g(l-U/Lmcqo) 


xqQ-g 


(51) 


rQ  and  r^  are  the  coefficients  of  the  residual  polynomial,  which  is  obtained 

after  q  and  q,  have  also  been  fixed.  Necessary  and  sufficient  conditions  for 
o  i 

stability  are  q  >0,  q.  >0,  r  >0,  i’  >0.  With  eq.  (49) 


k2 


%qlmci_U8ql/Lqo 


V"§ 


(52) 


k3  = 


qlS 


^g 


-im  (q  +  — )  +  (q  i-g  + - : — )  +  (m  +m_  )g 

c  o  q  4-g'  Lq  o  &  q  £-g'  c  l/° 


4.2.3.  Sensitivities,  Incremental  Stabilization 

The  influence  of  a  coefficient  p^  of  the  characteristic  polynomial 
on  k',  given  the  other  p  ,  follows  from  eq.  (32)  as 
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ok' 

■r —  =  e'A  . 

dp. 


(53) 


The  influence  of  an  eigenvalue  X^  on  k',  given  the  other  X  ,  is  by  eq.  (40) 


ok' 

dT 


-  -e'CA-X^)  ...  (A-Xi_1I)  (A-Xi+1I)  ...  (A-XnI), 


(54) 


For  complex  conjugate  eigenvalues  quadratic  factors  in  P(X)  are  more  convenient 
Let  P(X)  = (a+bX+X2)R(X),  then 


k'  =  e' (al +  bA  +  A  )*R(A) 


Ok1 

da 


=  e'R(A) 


dk* 

^-  =  e'AR(A), 


(55) 

(56) 


Example  4.5;  (Incremental  stabilization,  global  robustness)  For  the  crane  of 

2  2  2 

eqs.  (2  6)  and  (35).  the  open  loop  characteristic  equation  is  s  (s  +uo  )  =0. 

Find  a  small  stabilizing  feedback  Ak1  with  the  least  number  of  required  sensors 


2  2  2 

P(s)  =  (s  +as+b)  (s  +cs+ujp+d)  with  small  a  >  0,b  >  0,c  >  0,  and  small  d 


k'  =  e'  (A2+aA+bI)  [A^+cA-Ku?+d)lj 

m  i  _  dk1  .  dk{  .  ok1  ,  ok' 

\k  —  —  •  a  4-  — —  *  b  +  t —  •  C  +  -r-: —  •  d 

3a  3b  dc  3d 


(57) 

(58) 


~l  =  |-  e'(A2+aA)(A2+oi2l)  =  e'(A3+w2A) 

03  b=c=d=0  °a  “  ~  P- 


dk' 

db 


a=c=d=0 


=  |e'(A  +bl)  (A^-Hjo^)  =  e'  (A  +x\L) 


P~ 


ok' 

dc 


a=b=d=0 


z  9  ?  ^  o 

—  e 1  A“  (A+cA-tVl)  =  e'A 
oc - -  -  p-'  - 


dk' 

ad 


a=b=c=0 


d  2  2  9  2 

=  e'A  [A  +(<u  +d)I]  =  e'A 


Ak' 


Ak 


1 


Ak4 


=  e ' 


.[(a+c)A  +  (b+d)A  +  auipA  +  btu^l] 


=  b(jUpimc/g  *  (mc-hnL).b 
2 

=  aoj^im^/g  =  (m^-hn^)  -  a 

=  bcu2i2m  /g  -  (b+d)im  = 
pc  c 

2  2 

=  aiu  &  m  /g  -  (a+c)im  = 
p  c  c 


i(mLb-mcd) 
i(mLa-mcc) . 


With  a>0,  b>0,  c>0  this  is  the  cone  of  stabilizing  directions  at  the  origin 
of  the  four  dimensional  \  space.  It  includes  the  directions  Ak^ = 0  and  Ak^  =  0, 
i.e.  no  feedback  of  the  rope  angle  is  necessary  with 


d 


(59) 


No  feedback  of  the  rope  angular  velocity  is  necessary  with 


ra 


c 

Ak^  and  Ak^  must  be  positive  for  stabilization,  i.e.  crab 
velocity  must  be  available  for  feedback.  Output  feedback 


(60) 

position  and  crab 


Ak*  *  (mc+m^)[b  a  0  0] 


(61) 


with  small  positive  a  and  b  then  stabilizes  the  system.  For  sufficiently  small 
a  and  b  the  characteristic  polynomial  is  arbitrarily  close  to 


2  ?  "v,  ?  ny 

P(s)  =  (s  +as+b)  (s“  H - as  +fx“  4 - b)  , 

m  pm 

c  c 


(62) 


If  the  physical  parameters  of  the  crane  are  unknown,  and  output  feedback 


Ak*  =  [Ak^  Ak0  0  0]  with  small  Ak^  >0, 


(63) 
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is  applied,  then  in  eq.  (61)  a  -  k2 / (ln/fn^)  and  b  =  kjV  (n^dn^)  are  unknown,  but 
positive,  i.e.  with  the  feedback  of  eq.  (63)  stability  is  robust  with  respect 
to  arbitrary  changes  of  load  mass  m^*  crab  mass  mc,  rope  length  l ,  and  gravi¬ 
tational  constant  g.  Thus  in  this  example  global  robustness  of  stability 
with  respect  to  perturbations  in  four  directions  is  achieved. 

Note  that  it  is  possible  to  destabilize  the  model  A  and  in  eq.  (26) 
with  feedback  (63),  if  small  changes  in  arbitrary  coefficients  of  A  and  b  are 
permitted.  The  general  assumption  in  this  paper  is  that  A  and  b  do  not  change 
arbitrarily  but  in  known  directions.  Only  by  this  assumption  does  it  become 
possible  to  accomodate  large  parameter  variations. 

4.3.  Mapping  of  Boundaries 

4.3.1.  Mapping  from  X  plane  to  X  space 
In  the  complex  plane  of  eigenvalues 

\  =  v  +  jw  (64) 

2  2 

boundaries  w  =w  (v)  (i.e.  symmetric  with  respect  to  the  real  axis)  are  of 
interest,  which  are  related  to  desirable  system  properties,  e.g.  stability, 
damping,  bounds  on  the  natural  frequency.  Examples  were  given  in  Figs.  4.1 
and  4.2. 

Assume,  due  to  a  change  of  the  state  feedback  gains  k',  a  real  root 
crosses  the  boundary  at  its  intersection  with  the  real  axis  at  ,\  =  vD.  Then  in 
i?  space  a  boundary 

P(A)  =  (\-vr).R(,\),  R(.\)  =ro  +  r]X+...  +rn_2,\n"2 +\n_1  (65) 

is  crossed,  which  is  linear  in  its  parameters  r  ,  r,,...,r  ,  i.e.  it  is  an 

r  o  1  *  n-2 
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(n-l)  dimensional  hyperplane.  Due  to  the  linearity  of  k'  =£*E.®[£'  I]E, 
this  boundary  maps  into  a  hyperplane  in  X  space,  which  is  crossed  by  the 
feedback  gain  vector  k'. 

Assume,  due  to  a  change  of  the  state  feedback  gains  k'  a  complex 
conjugate  pair  crosses  the  X-boundary  at  v+jw.  Then  in  9  space  a  boundary 

P(X)  = (X-v-jw) (X-v+jw) *R(X)  =Q(X)R(X)  (66) 

where 

Q(X)  =  X2  -  2vX  +  v2  +  w2  (v) 

and 

R(X)  =rQ  +  r1X+  ...  +rn_3Xn  3+Xn"2 

is  crossed.  For  a  fixed  v,  i.e.  a  fixed  pair  of  eigenvalues  on  the  boundary, 

this  is  a  (n-2)  dimensional  hyperplane  in  and  X  space.  However,  for 

different  values  of  v  different  hyperplanes  are  obtained.  For  fixed  R(X)  the 

boundary  line  in  and  X  space  is  obtained  by  moving  a  pair  of  conjugate  eigen- 

2 

values  along  the  boundary.  In  this  case  the  form  of  the  boundary  w  (v)  in  eq. 
(66)  determines  the  shape  of  the  boundary  line  in  Q  and  X  space. 

Some  boundaries  of  particular  interest  are: 

1)  Imaginary  axis,  stability  boundary  in  s-plane,  v = 0,  Q( 

9  2  2 

Q(X)  = X~+w  ,  boundary  linear  in  w  . 

2)  Parallel  to  imaginary  axis 

.2  2  2  2 

v  =  v^,  Q(a .)=/.  -2v^\  +  v^+w  ,  linear  in  w  . 

3)  Conic  section  symmetric  to  the  reaL  axis,  i.e. 

2  2 

w=c+cv  +  c'/.  (6/ ) 

o  1  2  v  ' 


SpeciaL  cases  are 
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c2<0  ellipse,  of  particular  interest  are  circles  c2  =  -1,  e.g. 
constant  natural  frequency  curves  in  s-plane,  stability 
limit  and  other  boundaries  in  z-plane,  see  Fig.  1. 
c2  =  0  parabola,  or  if  also  0^  =  0,  cq  >0  straight  line  parallel  to 
the  real  axis.  For  co=c^  =  c2=0  boundary  between  real  and 
complex  eigenvalues. 

2  2 

c2  >0  hyperbola,  in  particular  2  straight  lines  for  w  =  c2(v-vq)  > 
c2  >0,  e.g.  constant  damping  lines  in  s-plane.  This 
boundary  is  frequently  combined  with  a  parallel  to  the 
imaginary  axis.  Here  it  is  more  convenient  to  use  a  hyperbola, 
which  guarantees  the  required  damping  and  minimum  negative 
real  part  of  the  eigenvalues,  see  Fig.  2. 

Substituting  eq.  (67)  into  Q(\)  from  eq.  (66)  gives 

Q(A.)  =\2  -  2v\+  (l+c2)v2  +CjV  +  co.  (68) 

The  boundary  is  quadratic  in  v.  It  becomes  linear  only  if  c^-l,  i.e.  for  a 
circular  boundary  in  \  plane.  In  other  words:  If  n-2  roots  in  R(\)  are  fixed 
and  the  remaining  two  roots  of  P(\)  move  as  a  conjugate  pair  along  any  circle 
in  the  X-plane  with  center  on  the  real  axis,  then  the  corresponding  point  in 
&  and  X  space  moves  along  a  straight  line.  This  is  the  reason  for  the  proposal 
of  a  family  of  circls  Tr  in  the  z-plane  in  Fig.  4.1.  Its  equation  is 

2  2 

(v-v  )  +  w  =  r 
o 

v  (v  -1)  =0.99r(r-l),  v  <0.5  for  r<l 
0  0  *  o 

v-0  for  r  2  1. 

o 


(69) 
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For  r  =  0  it  is  the  deadbeat  solution  with  all  eigenvalues  at  z  =  0.  With 
increasing  r  the  center  vq  of  the  circles  moves  to  the  right  until  it  reaches 
0.45  for  r»0.5,  it  then  goes  back  to  zero  to  produce  the  unit  circle  for 
r = 1.  If  boundaries  in  the  unstable  region  are  needed,  concentric  circles 
with  radius  r  may  be  used.  We  may  begin  with  a  radius  r  such  that  all  open 
loop  eigenvalues  are  enclosed  by  the  circle,  and  design  the  feedback  such  that 
r  is  reduced  to  a  radius  smaller  than  1.  In  the  further  reduction  of  r  at 
tradeoff  with  the  required  gains  must  be  made.  For  continuous -time  systems 
the  family  of  hyperbolas  r  of  Fig.  4.2  in  the  s  (=a+juj)-plane  may  be  used. 

Its  equation  is 


2  2  2  2 

oj  = -p  +c  / p  for  a< 0 
a  *  -p  for  cr  >0. 


(70) 


For  large  P  an  extremely  fast  solution  is  obtained,  P  = 1  gives  the  l/y2  damping 

line  as  asymptotes,  for  p  -  0  it  goes  to  the  imaginary  axis.  Negative  p 

represent  parallels  to  the  imaginary  axis  in  the  right  half  plane.  Beginning 

with  a  sufficiently  negative  p  to  have  all  open  loop  eigenvalues  to  the  left 

of  the  boundary  T  ,  p  may  be  increased  by  the  feedback  to  positive  values, 

P 

where  again  a  tradeoff  with  the  required  gains  must  be  made. 

Besides  the  real  and  the  complex  root  boundary  there  is  the  third 
possibility  of  a  root  leaving  the  region  through  infinity.  This  can  be 
avoided  by  closing  the  contour  in  the  s-plane  by  an  arc  of  a  circle  with  large 
radius.  Practically  this  may  be  a  circle  corresponding  to  the  design  bandwidth. 
This  is  of  particular  importance  if  we  need  a  40  db/decade  "roll  off"  band¬ 
width  limitation,  e.g.  to  avoid  interaction  with  unmodeled  modes  at  higher 
frequencies . 
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Typical  desirable  regions  for  the  eigenvalue  location  in  the  s-  or 
z-plane  are  connected  and  have  two  intersections  with  the  real  axis.  In  this 
case  there  are  two  real  root  boundaries  and  a  (possibly  piecewise  defined) 
complex  root  boundary  in  9  and  X  space. 

Equations  (65)  and  (66)  show  that  the  mapped  boundaries  in  9  and  X 
space  represent  the  conditions  under  which  the  number  of  eigenvalues  inside 
and  outside  a  \-region  can  change.  The  boundaries  partition  the  9  and  X 
space  into  regions,  each  of  them  corresponds  to  a  fixed  number  of  eigenvalues 
inside  the  \  region,  and  it  must  be  decided  for  which  9  or  ?C  region  all  eigen¬ 
values  are  inside  the  \  region.  For  closed  contours  in  the  X,-plane  the  X 
region  is  bounded,  since  by  eq.  (32)  no  k^  can  go  to  infinity.  If  there  are 
several  bounded  regions,  a  simple  test  is  to  check  the  eigenvalues  for  an 
arbitrary  K'  in  the  considered  X  region.  An  alternative  are  Siljak's  "shading 
rules"  for  the  boundaries  [49] . 

For  second  and  third  order  systems  it  is  possible  to  visualize 
regions  in  X  space  graphically.  This  is  done  in  the  following  for  the  unit 
circle,  i.e.  the  stability  region  of  discrete  systems  in  X  space  is  determined. 
Figure  4.3  shows  the  regions  in  the  k^-k^-plane  for  a  second  order  system.  The 
two  real  root  boundaries  are  two  infinite  straight  lines  intersecting  at  B. 

Thus  B  can  be  obtained  by  placing  one  pole  at  z  =  -l  and  one  at  z  =  1.  The 
complex  root  boundary  is  the  straight  line  AC.  A  is  obtained  by  placing  a 
double  pole  at  z  =  -l  and  C  by  a  double  pole  at  z  =  1.  Thus  Fig.  4.3  is  completely 
determined  by  three  pole  placements 


'  • 
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[1 

2 

1]E 

— B  = 

[-1 

0 

1]E 

(71) 

[1 

-2 

1]E. 

The  boundaries  partition  the  k^-^-plane  into  five  regions  with  the  properties 
1)  both  poles  inside  the  unit  circle,  2)  one  left,  one  inside,  3)  one  left, 
one  right,  4)  one  inside,  one  right,  5)  complex  outside  or  both  left  or  both 
right  (a  distinction  between  these  three  cases  in  region  five  would  require  a 
further  boundary  distinguishing  real  and  complex  roots).  Usually  only  the 
stability  region  1  is  of  interest. 

Only  the  stability  region  will  be  determined  now  for  third  order 
systems.  It  is  shown  in  Fig.  4.4.  The  two  real  root  boundaries  are  the  two 
planes  in  which  the  triangles  ABC  and  BCD  are  contained.  They  intersect  along 
the  straight  line  BC.  B  is  obtained  by  placing  two  poles  at  z  =  -l  and  one  at 
z  *  1,  C  corresponds  to  one  pole  at  z  =  -1  and  two  at  z  =  1.  For  any  fixed  real 
eigenvalue  and  the  two  others  moving  along  the  unit  circle  a  straight  line  is 
obtained.  Thus  the  complex  root  boundary  may  be  visualized  as  being  generated 
by  a  moving  straight  line  from  a  point  on  AB  to  a  point  on  CD.  It  moves  as 
the  real  eigenvalue  moves  from  -1  to  +1.  A  corresponds  to  a  triple  eigenvalue 
at  z  =  -l  and  D  to  a  triple  eigenvalue  at  z=+l.  Thus  the  vertices  of  the 
stability  region  are  obtained  as 


k^[l 

3 

3 

1]E 

k'  =  [-l 

-1 

1 

1]E 

i— ( 

ua 

il 

zX 

-1 

-1 

1]E 

i^D  =  [-l 

3 

-3 

1]E. 

(72) 


Apparently  the  tetrahedron  ABCD  is  a  convex  hull  for  the  stability  region. 

For  the  £  space  Fam  and  Meditch  [38]  showed  that  this  property  generalizes 
to  arbitrary  degree  n  of  the  characteristic  polynomial. 

Theo  rem  2  (Fam,  Meditch):  For  an  n-th  order  discrete  system  a  convex  hull 
of  the  stability  region  in  &  space  is  a  polyhedron  whose  vertices  correspond 
to  the  n+1  polynomials  with  zeros  in  the  set  {-1,  l}. 

Proof:  See  [38]. 

Corollary:  If  the  unit  circle  is  replaced  by  a  circle  with  center  vq  and  radius 
r  in  the  z-plane,  which  intersects  the  real  axis  at  v,  =v  -r  and  v„=  v  +r, 

then  the  vertices  of  the  convex  hull  of  the  corresponding  region  in  @  space 

correspond  to  the  n+1  polynomials  with  zeros  in  the  set  [v^,v^j.  This  may  be 
shown  by  reducing  this  problem  to  the  previous  one  via  z'  = (z-vo>/r. 

Theorem  3:  A  necessary  condition  for  all  roots  of  det (z]>A+b  k' )  =  0,  (A,_b 
controllable)  to  be  inside  a  circle  with  real  center  and  real  axis  inter¬ 
sections  at  z=v^  and  z  » v^  is  that  k1  is  in  a  polyhedron  in  X  space,  whose 
vertices  are  obtained  by  pole  placement  of  the  n+1  polynomials  with  zeros  in 
the  set  £ , v2 } * 

Proof :  Follows  from  Theorem  2  and  the  linearity  of  the  map  k1  =p*E. 

If  the  circle  in  z-plane  is  deformed  to  a  different  closed  contour 

with  the  same  real  axis  intersections  at  v^  and  ,  then  the  certices  of  the 

region  in  \  space  and  the  two  real  root  boundaries  remain  unchanged,  only  the 
hypersurface  for  the  complex  root  boundary  is  deformed.  It  is  an  open  question, 
how  far  the  region  in  X  space  can  be  extended  in  the  k-plane  such  that  k' 
remains  in  the  previously  described  polyhedron.  The  inverse  problem  (given  a 
region  in  X,  which  is  the  corresponding  region  in  the  \-plane?)  leads  to 
sufficient  conditions  on  k'  to  place  all  eigenvalues  in  the  resulting  region. 
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This  problem  was  studied  by  Marden  [52,  Theorem  8.2],  However,  these  regions 
in  X-plane  are  not  nicely  related  to  desirable  dynamic  properties.  For  the 
polyhedron  discussed  above  the  region  in  X-space  is  the  union  of  two  circles 
with  complex  conjugate  centers  and  real  axis  intersections  at  v^  and  v^. 

Marden  showed  that  this  region  cannot  be  reduced.  Note  that  in  this  direction 
from  the  n-dimensional  ^-boundary  to  the  two-dimensional  X-boundary  no  necessary 
and  sufficient  conditions  can  be  obtained  for  an  arbitrary  defined  X  region. 

This  problem  is  overdetermined.  Thus  it  is  advisable  to  assume  a  region  in 
X-plane  and  to  determine  the  necessary  and  sufficient  conditions  on  k'. 

4.3.2.  Mapping  to  a  Subspace  of  X 

Some  gains  may  be  fixed  to  zero,  like  in  output  feedback  or  under 
sensor  failures,  or  to  some  other  values  like  in  Example  4.4.  This  means  that 
we  are  looking  for  a  solution  in  a  subspace  of  Such  a  solution  may  not 
exist;  take  for  example  Fig.  4.3  and  fix  k^  to  be  bigger  than  k^(A)*  Then 
there  does  not  exist  a  stabilizing  k0 .  The  set  of  admissible  solutions  may 
also  become  disconnected,  even  if  it  was  connected  in  X  space;  take  for  example 
Fig.  4.4  and  fix  ^=0  such  that  the  plane  k£  =  c  intersects  the  two  tips  of  the 
stability  region. 

Example  4.6:  (Disconnected  stability  regions  in  a  subspace  of  X) 


0 

1 

0 

o" 

x(k+l)  = 

0 

0 

1 

x(k)  + 

0 

_0.  6 

_2 

2.1_ 

1 

u(k). 
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The  system  is  open  loop  unstable  (eigenvalues  z^  =  0.5,  z0  ^=0.8+j-v  0.5b). 
Fix  k9  = 0  (output  feedback)  and  find  the  set  of  stabilizing  gains  in  the 
k^-k^-plane.  The  real  root  boundaries  are  the  straight  lines 


for  z  =  1 


for  z  =  -1 


k3+  =  "kl  ‘O'3 


=  -k^  +  5.7 


and  the  complex  root  boundary  is  the  hyperbola 


k3c  =k1  +  1.5  +  l/(k1-0.6). 


Figure  4.5  shows  the  three  boundaries  and  the  two  disconnected  stabilizing 
regions.  Its  vertices  are 


0.1  -0.4 


Nonconvex  and  disconnected  solution  sets  like  in  this  example  lead  to  difficulties 
in  numerical  algorithms.  Sirisena  and  Choi  [53]  formulate  the  problem  of 
placing  poles  in  a  specified  region  by  output  feedback  as  minimization  of  a 
function  J,  which  becomes  zero,  if  a  solution  is  found.  Their  conclusion  from 
computational  experience  is:  "If  however  a  local  (nonzero)  minimum  of  J  is 
reached,  the  algorithm  should  be  restarted  with  a  different  initial  value  of 
the  feedback  matrix.  Repeated  failure  to  reduce  J  to  zero  would  indicate  the 
absence  of  a  solution".  If  we  want  to  find  the  set  of  admissible  solutions,  a 
systematric  or  random  search  in  the  appropriate  subspace  may  be  necessary.  In 
order  to  restrict  the  search  to  a  promising  region,  necessary  conditions,  like 
the  one  provided  in  Theorem  3,  are  very  useful.  In  Example  4.6  this  is  the 
quadrangle  EFHG,  more  general  the  polyhedral  cross-section  of  the  subspace  with 


the  polyhedron  in  the  n  dimensional  space.  If  no  such  cross-section  exists, 
it  can  be  immediately  concluded  that  no  solution  exists.  The  example  indicates 
that  points  near  the  real  root  boundaries  are  promising  candidates. 


For  a  fourth  order  system  there  exist  gains,  for  which  two  complex 
pairs  of  eigenvalues  cross  the  boundary  simultaneously.  Here  the  complex 
root  boundary  in  X  space  intersects  itself.  If  the  two  crossing  points  on 
the  boundary  approach  each  other  and  eventually  become  a  crossing  of  a  double 
pole,  then  the  bow  in  the  boundary  becomes  a  cusp. 

Example  4.7:  (Partial  gain  fixing)  For  the  crane  of  Example  4  let  m  = 1000  kg, 

2 

l  = 10  m,  g = 10  m/sec  ,  U  = 5000  Newton,  L = 10  m.  The  load  mass  m^  is  unknown. 
Example  2  showed  that  only  k^  = 


can  be  determined.  For  k^=U/L  =  500  and 
k^  = 0  find  the  region  in  the  k^-k^Q-plane,  for  which  all  eigenvalues  are  left 
of  the  hyperbola 

u  =  (2o)2  -  1/22  (74) 


knowledge  of  m^  only  =  k^-lOm^ 


k^Q+^g  depends  on  the  load  mass  m^.  Without 


in  the  s-plane.  Then  for  the  complex  root  boundary  from  eq.  (66) 

qQ  = 5ct2  -  0.25,  q L  =  -2c . 

and  by  eqs.  (51)  and  (52) 


r 

o 


1/2 V  ri 


qL(l-l/2qo) 

v1”- 


lOOOq 

-^r(vl/2V 

o 


k30c  =  L0000(l/2qo-l)[qo-l+q2/(Vl)] 


(75) 


(76) 


(.77) 


(78 
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The  nice  stability  region  will  be  constructed  in  the  l^-k^Q-plane.  The  complex 
root  boundary  k2(a),  k^Ca)  is  obtained  by  substituting  values  ct£-0.25  into 
eq.  (75)  and  qQ  and  q^  into  eqs.  (77)  and  (78).  The  real  root  boundary  at 
ctr  =  -0.25  follows  from  eq.  (65)  with  ^  =  500,  k^  =  0  as  the  straight  line 

k3R  "  k30R  +  10  V  k30R  =  95625  -42.5  V  (79) 

Both  boundaries  are  shown  in  Fig.  4.6.  For  cr  =  -0.25  the  complex  root  boundary 
starts  at  point  A.  With  increasing  a  it  goes  through  point  B  and  for  cr-*-0.5} 
(i.e.  q^  1)  to  infinity.  In  general  this  singularity  occurs  at  q^  =  g!  1.  For 
a  <-0.5  the  complex  root  boundary  returns  from  the  opposite  side  to  intersect 
the  real  root  boundary  at  C  and  itself  at  B. 

Note  that  the  characteristic  polynomial  is  obtained  by  eqs.  (75)  and 
(76)  in  factorized  form.  Thus  the  determination  of  the  eigenvalues  is  easy. 
They  are  given  together  with  the  k2  and  k^^  coordinates  in  the  following  table. 


k2 

k30 

Eigenvalues 

A 

4233 

-84292 

s1  2  = -0.25,  s3  4  = -1.867+j2.125 

B 

2367 

-35012 

s^  ^  = -0.275+J0.231,  s^  4 = -0. 908+j 1. 746 

C 

2769 

-22056 

s :  =  -0.25,  s2  =-1.337,  s3  4  =  -0.591+j  1.071 

At  A  the  real  and  complex  root  boundary  intersect,  i.e.  there  is  a  double  pole 
at  s.  1  =-0.25.  At  B  the  complex  boundary  intersects  itself,  i.e.  here  we 
have  two  complex  pairs  of  eigenvalues  crossing  the  boundary  simultaneously. 

At  C  a  real  root  at  s  = -0.25  crosses  simultaneously  with  the  complex  pair  4 
The  shaded  region  with  vertices  A,  B  and  C  corresponds  to  eigenvalues  to  the 
left  of  the  hyperbola  in  s-plane. 
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4.4.  Robustness  with  Respect  to  Large  Parameter  Variations 

Most  existing  methods  for  the  design  of  parameter  insensitive  control 
systems  try  to  achieve  robustness  of  a  system  property,  like  stability,  with 
respect  to  small  perturbations  of  plant  parameters  in  unknown  or  conveniently 
analyzable  directions,  e.g.  gain  or  phase  variations.  Such  methods  lead  to 
conservative  results  when  applied  to  problems  with  large  perturbations  in  known 
directions,  e.g.  for  the  crane  with  widely  varying  load  mass  or  rope  length  or 
for  an  aircraft  with  widely  varying  altitudes  and  speeds.  Let  the  plant  model 
in  sensor  coordinates 


x  =  Alj^x  +  B(9)u 


(80) 


be  given  for  several  typical  values  of  the  physical  parameter  vector  d_,  i.e. 
Aj=A(9j),  J3j=b(9j),  j=l,2,...,J.  A  fixed  state  feedback  k'  is  sought,  such 
that  all  eigenvalues  of  (A^-b^k')  are  located  in  a  specified  region  F  in  ^-plane. 

For  each  pair  ,b ^  we  obtain  a  different  matrix  E  ^ ,  which  maps  a 

desired  region  from  the  canonical  parameter  space  into  the  corresponding 
region  in  ?C-space  v*-a  k1  =jE*JF.  The  set  of  solutions  to  the  above  problem, 
if  it  exists,  is  the  intersection  of  ail  regions  in  X.  space.  If  no  inter¬ 
section  for  all  j=l,2,...,J  exists,  then  it  can  be  tested  whether  at  least  a 
group  of  plant  models  can  be  nicely  stabilized  with  one  gain,  and  it  may  be 
necessary  to  switch  to  a  different  gain  for  a  different  group  of  plant  models. 
Example  4.8:  (Maximum  parameter  variation,  gain  scheduling)  Let  the  mass  m^ 
of  the  crane  load  be  an  unknown  constant  between  the  weight  of  the  empty  hook 
and  the  maximum  load  for  which  the  crane  is  designed.  Under  the  constraints 
of  Example  4.7  find  the  fixed  gain  controller  which  accomodates  the  largest 


load  variation. 
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The  load  mass  enters  only  into  k^  = 
origin  of  the  k-j-axis  is  identical  with  k3Q  =  0  for  ^  =  0.  With  increasing 
load  mass  the  shape  of  the  region  of  nice  stability  is  unchanged,  but  it  is 
moved  upwards  by  10  in  the  k^k^-plane,  or  equivalently  the  origin  of  the 
k^-axis  is  moved  downwards  by  10  m^  in  the  k^-k^Q-plane.  Thus  for  load  varia¬ 
tions  of  cranes  it  is  not  necessary  to  plot  the  shifted  diagrams  in  order  to 
find  the  intersection.  The  largest  load  variation  can  be  accommodated  at  the 
largest  extension  of  the  nice  stability  region  in  k^  direction,  this  is  between 
C  and  D.  D  has  the  coordinates  ^=  2769,  k^Q  =  -45503  and  corresponds  to  the 
eigenvalues  s^  2  = -0.267+j0. 680  and  s^  ^  =  1. 118+j 1.872 .  Thus  k2  is  chosen  as 
2769.  This  results  in  an  admissible  load  variation  *  (-2205(H-45503)/10 
=  2344.7  %  2345kg.  Assuming  the  weight  of  the  empty  hook  is  50  kg,  then  ^=-21556 
puts  the  eigenvalues  for  11^  =  50  kg  at  s^  =  -0.25,  s2  *-1.337,  s^  ^  =  -0. 591+j  1. 071, 
where  s^  and  s^  ^  are  on  the  boundary  T.  For  1^  =  2395  kg  the  eigenvalues  are 
at  s ^  ^  =  -0.267+jO.  680  and  s^  ^  =  1. 118+j 1.872 ,  where  s^  2  is  on  the  boundary. 

In  summary:  The  solution 

k'  =  [500  2769  -21556  0]  (81) 

gives  the  following  properties  of  the  control  system. 

a)  Initial  peak  in  the  force  u  limited  to  500L,  where  L  is 
the  required  load  displacement. 

b)  No  measurement  or  estimation  of  the  rope  angular  velocity 
x^  = y  required. 

c)  Under  the  constraints  a)  and  b)  maximum  possible  load 

•?  2  ° 

variation.  The  eigenvalues  are  left  of  of"  =  (2c)“  -  1/2“  if 
and  only  if  50  kg<mr  <2  395  kg. 


k^Q  +  lO  m^.  In  Fig.  4.6  the 
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Now  assume  Chat  the  crane  is  designed  for  a  maximum  load  of  3500kg,  i.e.  a 
gain  scheduling  is  necessary.  The  second  load  range  may  be  chosen  as 
1155  kg  Crn^  <3500  kg,  i.e.  ^  =  -10506.  Then  for  50  kg<m^<1155  kg,  ^=-21556 
must  be  used  and  for  2395  kg  <m^  <3500  kg,  k^ = -10506.  For  the  overlapping 
range  1155  kg  <m^  <2395  kg  either  gain  is  good,  such  that  the  crane  operator 
can  switch  between  high  and  low  load  based  nn  his  very  crude  load  estimate, 
which  may  be  135%  wrong.  This  wide  overlap  provides  robustness  of  the  gain 
scheduling  scheme. 

If  the  rope  length  of  the  crane  is  varied,  the  shape  of  the  nice 
stability  region  in  Fig.  4.6  changes  and  an  intersection  of  various  regions 
must  ue  found. 


For  different  values  of  a  physical  parameter  vector  different 
regions  in  the  %-plane  may  be  given  and  the  intersections  of  the  corre¬ 
sponding  Te'-space  regions  may  be  found.  This  is  particularly  useful,  if  the 
plant  is  slow  for  some  parameter  values  and  fast  for  others  like  in  aircraft 
control.  A  general  recommendation  for  the  design  of  robust  control  systems 
with  input  constraints  is;  do  not  try  to  make  a  slow  plant  fast  or  a  fast 
plant  slow  by  feedback. 

The  graphical  determination  of  intersections  is  limited  to  two 
parameters  at  a  time.  Intersections  of  three  dimensional  regions  may  be  made 
visible  by  computer  graphic  methods.  In  situations  with  more  free  controller 
parameter  the  design  may  proceed  iteratively,  where  in  each  design  step  n-2 
feedback  gains  are  fixed  and  admissible  regions  in  the  plane  of  the  remaining 
two  feedback  gains  are  determined.  The  results  of  Section  4.3.1  on  mapping 
circular  boundaries  give  some  additional  insight,  which  is  useful  for  a  fully 
computerized  search  for  intersections  of  regions.  First  it  is  important  to 
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note  that  the  nice  stability  regions  are  not  convex  and  thus  their  intersection 
may  be  disconnected.  Consider  for  example  the  stability  region  for  a  third 
order  discrete-time  system  in  Fig.  4.4.  Assume  that  for  a  different  value  of 
the  parameter  vector  the  stability  region  is  turned  around  by  180  degrees 
such  that  the  two  tips  of  the  stability  regions  intersect.  Then  the  set  of 
solutions  is  disconnected,  even  for  full  state  feedback.  In  this  situation 
a  search  in  X  space  may  be  made.  Bounds  for  the  search  region  are  given 
by  the  following. 

J 

Theorem  4:  A  necessary  condition  for  all  roots  of  .tt  det  (zl-A .  +  b  .k1 )  =  0  to 
-  J  J=1  - J  “Ji¬ 
be  on  or  inside  a  circle  with  real  axis  intersections  at  z  =  v^  and  z  =  v^  is, 

that  k'  is  on  or  inside  the  intersection  of  J  polyhedra.  The  vertices  of  the 

j-th  polyhedron  are  obtained  by  assigning  all  (n+1)  polynomials 

P(z)  =  det(zI.-A.+b.k.')  with  zeros  in  the  set  [v.,v„j. 

J  J  ^ 

Proof:  Follows  from  Theorem  3. 

The  intersection  of  polyhedra  is  a  polyhedron  itself,  its  vertices 
are  promising  candidates  in  the  search  for  points  which  also  meet  the  sufficient 
conditions.  In  order  to  define  a  rectangular  grid  for  the  search  it  is 
convenient  to  put  each  polyhedron  into  the  smallest  box  with  surfaces  parallel 
to  the  axes  and  to  restrict  the  search  to  the  intersection  of  the  boxes. 


4.5.  Robustness  with  Respect  to  Sensor  Failures 

Sensor  failures  are  assumed  to  occur  in  the  form  that  the  sensor 
output  is  no  longer  correlated  with  the  measured  variabLe.  As  far  as  the 
characteristic  equation  is  concerned,  this  is  equivalent  to  having  a  sensor 
output  zero.  There  may  be  a  bias  or  other  noise  tern  introduced  by  the 
failed  sensor.  This  noise  term  can  be  considered  as  an  external  input.  This 
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may  require  that  the  failure  is  detected  and  the  failed  sensor  is  removed 
from  the  control  system.  Then  also  the  control  law  may  be  changed.  However, 
for  this  latter  decision  there  should  be  sufficient  time  to  come  to  a  reliable 
decision  without  false  alarms.  This  requires  that  after  the  failure  the  system 
at  least  remains  stable  with  some  stability  margin.  In  other  applications  it 
may  suffice  to  be  able  to  continue  the  mission  after  a  sensor  failure  without 
removal  of  the  failed  sensor,  e.g.  to  drive  an  automobile  safely  to  a  service 
station  to  get  a  broken  sensor  replaced,  such  that  optimal  fuel  economy, 
emission  control,  acceleration,  etc.  is  regained. 

The  robustness  problem  is:  Consider  M  failures  of  a  sensor  or 
combinations  of  sensors  leading  to  the  crippled  feedback  vectors  k^, 
m  =  1,2 , . . . ,M,  in  which  the  appropriate  elements  of  k'  are  replaced  by  zero. 

Find  k'  such  that  all  zeros  of 

M  J 

",  "  det(Xl-A. +b.k')  -  0  (82) 

m=l  j=l  - j  -j-m 

lie  in  an  "emergency  region"  T  in  \-plane.  The  emergency  specification  is 
robust  with  respect  to  a  failure  or  sensor  i  if  and  only  if  in  ?(-space  the 
projection  of  k*  into  the  subspace  k^ = 0  is  in  the  intersection  of  all  J 
emergency  regions. 

Figure  4.7  shows  an  example  of  the  intersection  of  emergency  and 
nominal  regions  in  the  k^-^-plane.  If  we  choose  k'  at  point  1,  then  the 
projection  on  the  k ^  axis  is  inside  the  emergency  boundary,  i.e.  f  is  robust 
with  respect  to  a  failure  of  sensor  1.  It  is,  however,  not  robust  with  respect 
to  a  failure  of  sensor  2,  since  the  projection  on  the  k^  axis  is  outside  the 
emergency  region.  Points  in  the  shaded  region  are  robust  with  respect  to 
failure  of  either  sensor.  For  no  k'  r  is  robust  with  respect  to  failures  of 
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sensors  1  and  2,  since  the  origin  k^ =  k2  = 0  lies  outside  the  emergency  region. 
Point  3  also  meets  the  nominal  specification  and  is  a  good  candidate  for  a 
robust  control  system.  Since  the  nominal  boundary  intersects  the  axis, 
an  alternative  to  the  robust  solution  3  is  to  eliminate  the  x^  sensor  and  to 
multiplex  the  x2  sensor.  This  would  maintain  the  nominal  specifications  under 
a  failure  of  one  of  the  sensors.  However,  it  requires  failure  detection 
with  at  least  three  sensors.  For  robustness  with  respect  to  sensor  failures 
a  dynamic  feedback  structure  like  in  eqs.  (24)  and  (25)  is  more  advantageous. 

In  Chapter  5  and  [54],  Franklin  designs  a  flight  control  system  with  dynamic 
feedback  such  that  emergency  conditions  are  robust  with  respect  to  an  accelero¬ 
meter  or  a  gyro  failure  in  different  flight  conditions. 

4.6.  Other  Features  of  ^-Space  Design 
4.6.1.  Input  Constraints 

Constraints  of  the  type  I  u.  ( t )  I  ^  U  for  all  t  or  |u(t)l-U  for  all 
t  can  be  indirectly  treated  in  X  space.  For  the  regulator  problem 

I u ( t )  |  =  Ik'x(t) I  £  ; ! x ( t ) !  1  (83) 

with  equality  for  the  worst  case  of  x(t)  (e.g.,  x = ck  for  some  c  ^ 0) .  Assuming 
that  all  state  variables  have  been  normalized  to  their  maximum  value,  the 
norm  ,]k:|  =  yk'k,  i.e.,  the  distance  from  the  origin  in  X  space  can  be  used  as 
a  measure  for  lul.  This  provides  a  criterion  for  the  selection  of  a  gain  from 
the  admissible  set:  Choose  the  point  closest  to  the  origin.  Similarly 
lu(t)l  =  Ik'x(t)  I  =  Ik'  (A-b  k'  )x(t)  I  and  ;k*  (A-b  k1 ) ,  can  be  used  as  a  measure 
for  lul. 
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4.6.2.  Short  Wordlength  Control  Law  Implementation 

The  feedback  control  law  may  be  implemented  approximately  in  a  short 
wordlength  microprocessor  as 


u  +  Au  =  (k1  +  4k' )  (x  +  Ax)  k'x  +  Ak*x  +  k' Ax. 


(84) 


For  small  x  the  dominant  term  in  Au  is  k'Ax,  i.e.  the  gains  should  be  not  too 
high.  For  large  x  the  dominant  term  is  Ak'x.  Robustness  with  respect  to  Ak' 
is  achieved  by  maintaining  a  distance  Ak^  from  a  boundary  in  each  direction 


Example  4.9: 


0  -4 

6/16 

x(k+l)  =  Ax(k)  +  bu(k),  A  = 

b  = 

! 

<*• 

T-* 

j*5/16 

u  =  -[k^  +  Ak  k2+Ak] 

Open  loop  eigenvalues  2=^‘  find  k^.k^  such  that  stability  is  achieved  for 
the  following  cases 


a) 

maximum 

admissible  Ak 

b) 

Ak  =  1, 

c) 

failure 

of  sensor  1, 

d) 

failure 

of  sensor  2, 

e) 

:k-'  =  k^ 

2 

+ k9  minimal. 

k '  =  £*E  =  [p_ 


1] 


5  6 

6  4 

4  -8  ; 


The  vertices  of  the  stability  triangle  ACC  in  the  k^-k^  plane  of  .g.  4.8  are 

determined  by  the  3  pole  placements  of  eq.  (52):  k'=[21  6],  k'  = [-1  -14], 

A  o 


a)  This  is  the  center  of  the  largest  square  inside  ABC  with 


sides  parallel  to  the  axis.  Here  k'  = [-0.4545  -10.7272]. 

It  admits  Ak = 1.4545  and  places  the  eigenvalues  at  z^ =0.132, 

=  0 • 686. 

b)  The  region  for  which  stability  is  robust  with  respect  to 

Ak  =  l,  is  the  triangle  A'B'C'  with  sides  parallel  to  those 

of  ABC,  with  a  distance  of  the  sides  of  a/2  under  +45°.  Note 

that  this  region  does  not  include  the  deadbeat  solution  at 
2 

k'  =e'A  =  [4  -8],  This  is  a  warning  that  points  with  the 

maximum  distance  from  the  stability  boundary  in  the  X-plane 
need  not  be  particularly  robust. 

c)  The  region  c  is  robust  with  respect  to  a  failure  of  sensor  1. 

d)  The  region  d  is  robust  with  respect  to  a  failure  of  sensor  2. 

e)  Point  E  with  k^  = [48/13  -72/13]  has  the  minimum  distance 

from  the  origin.  This  minimum  norm  solution  puts  the  eigen¬ 
values  at  s^  2  = 0.442+J0.897  on  the  unit  circle.  For  stability 
k^  can  be  increased  or  k2  decreased  by  an  arbitrarily  small  e. 

4.7.  Multi-Input  Problems 

4,7.1.  Characteristic  Parameter  Assignment 

In  the  single-input  case  it  was  convenient  to  have  the  canonical 


parameter  space  9  as  an  intermediate  step  between  the  \-plane  and  /(-space,  I 
aLlowed  studying  the  shape  of  stability  regions  without  reference  to  a 
particular  plant.  By  the  linearity  of  k'  =  jr-'-E  linear  properties  in  9  .pace 
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were  preserved  in  the  X  space.  X  and  ^  space  had  the  same  dimension  n  and 
the  mapping  was  one  to  one.  In  the  case  of  p  inputs  the  feedback  gain 
matrix  K  has  pxn  free  parameters.  Thus  a  pxn  dimensional  parameter  space 
with  the  elements  of  K  as  coordinates  is  naturally  defined.  The  question  here 
is  whether  there  exists  a  pxn  dimensional  canonical  parameter  space  ^  , 
which  is  linearly  related  to  X  hut  independent  of  the  particular  plant.  Another 
question  is:  If  it  exists,  is  there  a  simple  relationship  with  the  n-dimensional 
■P  space,  in  which  stability  or  nice  stability  is  defined  as  before.  The 
answer  to  both  questions  is  a  conditional  yes.  Some  results  are  available  in 
[56],  which  will  be  reformulated  and  used  to  design  simple  systems  for  robust¬ 
ness  with  respect  to  actuator  failures.  Before  the  main  result  can  be  formu¬ 
lated  we  have  to  make  some  additional  assumptions  and  to  introduce  some  notation. 

In  the  single-input  case  the  implicit  assumption  was  that  changes 
in  physical  parameters  do  not  cause  changes  in  the  order  n  of  the  system 
(which  is  true  for  all  cranes,  aircraft  ,  and  other  examples).  In  the  generali¬ 
zation  to  the  multivariable  case  it  is  assumed  here  that  the  controllability 
structure,  as  defined  by  the  Kronecker  indices,  is  unchanged  by  physical 
parameter  variations.  For  a  controllable  system 

x  =  A  x  +  Bu,  B  =  [b^  b^ . .  .bp]  (85) 

the  Kronecker  indices  -  ,  i=l,2,...,p,  are  the  smallest  integers  k  such  that 

r .  .  =  r .  ,  ,  .  where 
i ,  k  i  -  1 ,  k  ’ 

r;  =  rank  [B,A3...A^  ^3,A^b  . . .  A%  .  ]  i  =  l,2,...,p 

L  ,  K  XL 

r  ,  =  r  .  . 

o,  k  p , k-  1 


(86) 


The  Kronecker  indices  satisfy  ^  +  •  •  •  +  H-p  =  n*  A  vector  A 
regular  if  r^  ic  =  1-+ri_l  i-e.  the  regular  vectors 


is  called 


^  l”l 

b^jA  ••••A  2 

b2  jA  b9  ...  A  b^ 


b  ,Ab 
~P - P 


A  P  b 


are  the  first  n  linear  independent  columns  of  [B,A_B,A  B  . . .]  .  By  the  defini- 
tion  (86)  A  b  is  not  regular  and  can  be  uniquely  expressed  as 

A  V  =  -[B,A  B ...  A  B]«  -[A  ^  ...  A  (88) 

where  all  elements  of  a  and  multiplying  nonregular  vectors  are  zero,  in 
particular  jL  =  [3.^  . . .  3^.^]  'has  3^  =  0  if  By  Popov's  theorem  on 

invariants  [10]  the  Kronecker  indices  h- .  and  the  ^-parameters  3^,  i=l,2,...,p 
constitute  a  complete  set  of  independent  invariants  for  (A,B)  under  all 
transformations 


u  =  M  u 


In  order  to  avoid  the  distinction  between  a-  and  3-parameters  in  the 
definition  of  a  canonical  parameter  space  ,  first  the  3-parameters  are 
made  zero  by  an  input  transformation,  i.e.  a  modified  system 


x=Ax+BMu=Ax+  3  u 


•nsidered  with  the  "normalized  input" 


u  =  -K  x  4-  r 


re  laced  to  the  original  variables  by 
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B  K  =  B  K  ,  B  =  BM  ,  K  =  M-1K  (93) 


Now  the  g-parameters  of  the  pair  (A,B)  can  be  made  zero  by  the  choice 


Then  with  B  =  B  M =  [b,  . . .  b  ] 

P,  A.  -1 

A  ‘b  =  - [BjA  B  ...  A  1  B]^ 

A  bL  =  -[BM~  ,  A  BM~  .  . .  A  1  BM  ]a.  (95) 

(-  .  _  i-1  - 1  __ 

A  ibi  =  - [ B , A  13  ...  A  1  BjaL. 

A  comparison  with  eq .  (88)  shows  that  now  the  P-parameters  multiplying 

f*i- 

A  . . .  A  b^  ,  do  not  appear  in  eq .  (95). 

By  Popov's  theorem  the  P-parameters  remain  zero  under  feedback, 
i.e.  the  closed  loop  system  with  input  £=M  ^£ 


x  =  (A -B  K)x+Br=  (A  -  BK)x  +  B  r  =  Fx  +  B  r 


has  no  ^-parameters  and  F  1fc^  can  be  expressed  as 


F  bi  =  -[B,FB...F  B]£i. 


The  n  elements  of  p^,  i=l,2,...,p,  i.e.  nxp  parameters  will  be  used  in  the 
following  as  coordinates  of  a  canonical  parameter  space  6^,  in  which  desired 
closed  system  properties  can  be  specified  without  reference  to  the  particular 
plant,  with  the  only  condition  that  the  closed  loop  must  have  the  same 
Kronecker  indices  u,.  as  the  open  loop.  In  the  single-input  case  p.  =  p,  i.e. 
p^  consists  of  the  n  coefficients  of  the  desired  characteristic  polynomial, 
in  the  multi-input  case  the  nxp  elements  of  will  be  called  charac¬ 

teristic  parameters.  Note  that 


KB  KAB  KAk_1B 


_  ~  v  ~  k.'"' 

[B,  A  B  .  .  .  A  B]  =  [B,  F  B  . . .  F  B] 


Both  sides  may  be  truncated  at  any  column  of  AKB,  always  the  second  factor  on 
the  rigut  hand  side  is  nonsingular  and  thus  F^b  i  is  regular  if  and  only  if 
A^b^  is  regular.  Due  to  the  particular  form  of  M  in  eq.  (94)  the  same 
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Thus  in  eq.  (97)  the  regular  vectors  are  in  the  same  locations  as  in  eq.  (96), 
i.e.  the  list  of  regular  vectors  is  the  same  as  in  eq.  (87)  with  A  replaced 
by  F  = A-B K  and  replaced  by  fc4  Then  eq.  (97)  may  be  rewritten  as 


r.  r  .  -  j. 

F  6  =  -[b  F  b  ...  F  1  b*b  ...  :: 
*  ^  I  A.  •  ’  c.  •  • 


F  **  b  ]  *i2 

“P 


(100) 


where  p  .  =  [p . p..  ]  '  is  a  p,.-  vector.  Now  define  the  pxn  coefficient 

l-J  1-JO  l-JP'j”*-  J 

matrix 


4  = 


(101) 


It  generalizes  =  [PQ  . . .  p  to  the  multi-input  case.  The  generalization 

of  p* = [p  ...p  ,  1]  is  the  "characteristic  matrix" 

c  Lro  n-1 


0 


4  = 


P21  0  K2 


0  4 


piP  0 

P2P  ° 


JE  1 
^PP  j 


(102) 


Now  eq.  (100)  for  i=l,2...p  may  be  written 


"l"1- 


tl  bL  ...  F  Pbp]  =-[bL  ...  F  bL  I  b. 


F"P  Lb  ]P' 

t  »  »  t  r  p  n. 


(103) 


which  also  can  be  expressed  as 


L 


78 


[bn  . . .  F  1b-  I  b,  ...  I  I  ...  F  Pb  JP*'  =  O.  (104) 

1  1  •  4  •  •  “"p  “ft 

In  the  single-input  case  the  characteristic  polynomial  is  obtained  from  eq. 
(104)  by  replacing  the  vector  F in  eq.  (104)  by  X  .  In  the  same  way  in 
the  multi-input  case  a  "characteristic  polynomial  matrix"  is  defined  as 

2X  =  diag(Y^).P*\  V’  =  [l  \  ...  X  X].  (105) 

It  is  related  to  the  characteristic  polynomial  of  F  by 

P(X)  =  det  P^  (106) 


as  was  shown  in  [56].  The  main  advantage  over  a  direct  calculation  of  the 
det(Xl_-A +  B  K)  with  general  elements  k_  is  that  eq.  (106)  does  not  involve 
A  and  B  and  is  done  only  once  for  all  systems  with  the  same  Kronecker  indices. 
The  multilinear  problem  of  solving  eq.  (106)  for  some  characteristic  coefficient 
becomes  a  linear  one  if  the  n  coefficients  of  one  row  of  P^  are  expressed  in 
terms  of  the  £^,  i=l,2,...,n  and  the  remaining  px(n-l)  free  parameters  of  the 
other  rows  of  P^.  We  are  now  able  to  formulate 
Theorem  5;  A-BK  has  the  characteristic  matrix  P1V  if 


k  =  mpJe 

where  M  is  given  by  eq.  (94)  and 


i 

e.  is  the  last  row  of  the  a.  xn  matrix  Q.  in 
— l  i  -H 


(107) 


(108) 
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M.  -1  r*  -1 

r  «  [b.  . . .  a  b.  :  b9 a  p  b  j . 

Proof:  Follows  from  [56]. 

Example  4.10: 


5  -1  2 

0 

1 

X  = 

-2  -2  6 

x  + 

1 

5 

-3  7 

1 

6 

u  =  A  x  +  B  u  . 


(109) 


(110) 


Find  an  output  feedback 


u  = 

kll 

0 

k13 

_ k21 

0 

k23 

if  it  exists,  which  places  poles  at 

=  ~1, 

s2=-2 

maximum  Ik.. I 
ij 

detfbj^  b2  A  bj  #  0,  i.e.  ^  =  2, 
For  H  we  need  the  3  parameter 


-3  and  minimizes 


=  1. 


the 


A  b0  =  ■ 

-6  . 

.A  b. 

-a,b ,  -  a„b„ 

21 - 1 

1_  ]_  9—2 

it  is 

— 

— i 

1 

-5  1 

P  =  -5 

21 

and  N 

= 

I 

l 

1  * 

« 

1 

0 

L_i 

The  controllability  vectors 

i 

e ^  and 

1 

—2 

are 

obtained  from 
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n 


to 


to 
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1 

0 


1  -1 

-1  1 


E  * 


r  • 

e. 

“l  1  -ll 

-1 

t 

e,A 

H 

O 

T— 1 
1 

— 1— 

e*A2 
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-1  -2  5 

— 1— 

1 

—2 

r— 1 

H 

1 

O 

e 'a 

6-1  1 

—2— 

_ 

=  [ 
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The  characteristic  matrix  has  the  form 


* 

P110 

Plll 

1 

P120 

0 

*A  = 

P210 

P2 1 1 

0 

P220 

1 

It  is  related  to  the  characteristic  polynomial  by 


2  3 

P(s)  =  PQ+  p +  P2S  +  s  =  det  p^  =  det 


P110  +  P111S+S  P120 


P210  +  P211S  P220+S 


Use  the  coefficients  of  the  second  row  as  free  parameters  Y’PU0- 
6  =  ~P2ii’  and  s  -  "P220* 


Po  YP120'£P110 


pl  =  6p12'£plll  +  p110 


(111) 


P2  Pm"6' 


With  eq.  (107)  the  output  feedback  condition  results  in 


=  MPT  = 

k  A  1  0 

k22  ll 


p110  ”  P120  2 


(112) 


e  -  Y  =  1. 


Substitute  into  eq.  (Ill) 


P„  =  (e-l)p1on-e(2  +  p  )  =  -2e  -  p. 


P1  6p120  “  ePlll  +  2 +P120  2  "  eplll  +  (6+1)p120 


P2  =  Pill"6- 

The  particular  choice  Pj^o =  ®  *-eaves  6  undetermined  and  Po  =  -2e,  p^  =  2-ep^p 

2  3 

p2=Plll~e’  and  wi-th  P(s)  =  (s+1)  (s+2)(s+3)  =  b+lls+6s  +s  ,  s  = -3,  Pj^  =  3, 
pU0  =  2,  Y  =-4. 

The  remaining  four  feedback  gains  are 


kll  k13 


k2l  k23 


=  MP*[Ttl  l3] 


P110"P111 


-Y-h5+6 


'P110+plll+5+p120 


y-6-e+l 


1  -5~f  f~-2  6 

I 

! 

0  1_[  10+5  -5 


-52  -  56  6  -  56 


10  +  6 


&  is  chosen  as  -3.8  in  order  to  minimize  maxik..l. 


In  principle  all  concepts  for  the  design  of  robust  control  systems 
in  X  space  carry  over  from  the  single-input  to  the  multi-input  case.  However, 
the  X  space  now  has  dimension  pxn,  such  that  graphical  methods  in  several 
two-dimensional  subspaces  require  many  iterations  and  can  become  feasible  only 
with  a  good  software  for  graphical  displays  and  interactive  design.  General 
results,  which  would  give  a  clearer  understanding  of  the  multivariable  X  space 
are  presently  not  available. 

Changes  of  physical  parameters  can  be  tackled  in  the  same  way  as 
in  the  single-input  case.  Sensor  failures  now  result  in  p  coefficients  of 
one  column  of  K  becoming  zero.  Actuator  failures  could  not  be  handled  at 
all  with  the  single-input  method.  Here  some  insight  into  the  geometry  of  the 
problem  and  possible  problem  formulations  can  be  gained  from  eq.  (107),  as 
will  be  illustrated  by 
Example  4.11: 

fi  i~|  j  i  i  i  u  (k7! 

x(k+l)  =  |  x(k)  +  j  j  .  |  .  (113) 

[o  ij  !__i  2_;  _u2(k)j 

The  state  feedback  u=Kx  may  have  three  configurations 


a)  nominal 
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b)  failure  of  actuator  1 


—2  = 


21 


22 


c)  failure  of  actuator  2 


*1  = 


11 

0 


a  2 
0 


Find  K  such  that  it  places  a  double  pole  at  z=0.4  and  the  eigenvalues  in  the 
two  failed  cases  are  in  the  smallest  possible  circle  in  the  family  of  circles 
rr  of  Fig.  4.1. 


In  eq.  (107)  M  =_I  and 


11 


*-■&-  I 


1^21 


1 

0 


12 


22 


0 

1 


2 

2 

-1 

-1 


“I 

1  i 
1 
0 


2Pu  +  2-Pl2 
j  2P2i  '  P22 


'Pll  +  1  +  P 


12 


’p21  +  P22 


PU+Z 


12 


P(z)  =  det  i 


21 


p22  +  z  ; 


2  +^Pll  +  P22^Z_r  ^P11P22  "  Pl2P2p  ' 


,2  2 


In  case  a)  p(a)=(z-0.4)  =  z  -0.8z  +  0. 16.  Use  a  =  “P2^>  5  =  -p99  as  free 

parameters,  then  =  (3-0.8  and  p^9  =  P(p)/a.  Eliminating  a  and  p  by  and 

22 


gives 
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(k9l  +  2k,,  +  1.4)" 


kll(k21,k22)  k2l  +  k22  +  1  2k2l  -  4k22  -  1.6 


(114) 


k12^k2l,k22^  kll^k2l,k22^  "  k21  "  2k22  "  3*2‘ 


These  two  equations  describe  a  two  dimensional  surface  in  the  four  dimensional 
X  space.  All  points  on  this  surface  give  the  desired  pole  placement. 

The  failed  cases  are  single-input  problems  for  which  a  circle  in 
z-plane  maps  into  a  triangle  in  the  subspaces  of  the  remaining  gains.  The 
problem  can  then  be  described  geometrically  as  follows:  For  a  given  circle 
radius  the  two  triangles  are  obtained.  Now  we  are  looking  for  a  point  in 
the  surface  (112)  which  has  projections  into  the  k^-k^-plane  and  the  k^-k.^- 
plane,  which  lie  in  the  respective  triangles.  To  check  whether  such  a  solution 
exists  we  could  find  the  regions  in  the  surface  (112),  which  have  such  pro¬ 
jections,  and  see  whether  the  two  regions  overlap.  It  is  more  convenient 
however,  to  use  one  of  the  planes,  say  the  k^-k^-plane ,  where  the  triangle 
is  one  of  the  regions  and  the  other  is  obtained  by  reflecting  the  triangle 
in  the  Ic^-k^-plane  at  the  surface  (112)  into  the  k^-kp -plane .  If  there 
exists  an  intersection,  then  the  radius  of  the  z-plane  circle  is  reduced  until 
the  set  of  admissible  solutions  shrinks  to  a  point. 

We  begin  with  r  =  0.5.  The  two  triangles  are  shown  in  Fig.  4.9a  and 
b.  The  vertices  D,E,F  are  now  mapped  by  eq.  (Ill)  into  D'jE'jF'  in  the 
k^-k^-plane .  The  sides  of  the  figure  D'E'F'  are  not  straight  lines,  since 
they  have  been  reflected  through  a  curved  surface.  However,  it  is  easy  to 
check  that  there  exists  a  set  of  admissible  solutions  between  F'  and  C,  thus  r 
may  be  reduced.  Figure  4.9c  shows  the  case  r  =  0.4,  where  no  solution  exists. 
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Going  to  r  =  0.41  as  shown  in  Fig.  4.9d  then  results  in 


the  solution 


K  = 


0.278 

0.020 


0.522 

0.190 


The  following  eigenvalue  locations  are  obtained 


nominal 

Zl,2 

=  0.4 

failure  of  actuator  1 

Zl,2 

=  0.8  on  circle 

0.41 

failure  of  actuator  2 

Z1 , 2 

=  0.6  +  0.344j  on  Tq  ^ 

both  failures 

*1,2 

=  1  open  loop  unstable 

4.8.  Conclusions 

The  design  of  control  systems  in  the  parameter  space  \  of  state 
feedback  gains  has  been  studied.  Conceptually  this  offers  the  following 
advantages : 


1.  Robustness  with  respect  to  large  parameter  variations  can 
be  achieved.  It  is  possible  to  design  the  feedback  such  as 
to  maximize  the  admissible  variations  in  known  directions. 

2.  Robustness  with  respect  to  sensor  and  actuator  failures  can 
be  achieved.  The  feedback  can  be  designed  such  that  in  the 
considered  failure  situations  at  least  some  emergency  specifi¬ 
cations  are  met  or  such  that  the  deviation  from  the  nominal 
behavior  is  minimized. 

3.  The  feedback  can  be  designed  such  that  for  the  worst  case 
initial  conditions  the  maxima  I  required  control  input 
lu(t) i  is  minimized  and  thus  saturation  can  be  avoided. 
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4.  The  feedback  system  can  be  designed  for  robustness  with 
respect  to  short  wordlength  implementation  or  other 
inaccuracies  of  the  feedback  law. 

5.  Static  output  feedback  or  fixing  some  gains  simplifies  the 
analysis  since  it  reduces  the  number  of  free  parameters. 

It  will  of  course  give  less  favorable  results. 

6.  Dynamic  feedback  can  be  tackled  by  the  same  methods.  It  is 
particularly  desirable  in  situations  with  sensor  failures. 

These  conceptual  advantages  are  apparent  in  situations  with  only  two  essential 
parameters.  Here  it  is  general  engineering  practice  to  present  and  analyze 
results  in  diagrams,  showing  boundaries,  regions  and  their  overlap,  etc. 
However,  in  most  cases  the  boundary  points  have  to  be  calculated  point  by 
point  in  more  or  less  involved  computations.  In  this  paper  desirable  dynamic 
properties  of  control  systems  are  specified  in  terms  of  regions  in  the  eigen¬ 
value  (X.)  plane.  A  particularly  simple  pole  placement  algorithm  is  introduced 
and  used  for  mapping  boundaries  from  the  \-plane  to  the  X  space.  The  mapping 
of  boundaries  point  by  point  becomes  very  simple.  In  fact  all  examples  of 
this  paper  were  done  by  pocket  calculator.  This  makes  it  promising  to  develop 
software  for  computer-aided  design  with  rapidly  changing  graphical  displays 
of  boundaries  in  various  subspaces. 

Due  to  the  simplicity  of  the  mapping  it  was  possible  to  obtain  a 
few  general  results  on  the  shape  of  boundaries  and  stability  regions  in 
space.  For  the  design  of  digital  control  systems  a  family  of  circular 
boundaries  in  the  z-planc  can  be  used  to  characterize  desirable  dynamic 
properties.  Circles  have  particularly  nice  mapping  properties  and  the  convex 
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hull  of  their  image  in  the  X  space  is  easily  determined.  For  arbitrary  order 
systems,  it  is  useful  as  a  necessary  condition  for  the  existence  of  various 
robust  solutions.  Further  research  will  be  necessary  to  obtain  more  general 
results  on  the  shape  of  boundaries,  sufficient  conditions,  etc.  also  for 
other  than  circular  boundaries  in  X-plane.  The  development  of  good  numerical 
algorithms  will  depend  on  such  insights  into  the  geometry  of  the  solution  sets. 

It  was  shown  that  already  in  very  simple  examples  this  set  may  be  disconnected. 

A  systematic  search  inside  the  convex  hull  of  the  X  space  region  may  be 
necessary. 

In  its  present  stage  the  design  in  X  space  is  already  a  useful  design 
tool.  It  may  be  used  for  example  in  conjunction  with  the  root  locus  method, 
which  visualizes  the  influence  of  one  gain  on  the  eigenvalue  location.  The 
present  method  visualizes  the  influence  of  two  gains  on  the  eigenvalue  loca¬ 
tions.  The  use  ot  this  tool  has  been  shown  in  this  paper  by  the  example  of  a 
crane.  In  Chapter  5,  it  will  be  used  to  design  a  dynamic  controller  for  the 
short  period  longitudinal  mode  of  an  F4-E  aircraft  with  canards,  which  is 
unstable  in  the  subsonic  flight  conditions.  A  solution  using  two  gyros  and 
one  accelerometer  was  found  which  meets  the  nominal  specifications  for  the 
unfailed  system  or  after  a  failure  of  any  single  sensor,  and  also  meets  the 
emergency  specifications  after  failure  of  any  two  sensors,  where  these  properties 
pertain  to  four  very  different  flight  conditions. 


CHAPTER  5 


APPLICATION  TO  THE  LONGITUDINAL  CONTROL  OF  A  FIGHTER  AIRCRAFT 

5.1.  Introduction 

The  purpose  of  this  chapter  is  to  apply  the  parameter  or  K  space 
design  procedure  described  in  Chapter  4  to  a  realistic  design  problem.  The 
system  to  be  considered  is  a  third  order  model  of  the  longitudinal  axis  of  a 
McDonnell-Douglas  F4-E  fitted  with  horizontal  canards. 


5.2.  System  Description  and  Design  Objectives 
5.2.1.  Model  formulation 

Airframe  dynamics.  The  example  chosen  to  illustrate  K-space  design 
is  control  of  the  longitudinal  axis  of  a  fighter  aircraft.  The  complete 
equations  of  motion  describing  the  dynamics  of  the  airplane  are  nonlinear 
and  too  complex  to  be  used  in  control  law  development.  Standard  procedure 
is  to  linearize  these  equations  about  typical  flight  conditions,  and  then  use 
these  linear  system  representations  to  design  the  control  system.  The 
linearization  decouples  the  dynamical  equations  into  two  separate  sets  of 
equations  called  the  longitudinal,  and  lateral-directional.  Typically  a 
separate  control  system  is  designed  for  each  set  of  equations.  For  fighter 
aircraft  two  dominant  modes  describe  the  longitudinal  rigid  body  motion.  One 
of  these  modes,  called  the  phugoid  mode,  is  usually  very  slow.  It  is  easily 
controlled  by  the  pilot,  and  therefore  is  not  included  in  the  control  law 
design.  The  other  mode,  called  the  short  period  mode,  is  the  mode  that  most 
effects  the  handling  qualities  of  the  aircraft,  and  is  the  only  airframe 


mode  considered  in  this  report. 


The  state  description  of  the  short  period  mode  depends  on  the 
individual  aircraft  and  flight  condition.  The  aircraft  chosen  for  this  paper 
was  a  special  F4-E  fitted  with  horizontal  canards.  Figure  5.1  shows  the 
F4-E  flight  envelope  and  those  flight  conditions  for  which  linearized 
aerodynamic  data  are  available. 


Mach  Number 


Figure  5.1.  Flight  envelope  and  operating  points  [57]. 

A  complete  description  of  the  aircraft  including  aero  data  is 
given  in  [57].  One  notable  feature  is  that  the  uncontrolled  short  period 
mode  is  unstable  for  all  subsonic  flight  conditions.  Table  5.1a  lists  the 
uncontrolled  short  period  eigenvalues  for  each  of  the  four  flight 


conditions. 


FC 

Mach 

Altitude 

Open  Loop  Short 

Period  Eigenvalues 

1 

.5 

5000* 

-3.07 

1.23 

2 

00 

• 

5000* 

-4.90 

1.78 

3 

.9 

35000' 

-1.87 

.56 

4 

1.5 

35000' 

-.87  +  J4 

.3 

Actuators .  The  two  major  control  surfaces  available  for  control 

of  the  short  period  mode  are  the  elevator  position  (5  )  and  the  canard 

e 

position  (6c)  ;  these  are  shown  in  Figure  5.2.  Simplified  models  for  the 
canard  and  elevator  actuators  were  used.  The  actuator  state  equation  used  was 

S  =  -a<5  +  aS  (115) 

com 

where  6  is  the  actuator  position,  a  is  the  actuator  eigenvalue  (time  constant 


or  equivalent  bandwidth)  and  <$CC)m  is  the  commanded  actuator  position. 

Sensors .  Typical  sensors  which  are  used  to  control  the  pitch  axis 
are  (for  definitions  see  Figure  5.2): 

1)  Inertial  sensors  which  measure  pitch  rate  (q)  and  normal 
acceleration  (Nz). 

2)  Air  data  sensors  which  measure  angle  of  attack  (a)  and  dynamic 
pressure. 

3)  Position  sensors  which  measure  the  elevator  and  canard  positions. 
Only  the  pitch  rate  and  normal  acceleration  are  assumed  to  be  available 
for  feedback.  Air  data  sensors  were  not  used  due  to  their  unreliability. 
Position  sensors,  although  reliable,  would  not  be  useful  without  an  estimate 


V=  velocity  vector 
We  weight 
ci=  angle  of  attack 
Gv  pitch  angle 
q*=G=  pitch  rate 

N'zHnormal  acceleration  at  sensor  (*) 


Figure  5.2.  Aircraft  nomenclature 
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of  the  equilibrium  (trim)  surface  position.  This  estimation  was 
undesirable.  A  failed  sensor  will  be  one  which  in  the  output  is  uncorrelated 
with  the  input. 

System  state  description.  The  aero  data  and  equations  of  motion 
given  in  [57]  were  transformed  to  a  state  space  description  resulting  in  the 
following  system  representation  for  the  airframe  and  actuators 


al 

a2 

a3 

a 

a5 

a6 

a7 

a 

0 

0 

-a 

e 

0 

0 

0 

0 

-a 

f 

C1 

C2 

C3 

c 

t 

0 

1 

0 

0 

q  0 

+ 

6  a 


(116a) 


(116b) 


where  a^~ag  and  c^-c^  depend  on  the  flight  condition  and  are  listed  in 
Appendix  I. 

This  two-input  representation  was  reduced  to  a  single  input  problem 
by  considering  the  canard  command  to  be  proportional  to  the  elevator  command. 
A  study  was  done  in  [57]  to  determine  this  proportionality  constant,  Kc,  for 
minimum  drag  flight  under  a  wide  range  of  conditions.  That  minimum  drag 
value,  Kc  =  -.7,  was  the  value  used  in  this  paper.  There  is  not  much  loss  of 
generality  in  this  assumption  of  dependent  inputs  since  the  ratio  of  control 
surface  effectiveness  (essentially  aD/a_  since  a_  and  a.  are  small)  between 
the  canard  and  elevator  does  not  change  much  as  flight  conditions  change. 

To  study  the  effect  of  sensor  failures  it  is  easiest  to  have  the 
system  equations  written  in  sensor  coordinates.  Equations  (116)  were 


4 


transformed  to  make  the  normal  acceleration  a  state  variable.  Transforming 


the  system,  using  K  “-.7,  and  assuming  a  *  a  resulted  in  the  following 

C  C  6 

system  description. 
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(117a) 


with  output 


(117b) 


where  <$e  now  represents  the  effective  control  position  of  the  elevator  and 
canard,  and  "a"  represents  the  effective  bandwidth  of  both  actuators.  This 
effective  bandwidth  was  assumed  to  be  14  rad/sec.  Appendix  I  gives  a  tabula¬ 
tion  of  the  a^  ,  b^  used. 

5.2.2.  Design  Objectives  and  Design  Specifications 

Using  normal  acceleration  and/or  pitch  rate  feedback,  the  basic 
design  objective  is  to  design  a  continuous  time,  fixed  gain  controller  which: 
1)  meets  certain  nominal  performance  requirements  at  all  four  flight  condi¬ 
tions  when  all  sensors  are  available  and  2)  meets  emergency  performance 
requirements  after  sensor  failure.  The  specific  requirements  to  be  met  are: 

1)  the  controlled  short  period  eigenvalues  must  be  in  the  range 
specified  by  military  standards  [58], 

2)  the  remaining  closed  loop  eigenvalues  resulting  from  actuator  and 


controller  dynamics  should  be  in  a  specified  region. 
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3)  time  responses  of  the  so  called  output  to  pilot  step  commands 
should  be  acceptable  in  the  sense  of  [57]. 

The  region  of  allowable  short  period  eigenvalue  location  is  given  in  [58] 
as  a  requirement  on  the  range  of  damping  and  natural  frequency  for  the  short 
period  mode.  For  the  short  period  mode  described  by 


2  2 
S+2£<dS+0J  =0 

sp  sp  sp 


(ns; 


the  restricted  range  of  Csp  and  ui  under  normal  operating  conditions  is 


.35  <  ?  <1.3 

sp  " 


0)  <  0)  <  0i, 

a  sp  b 


(119a; 

(119b; 


and  for  emergency  conditions  is 


.15  <  5 


sp 


oj  <  w 
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(120b; 


where  u>  ,  ai  ,  and  u  depend  on  flight  condition. 

3  D  C 

Table  5.2  lists  the  frequency  range  for  each  flight  condition  and 

Figure  5.3  shows  the  region  defined  by  equations  (119).  Since  the  short 

period  damping  can  be  greater  than  1,  a  single  real  eigenvalue  is  permitted 

outside  the  circle  of  radius  w,  or  inside  the  circle  of  radius  w  .  As 

b  a 

discussed  next,  the  single  real  eigenvalue  region  S > u  would  overlap  with 

b 

the  permissible  region  for  non-short  period  eigenvalues  and  thus  make  a 
distinction  between  the  two  types  of  eigenvalues  impossible.  To  prevent 
this  overlap,  the  simplified  region  shown  in  Figure  5.3  will  be  used  as  the 
short  period  eigenvalue  permissible  region.  Since  all  real  pairs  of  eigen¬ 
values  inside  the  simplified  region  result  in  a  damping  less  than  1.3,  all 
eigenvalues  in  the  simplified  region  meet  military  specifications. 


Figure  5.3.  Allowable  short  period  eigenvalue  locations. 


As  described  in  [57]  this  aircraft  has  several  lightly  damped 


structural  modes  which  were  not  modeled  in  this  paper.  The  control  bandwidth 
should  be  less  than  the  lowest  structural  frequency,  which  is  85  rad/sec,  so 
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the  upper  limit  on  all  eigenvalues  was  set  at  70  rad/sec.  This  is  a  high 
limit,  and  in  the  design  process  an  attempt  will  be  made  to  lower  this  value. 
The  requirements  on  the  non-short  period  eigenvalues  (requirement  2)  will  be 
that  they  lie  in  the  region  defined  by 

to.  <  a)  <  70  rad/sec  (121) 

D 

with  a  minimum  damping  of  .35  (see  Figure  5.4).  In  the  emergency  situations 
no  distinction  between  eigenvalues  will  be  made.  Figure  5.5  shows  the 
required  region  for  all  eigenvalues  under  emergency  conditions. 

Well  placed  eigenvalues  do  not  guarantee  good  time  responses. 
Requirement  3  will  ensure  well  behaved  transient  response.  The  response  of 
most  interest  will  be  the  C*  response.  As  discussed  in  [57],  C*  is  a  linear 
combination  of  normal  accleration  and  pitch  rate  given  by 

C*  =  (Nzp  +  12.43q)/Kc*  (122) 

where  Nzp  is  the  normal  acceleration  at  the  pilot's  location  (same  as  Nz 
for  this  plane)  and  the  stationary  value,  Kc*,  is  used  for  normalization.  The 
C*  response  to  a  pilot  step  command  should  fall  in  the  region  shown  in 
Figure  5.6. 


5.3.  Design  Using  Static  Output  Feedback 

5.3.1.  Robustness  with  Respect  to  Changing  Flight  Conditions 

The  first  design  objective  will  be  to  design  a  controller  which 
meets  the  nominal  requirements  at  all  flight  conditions.  We  initially  assume 
a  controller  structure  using  only  static,  or  unEiltered,  feedback  of  the  two 
states;  normal  acceleration,  and  pitch  rate.  Figure  5.7  shows  the  system 


Time  ■  Seconds 


I 

Figure  5.6. 


C  response  envelope 


structure  where  equations  describing  the  F4-E  dynamics  were  given  in 
Section  5.2,  and  the  pre-filter  used  to  shape  the  step  responses  is  the 
same  as  in  [57]. 

Computer  programs  using  the  algorithms  of  Chapter  4  were  developed 

to  perform  the  mapping  of  the  eigenvalue  constraints  given  in  Section  5.2. 

Boundaries  in  the  k^z,  k^  plane  were  found  for  each  flight  condition.  A 

typical  boundary  (for  flight  condition  2),  is  shown  in  Figure  5.8a.  For 

gains  in  the  region  Rnom2  the  closed  loop  eigenvalues  will  all  be  in  the 

region  rnom2*  which  is  shown  in  Figure  5.8b.  Each  section  of  the  K-space 

boundary  is  described  in  Figure  5.8b. 

The  regions  ^nom2-^nom4  were  f°und  by  mapping  the  eigenvalue 

constraints  for  each  flight  condition  (T  -T  .).  The  intersection  of 

°  noml  nom4 

these  regions,  R  ,  is  shown  in  Figure  5.9.  For  any  gain  chosen  in  R  the 
nora  °  b  non 

ith  flight  condition  will  have  closed  loop  eigenvalues  in  I\,  i  =  l,2,3,4. 
Therefore,  the  requirement  that  the  system  be  robust  with  respect  to  changing 
flight  condition  can  be  met  by  using  static  output  feedback.  Since  Rnom  does 
not  intersect  either  axis,  no  robustness  with  respect  to  either  sensor  failure 
can  be  achieved  by  static  output  feedback. 

5.3.2.  Selection  of  a  Gain  in  the  Permissible  Region 

Any  gain  choice  from  R  would  meet  the  nominal  eigenvalue 

nom 

requirements.  Several  alternative  methods  are  available  to  aid  in  the 
selection  of  a  specific  design  point.  One  method  could  be  analysis  of  the 
eigenvalue  locations  corresponding  to  different  points  in  K-space  as  described 
in  [49].  Every  point  in  the  k^,  k^  plane  represents  a  set  of  closed  loop 
eigenvalues  for  each  flight  condition.  Using  the  mapping  algorithms  of 


F°r  kN2'  kq  °n 

a'-b'  the  short  period  eigenvalues  are  on  a-b  ( oj 2=3 . 

b'-c'  the  short  period  eigenvalues  are  on  b-c 

c'-d'  the  actuator  eigenvalue  is  at  d  (a 2~~' 

d'-e'  a  short  period  eigenvalue  is  at  e  (c2=-! 

e'-a'  the  actuator  eigenvalue  is  at  e  ( ^r2==— ' 

Figure  5.8b.  Region  in  s-plane  and  description  of  K-space  boundary. 


Chapter  4  families  of  constant  damping,  constant  frequency  and  single  real 
root  curves  may  be  obtained.  Figure  5.10  shows  such  families  for  flight 
condition  2. 

Another  method  would  be  to  decrease  the  allowable  eigenvalue  regions 

and  thereby  decrease  the  size  of  R  .  Figure  5.11  shows  the  region 

where  the  high  frequency  limit  has  been  lowered  from  70  to  50  rad/sec,  the 

minimum  damping  has  been  increased  to  .5,  and  the  minimum  short  period 

frequency  has  been  increased  by  fifty  percent  for  each  flight  condition.  Any 

gain  chosen  from  R^^  would  meet  these  tighter  requirements  at  all  four 
nom 

flight  conditions. 

A  further  technique  is  to  use  gains  which  will  require  smaller 
control  inputs.  In  [34]  Ackermann  shows  for  systems  where  the  state  variables 
have  been  normalized  to  their  maximum  values,  the  distance  from  the  origin  in 
K-space  can  be  used  as  a  measure  of  the  maximum  control  needed.  Where 

| u |  =  |k’x|  <  Ilk II  *11x11  .  (123) 

Even  though  Nz  and  q  have  not  been  normalized  to  their  maxima,  this  principle 
can  easily  be  demonstrated.  For  gains  g^,  g£,  g^  of  Figure  5.11,  the  6^  and 
C*  responses  to  a  step  command  for  flight  condition  1  are  shown  in  Figure  5.12. 
As  seen  from  the  figure,  the  control  (essentially  6^)  is  less  for  the  smaller 
gains,  as  is  the  control  rates. 

Using  the  above  methods  as  guides,  point  of  Figures  5.8-5.10 
was  selected  as  a  trial  design  point.  Using  these  trial  gains,  the  step 
responses  shown  in  Figures  5.13a,b  were  obtained.  The  responses  are  seen 
to  be  well  within  the  required  boundaries  for  all  flight  conditions. 

Table  5.3  lists  the  closed  loop  eigenvalue  locations  using  the  gains  of  point 


s 


Q^.  The  notation  (?>“n)  refers  to  the  damping  and  natural  frequency  of  the 
complex  pair  of  eigenvalues.  All  eigenvalues  are  well  within  the  limitations 
and  point  appears  to  be  a  good  choice  for  the  controller  gains . 


Table  5.3. 


Closed  Loop  Eigenvalues  k^zs*-,115 


k 

q 


-.8  (point  Q  ) 


F.C. 

Mach 

Altitude 

Closed  Loop  Eigenvalues 

1 

.5 

5000' 

(.94, 

4.68) 

-19.31 

2 

.85 

5000' 

(.61, 

9.18) 

-37.29 

3 

.9 

35000’ 

(.79, 

4.63) 

-17.78 

4 

1.5 

35000' 

(.55, 

8.11) 

-27.04 

5.3.3.  Robustness  with  Respect  to  Sensor  Failure 

An  additional  method  to  help  choose  a  final  gain  from  R  would 

nom 

be  to  consider  only  those  gains  which  will  meet  the  emergency  specifications 

after  sensor  failure.  To  determine  these  gains  the  emergency  eigenvalue 

regions  were  mapped  into  the  K^,  plane  for  each  flight  condition.  The 

intersection  of  these  four  regions,  Rgm,  is  shown  in  Figure  5.14  along  with 

the  region  R  redrawn  from  Figure  5.9. 
nom 

For  the  controller  structure  of  Figure  5.7,  failure  of  either 
sensor  results  in  the  corresponding  gain  going  to  zero.  Since  R  does  not 
intersect  either  axis,  there  is  no  gain  choice  which  is  robust  with  respect 
to  failure  of  either  sensor.  Additional  sensors  or  dynamic  feedback  are 
therefore  needed  to  meet  the  requirement  that  the  system  meet  the  emergency 
eigenvalue  requirements  after  sensor  failure. 


For  a  controller  with  m  identical  sensors  with  summed  outputs,  a 


sensor  failure  reduces  the  total  feedback  gain  by  only  a  fraction  ((m-l)/m). 

For  the  controller  structure  shown  in  Figure  5.15  (m=2) ,  any  k.,  ,  k  chosen 

Nz  q 

from  the  region  R  of  Figure  5.16  will  for  all  flight  conditions: 

1)  meet  the  emergency  requirements  after  failure  of  any  single  sensor 

2)  meet  the  emergency  requirements  after  the  combined  failure  of  an 

accelerometer  and  gyro,  and 

3)  meet  the  nominal  eigenvalue  requirements  when  no  sensor  has  failed 

After  an  accelerometer  failure  the  Nz  feedback  is  havled.  For  design  point 
Q^,  shown  again  in  Figure  5.16,  this  failure  would  result  in  point  A. 

Similarly  a  single  gyro  failure  would  result  in  point  B.  A  combined  failure 

of  an  accelerometer  and  a  gyro  would  result  in  point  C. 


5.4.  Design  Using  Dynamic  Output  Feedback 
5.4.1.  Search  for  Filters 

As  shown  in  Section  5.3,  a  controller  using  static  output  feedback 
requires  two  sensors  to  meet  the  nominal  specifications,  and  four  sensors  to 
meet  the  emergency  specifications  after  failure  of  any  single  sensor.  The  us 
of  dynamic  feedback  of  the  outputs  may  reduce  the  minimum  number  of  sensors. 
For  example,  if  an  estimate  of  the  normal  acceleration  was  obtained  from  the 
gyro  signal,  the  accelerometer  may  no  longer  be  necessary.  A  preliminary 
problem  is  to  find  a  candidate  dynamic  feedback  structure. 

At  first  it  may  seem  all  that  need  be  done  is  to  use  a  Luenberger 
observer  to  construct  these  estimates.  Figure  5.17  shows  how  an  observer 
might  be  implemented.  Now,  when  the  pitch  gyro  fails,  and  the  dashed  branch 


A 


Figure  5.17.  Possible  observer  structure. 

is  no  longer  available,  the  estimate  q  still  remains.  There  are  two  serious 
problems  with  this  method.  First,  the  system  dynamics  change  dramatically 
with  flight  condition  and  finding  a  non-adaptive  observer  that  gives  reason¬ 
able  estimates  of  q  and  Nz  would  be  difficult,  possibly  impossible,  to  find. 

As  seen  from  the  data  in  Appendix  I  the  control  effectiveness,  essentially  the 
open  loop  gain,  changes  by  a  factor  of  three  with  changing  flight  condition, 
and  the  open  loop  eigenvalues  vary  as  shown  in  Table  5.1. 

The  second  undesirable  feature  is  when  one  sensor  fails,  the 
observer  connected  to  it  still  is  driven  by  the  control  u.  Therefore  the 
observer  still  has  an  effect  on  the  closed  loop  characteristic  equation.  Even 


for  adaptive  observers  which  closely  match  the  system,  this  effect  would 
complicate  the  design  since  the  separation  principle  would  no  longer  hold 
after  sensor  failure.  For  these  reasons  observer  structures  were  not 


considered . 

In  [48]  G.  Kreisselmeier  discusses  the  use  of  inverse  filters  for 
robust  control,  which  are  applicable  to  minimum  phase  systems.  Figure  5.18 
shows  the  structure  of  the  controller. 
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Figure  5.18.  Structure  of  the  inverse  filter. 


If  the  filter  elements  and  are  chosen  equal  to  and  N0,  the  variabl 
V.z,  q  would  be  estimates  of  Nz  and  q.  For  this  case  by  choosing 


the  feedback  is  the  same  as  that  shown  in  Figure  5.15.  Since  and  change 
with  flight  condition  a  non-adaptive  controller  will  not  generate  true 
estimates  of  Nz  and  q,  but  only  signals  which  hopefully  are  closely  enough 
related  to  Nz  and  q  in  magnitude  and  phase  to  help.  If  as  flight  condition 
changes,  K^/K^  and  the  roots  of  and  N^  remain  close,  the  inverse  filter 

idea  may  succeed.  Table  5.4  lists  the  critical  values. 


Table  5.4. 

Open  Loop 

Zeros  and  Gain 

Ratio 

Mach 

Altitude 

q-Zero 

Nz-Zeros 

Vkq 

.5 

5000’ 

-  .884 

-.542+j5. 33 

.527 

.85 

5000' 

-1.57 

- .929+j9 . 12 

.536 

.9 

35000' 

-  .637 

-.392+j5.67 

.537 

1.5 

35000' 

-  .826 

-  .481+j8 .05 

.577 

Averaged  values 

-  .98 

0.586  j 7 .04 

.543 

The  averaged  values  were  used  to  construct  two  filters 


_q_  =  1 . 84  (s  +  .98) 

Nz  s2  +  1.172s  +  49.9 


Filter  1 


Nz  =  s  +  1.172s  +  49.9 
q  1.84(s+  .98)<s+ 10) 


Filter  2 


where  the  pole  at  -10  was  added  to  make  Filter  2  realizable. 


The  poles  of  Filter  1  are  only  slightly  damped,  and  do  not  cancel 


well  with  the  system  zeros.  When  it  was  tried,  highly  oscillatory  step 
responses  resulted  and  this  filter  was  discarded.  When  Filter  2  was  imple¬ 
mented,  reasonable  step  responses  resulted  using  typical  gains  chosen  with  the 
aid  of  the  results  from  subsection  5.3.1.  This  filter  was  retained  for 
further  investigations.  The  pole  at  -.98  is  weakly  controllable  since 
essentially  it  is  cancelled  with  the  pitch  rate  zero.  The  system  eigenvalue 
resulting  from  this  pole  will  be  exempt  from  the  eigenvalue  requirements  gives 
in  Section  5.2.  The  system  time  responses  will  give  a  measure  as  to  the 
validity  of  this. 

The  system  of  Figure  5.19  was  used  to  study  different  types  of 
dynamic  feedback  of  Nz.  Several  filters  were  tried.  One  that  appeared 
promising  was 


Figure  5.19.  Acceleration  feedback  only. 
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Nz 


s 

s  +  15 


(126) 


As  discussed  in  [59]  Nz  is  a  close  approximation  to  a  which,  if  used  as  a 
feedback  signal,  can  increase  the  short  period  damping.  For  low  frequencies 
this  filter  acts  as  a  differentiator. 

Figure  5.20  shows  the  structure  of  the  controller  to  be  investigated, 
which  may  be  represented  in  state  form  as 
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(127b) 


5.4.2.  Dynamic  Feedback  of  Pitch  Rate 

4  To  check  if  the  design  specifications  can  be  met  using  only  pitch 

rate  feedback,  the  system  of  Figure  5.20  was  used  with  and  both  zero. 

Each  nominal  (emergency)  eigenvalue  region  was  mapped  into  the  k^,  k^  plane 

•  The  intersection  of  the  four  regions,  R  (R  ) ,  is  shown  in  Figure  5.21. 

°  nomq  emq 

The  nominal  eigenvalue  specifications  can  be  met  by  choosing  any  gain  set 
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Figure  5.20.  Structure  of  controller  using  dynamic  output  feedback. 


from  Rnomq*  Similarly  the  emergency  specifications  can  be  met  by  use  of  any 

gains  in  R 
a  emq 

Gyro  failure  now  corresponds  to  a  simultaneous  reduction  of  k^  and 
k^.  For  a  system  using  m  parallel  gyros  with  identical  feedback  gains,  the 
failure  of  f  gyros  reduces  both  k£  and  k^  by  a  factor  (m-f)/m.  For  ra  =  2  and 
f=l,  any  gain  in  of  Figure  5.22  will  meet  the  nominal  specifications  for 
the  unfailed  system  and  also  meet  the  emergency  specifications  after  failure 
of  either  gyro.  Similarly  for  m=3,  f  =  1  or  2  any  gain  in  will  meet  the 
nominal  specifications  without  failure  and  still  meet  the  emergency  specifics 
tions  after  any  two  sensor  failures.  This  is  an  especially  nice  result  since 
reliability  criteria  often  stipulate  that  the  control  system  be  capable  of 


handling  a  failure  of  any  two  sensors.  Specific  eigenvalue  locations  and 
typical  C*  responses  are  shown  in  the  last  section  of  this  chapter. 

5.4.3.  Dynamic  Feedback  of  Normal  Acceleration 

The  system  of  Figure  5.19  was  studied  using  the  filter  of  equation 
(126).  The  nominal  eigenvalue  regions  were  mapped  into  the  k^,  plane. 

There  was  no  intersection  of  these  regions.  Thus,  for  the  controller  structure 
of  Figure  5.20  the  gyro  is  necessary.  The  intersection  of  the  emergency 
regions  RemNz  is  shown  in  Figure  5.23.  Any  gain  chosen  from  RemNz  will  meet 


the  emergency  specifications  at  all  flight  conditions. 
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5.4.4.  Dynamic  Feedback  of  Both  Outputs 

Now  the  total  system  of  Figure  5.20  will  be  considered  and  the 

complete  regions  in  K-space  will  be  four  dimensional.  First  the  system  using 

a  single  gyro  and  a  single  accelerometer  will  be  considered. 

For  the  system  to  be  robust  with  respect  to  either  sensor  failure, 

k0  and  k_  should  lie  in  the  region  R  of  Figure  5.24  (partially  shown  in 
2  5  emq 

Figure  5.21),  and  k^,  k^  should  lie  in  the  region  RemNz  of  Figure  5.23.  Sinc< 
graphical  representations  are  limited  to  two  dimensional  subspaces  of  K-space 
two  of  the  four  gains  will  have  to  be  at  least  temporarily  fixed.  As  k^  and 
k^  range  throughout  the  region  RemNz  the  closed  loop  eigenvalues  do  not  move 
much.  It  is  therefore  natural  to  select  k^  and  k^  as  the  gains  to  be  fixed. 
They  were  fixed  as 

(kr  k4)  =  (-.034,  -.1)  (128) 

designated  in  Figure  5.23.  There  are  two  free  gains  to  be  determined, 
and  k^.  The  four  dimensional  problem  has  been  reduced  to  looking  in  the  k^,  1 
plane  for  a  solution. 

As  before,  the  nominal  eigenvalue  boundaries,  V  .,  are  mapped 

°  nomi 

into  the  k^,  k^  plane  using  the  algorithms  of  Chapter  4.  A  typical  mapping 

is  shown  in  Appendix  II.  Figure  5.23  shows  the  intersection  of  all  four 

regions,  labeled  R  .  For  any  k_ ,  k„  chosen  from  R  of  Figure  5.24  the 

nominal  eigenvalue  requirements  will  be  met  at  all  flight  conditions.  Since 

k^,  k^  were  chosen  from  Rem^z  of  Figure  5.23,  the  emergency  requirements  are 

met  after  gyro  failure.  If  k„ ,  k„  are  also  chosen  from  R  the  system  will 

23  emq 

also  be  robust  with  respect  to  accelerometer  failure.  All  gains  in  the 
region  R  of  Figure  5.24  meet  the  above  requirements  for  k^  and  fixed  as  in 
equation  (123) . 


The  eigenvalue  requirements  for  the  nominal  and  emergency  conditions 
can  be  met  with  one  accelerometer  and  one  gyro,  or  with  two  gyros.  With  the 
assumed  filter  structure  the  gyro  is  necessary.  Step  responses  were  obtained 
for  several  robust  gain  combinations.  The  gains  chosen  were  based  on  the 
parameter  plane,  tightening  constraints  and  minimal  control  magnitude 
techniques  discusses  in  subsection  5.3.2.  The  responses  behaved  as  antici¬ 
pated  with  the  exception  of  some  responses  for  flight  condition  2.  For  this 
flight  condition  the  pitch  rate  zero  was  the  farthest  from  the  "cancelling" 
pole  in  the  Nz  filter  (see  Table  5a).  Figure  5.25  shows  typical  C*  responses 
at  flight  condition  2  using: 

1)  gains  from  Chapter  4,  k^z  =  -.115  k ^  -  -.8, 

2)  gains  k^,  k^,  k^,  k^  =  (-.034,  -1.5,  -1,  -.1), 

3)  gains  (0,  -1.5,  -1.0)  (no  accelerometer). 

The  filter  pole  at  -.98  was  exempted  from  the  eigenvalue  constraints  under  the 
assumption  that  it  would  have  little  effect  on  time  response.  As  seen  in 
Figure  5.25  the  more  acceleration  feedback  (and  less  use  of  Nz) ,  the  better 
this  assumption  is.  An  accelerometer  will  therefore  be  included  in  the  trial 
design  point. 

To  ensure  the  system  is  robust  with  respect  to  any  two  sensor 
failures,  two  gyros  will  also  be  included.  The  gains  chosen  were  those  repre¬ 
sented  above  by  point  shown  in  Figures  5. 5-5. 8.  The  final  configuration 
is  shown  in  Figure  5.26.  Figures  5.  a,b  show  the  response  to  pilot  step 

commands  for  the  unfailed  svstem.  The  C*  responses  are  well  within  the 

N 

required  boundaries  at  all  flight  conditions.  Table  5.5  lists  the  eigenvalue 
locations  for  the  unfailed  system  and  after  all  possible  combinations  of  sensor 


Flight 

Condition 


No 

Failures 


All  Sensors 
Failed 


One  Gyro 
Failed 


-.89 

1.23 

-.89 

(.64, 

4.85) 

-.98 

(.52, 

4.29) 

(.78, 

12.9) 

-10.0 

(.89, 

9.88) 

-27.6 

-14.0 

-31.76 

-15.0 

-1.43 

1.78 

-1.35 

(.84, 

5.66) 

-.98 

(.94, 

6.28) 

(-66, 

15.5) 

-10.0 

(.53, 

10.3) 

-48.13 

-14.0 

-55.5 

-15.0 

-.67 

.56 

-.70 

(-61, 

4.98) 

-.98 

(.49, 

4.4) 

(.79, 

12.3) 

-10.0 

(.9,  9 

>.76) 

-26.5 

-14.0 

-15.0 

-30.1 

-.88 

-.98 

-.91 

(.77, 

6.4) 

(.20,  4.4) 

(.91, 

7.84) 

(.66, 

13.2) 

-10.0 

(.45, 

8.33) 

Flight 

Accel. 

Accel,  and 

Two  Gyros 

-  '  Condition 

Failed 

Gyro  Failed 

Failed 

-.87 

-.84 

-.98 

(.65, 

4.38) 

(.59, 

3.07) 

(.80,  1.77) 

-15.0 

-15.0 

-5.78 

(.56, 

18.1) 

(.82, 

13.7) 

-10.0 

-35.33 

-1.62 

-1.81 

-.98 

(.72, 

4.62) 

(.66, 

3.83) 

-3.54 

-15.0 

-15.0 

(.27,  6.73) 

(.35, 

28.3) 

(.59, 

17.8) 

-10.0 

-61.43 

-.62 

.59 

-.98 

(.62, 

4.46) 

(.51, 

3.23) 

(.43,  2.34) 

-15.0 

-15.0 

-6.44 

(.60, 

16.9) 

(.85, 

13.2) 

-10.0 

-33.27 

-.86 

-.89 

-.98 

(.68, 

5.26) 

(.53, 

5.56) 

(.16,  6.27) 

-15.0 

-15.0 

-5.35 

(.42, 

22.3) 

(.69, 

14.4) 

-10.0 

n  (s,<ju  )  refers  to  the  damping  and  natural 
f  a  comple  pair  of  eigenvalues. 


failure.  All  eigenvalues,  of  course,  are  in  the  required  regions.  Figures 


5.28a,b,c  show  C*  responses  for  all  possible  sensor  failures.  In  these 
figures  the  normalizing  factor  (Kc*  of  equation  (122))  was  taken  to  be  the 
value  of  C*  at  3  seconds,  which  is  sufficiently  close  to  the  stationary  value 
for  the  closed  loop  responses  and  finite  for  the  unstable  open  loop  responses. 


5.5.  Summary  and  Conclusions 

A  fixed  gain  controller  has  been  designed  for  the  short  period  mode 
of  a  F4-E  aircraft  which  is  destabilized  by  horizontal  canards.  The 
uncontrolled  mode  is  unstable  at  all  subsonic  flight  conditions  and  insuf¬ 
ficiently  damped  in  supersonic  flight.  Therefore,  the  control  system  is- 
necessary  at  all  flight  conditions  in  order  to  achieve  handling  qualities  as 
specified  by  the  military  requirements.  The  problem  of  sensor  failures  is 
usually  solved  by  use  of  redundant  sensors  and  failure  detection.  In  this 
paper  controller  schemes  not  requiring  failure  detection  were  studied.  The 
control  system  was  designed  using  graphical  techniques  based  on  the  K-space 
methods  of  Chapter  4.  The  main  feature  of  this  graphical  design  method  is  the 
ability  to  determine  the  effect  of  two  controller  parameters  on  the  systo... 
eigenvalues.  First  the  unfailed  system  using  one  gyro  and  one  accelerometer 
was  studied.  This  is  a  third  order  system  with  the  two  free  parameters  being 
the  two  feedback  gains.  The  set  of  admissible  gains  was  determined  for  which 
the  nominal  requirements  are  met  at  all  flight  conditions.  No  gains  from  this 
set  are  robust  with  respect  to  either  sensor  failure.  If  both  sensors  are 
paralleled  by  an  identical  sensor,  the  system  can  be  made  robust  with  respect 


to  any  single  sensor  failure,  corresponding  to  a  fifty  percent  gain  recuction. 


In  order  to  reduce  the  number  of  sensors  needed  to  three  and  to 


I 
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S 
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achieve  robustness  for  double  sensor  failures,  various  dynamic  feedback  con¬ 
figurations  were  studied.  In  order  to  avoid  the  increased  dimensionality  of 
a  general  dynamic  feedback  configuration,  preliminary  studies  were  made  to 
find  good  candidates  for  two  filters,  one  for  each  input,  such  that  after 
failure  of  either  sensor  there  reamin  two  independent  variables  for  feedback. 
Using  the  most  promising  configuration,  it  turned  out  that  without  gyros  only 
the  emergency  specifications  could  be  met.  The  gains  for  the  accelerometer  and 
its  filter  were  thus  fixed  to  guarantee  this  property.  In  the  plane  of  the  two 
remaining  gains  for  the  gyro  and  its  filter,  the  set  of  admissible  solutions 
for  both  the  nominal  and  emergency  specifications  was  determined.  A  solution 
using  two  paralleled  gyros  and  one  accelerometer  was  found  which  met  the 
nominal  specifications  with  no  failure  or  a  failure  of  any  single  sensor,  and 
met  the  emergency  specifications  after  failure  of  any  two  sensors.  These 
robustness  properties  pertain  to  all  flight  conditions. 

Like  all  graphical  methods  this  design  technique  is  very  intuitive 
for  two  dimensional  problems,  but  not  as  well  suited  for  higher  dimensional 
problems,  where  the  dimension  here  refers  to  the  number  of  controller  para¬ 
meters  being  considered  in  a  particular  design  step.  To  make  the  problem 
tractable,  the  designer  has  to  break  the  problem  into  a  series  of  two  dimensional 
problems  by  fixing  the  additional  parameters.  While  in  this  study  the  proper 
sequence  was  clear,  for  higher  dimensional  problems  it  will  in  general  be  more 
difficult.  The  multi-input  and  general  dynamic  feedback  problems  are  particu¬ 
larly  difficult  because  of  the  increased  dimensionality.  As  seen  from  the 
results  of  this  Chapter,  however,  h -space  techniques  are  useful  tools  in  the 
design  of  control  systems. 
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CHAPTER  6 

AN  OPTIMIZATION  TECHNIQUE  FOR  ROBUST  CONTROL  SYSTEM  DESIGN 

6.1.  Introduction  ] 

This  chapter  formulates  an  alternative  design  technique  for  the 
control  of  systems  which  are  subject  to  large,  structured  perturbations.  The 
problem  formulation  is  the  same  as  that  of  Chapter  4.  It  is  assumed  that  a 
finite  number  of  operating  points  can  be  used  to  accurately  characterize  the 
perturbations.  Design  objectives  and  constraints  are  modeled  as  constraints 
on  the  locations  of  the  closed  loop  eigenvalues  in  the  complex  plane.  The 
objective  is  to  choose  a  fixed  gain  control  system  which  satisfies  the 
constraints  on  the  pole  locations. 

The  approach  of  this  chapter  is  to  use  a  quadratic  cost  functional  , 

to  represent  trade-offs  between  possible  design  points.  The  problem  is  refor¬ 
mulated  as  an  optimization  problem  over  the  free  parameters  of  the  control 
system.  This  optimization  is  then  solved  using  an  augmented  Lagrangian 

approach.  ' 

' 

The  outline  of  this  chapter  is  as  follows.  In  Section  6.2,  a  precise 
mathematical  formulation  of  the  problem  is  given.  Section  6.3  discusses  the 
nonlinear  programming  method  which  was  used  to  solve  the  problem.  Section  6.4 
derives  the  gradients  which  are  necessary  to  solve  the  problem.  In  Section  ' 

6.5,  a  second  order  numerical  example  is  presented.  The  purpose  of  presenting  < 

this  example  is  to  discuss  some  of  the  problems  involved  in  implementing  this 
method.  Finally,  this  approach  to  controller  design  is  applied  to  the  problem 
considered  in  Chapter  5  of  designing  a  fixed  gain  controller  for  the  linearized 
longitudinal  flight  dynamics  of  an  F4  aircraft.  The  resulting  design  is  compared 
to  the  corresponding  design  of  Chapter  5.  ! 

--  -  -  - -  -  ■  ■  -  -  -  -  -  -  i  -  ■  ■  . . -  - ^ 
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6.2.  The  Problem  and  Its  Reformulation 

The  purpose  of  this  section  is  to  present  a  precise  mathematical 
formulation  of  the  problem  discussed  in  the  Introduction  of  this  thesis  and  - 

UN 

then  to  reformulate  the  problem  in  a  form  that  is  computationally  easier  to 

work  with.  As  mentioned  in  the  previous  section,  this  chapter  deals  with  the 

problem  of  using  output  feedback  to  control  a  fixed  structure  system.  The 

design  problem  is  to  choose  the  constant  feedback  gains  which  are  best  with 

respect  to  some  cost  function  such  that  the  closed  loop  system  satisfies 

certain  design  specifications.  It  is  assumed  that  these  design  specifications  > 

can  be  represented  most  naturally  as  regions  in  the  complex  plane  where  the 

eigenvalues  of  the  closed  loop  system  must  be  located.  The  system  is  also 

assumed  to  be  linear  time  invariant. 

The  precise  problem  formulation  is  as  follows: 

min  J  =  ‘/2  E{ /  [xT(t)Qx(t) +  uT(t)Ru(t) ]dt} 
kES  o 

subject  to 

x(t)  *  Ax(t)  +  Bu(t);  x(tQ)  * 

E{x  }  =  0;  E{x  xT}  =  X 
o  0  0  o 

v_(t)  =  Cx(t) 

u(t)  *  -Kv(t) 

g.(\)  <0  i  -  1 , . . . , N 


4 


(129) 

(130) 

(131) 

(132) 

(133) 
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where  xCt)  £  R  ,  u(t)€R  ,  and  y(t)€R  .  S  is  the  space  of  permissible  feed¬ 
back  gains.  The  expectation  of  the  integral  is  used  so  that  the  cost,  J,  is 
independent  of  any  particular  initial  state  of  the  system,  but  depends  instead 
on  an  average  initial  condition  of  all  the  possible  initial  states.  Q  and  R 
are  nonnegative  definite  constant  matrices  chosen  so  that  given  some 
predetermined  criteria,  by  minimizing  J,  one  is  improving  the  closed  loop 
behavior  of  the  system.  For  example,  if  Q  is  the  zero  matrix  and  R  is  the 
identity,  J  represents  the  total  energy  used  to  control  the  system.  By 
minimizing  J,  one  is  minimizing  the  total  energy  used.  The  functions,  g^(\), 
represent  constraints  on  the  location  of  the  eigenvalues,  X,  of  the  closed 
loop  system  in  the  complex  plane. 

The  problem  (129)-(133)  looks  computationally  difficult  to  solve; 
however,  it  can  be  reformulated  as  follows: 


subject  to 


where 


min  J  =  ‘/jtrtMQOP} 

kes 

gi(A>  <  0  i  *  1, . . .  ,N 
S  *  {k/(A-BKC)  is  asymptotically  stable} 


(134) 


(135) 


(136) 


M(k)  -  Q  +  CTKTRKC 


(137) 


A(k)  =  A  -  BKC 


(138) 


X  =  E{x(0)xT(Q) } 

o  —  — 


(139) 


A(k) P  +  PA'  (k)  =  -X  . 
—  —  o 


(140) 


Here  ?  is  a  constant  positive  definite  matrix;  k  is  a  vector  comprised  of  the 
elements  of  the  matrix  K;  x(0)  is  the  initial  state  of  the  svstem;  and  A(k)  is 


the  closed  loop  system  for  a  given  k.  Given  this  formulation,  the  cost  J  is 
easily  found.  To  solve  for  J  directly  from  (129),  one  would  first  have  to 
calculate  3C ( t )  for  each  k;  whereas  to  solve  (134),  one  has  to  solve  the  Lyapunov 
equation  (140)  for  each  k  and  then  perform  a  few  simple  matrix  operations. 

The  derivation  of  (134)-(140)  is  a  direct  consequence  of  the  results  in  [60]. 

Problem  (134)-(140)  is  a  mathematical  representation  of  the  problem 
(described  in  the  introduction  of  this  chapter)  of  choosing  constant  output 
feedback  gains  for  a  linear  time  invariant  system  subject  to  certain  design 
criteria  which  are  represented  by  a  cost  function  and  constraints  on  the 
locations  of  eigenvalues  in  the  complex  plane.  The  question  of  how  to  solve 
the  problem  (134)-(140)  remains.  Problem  (134)-(140)  is  a  nonlinear  constrained 
minimization  over  a  finite  dimensional  space.  The  next  section  will  discuss 
methods  of  solving  such  a  problem. 

6.3.  Nonlinear  Programming  Solution  Procedure 

There  are  several  ways  to  solve  a  nonlinear  problem  of  the  form 
(134)-(140).  As  stated  previously,  the  basic  problem  is 

min  J(k)  (141) 

k€S 

subject  to 

3i(k)  <0  i =  1 , . . . ,N  <142> 

where  the  exact  form  of  the  cost  function  is  defined  explicitly  by  equations 
(134),  (139),  and  (140)  in  Section  6.2.  ?w°  of  the  most  common 

me  chads  of  solving  a  problem  of  this  type  are  penalty  function  and  Lagrange 
multiplier  mechods  [61 j.  Each  of  these  methods  has  certain  problems  (to  be 

inherent  to  the  method. 


described  below)  which  are 


However,  by  using  a 


combination  of  both  methods,  these  problems  can  be  avoided  and  a  better 
approximation  to  the  solution  for  (141)-(142)  can  be  obtained  [62]. 


In  its  simplest  form  the  Lagrange  multiplier  method  solves  the 
following  problem  [61] 


max  ^(d.) 


(143) 


subject  to 


d^  >  0  i*l,...,N 


(144) 


N 

where  d€  R‘  and 


N 

'P(d)  =  min  J(k)  +  E  d.g.(k).  (145) 

kES  i=1  1  1 

This  problem  is  often  easier  to  solve  since  the  nonlinear  constraints,  g^k), 
have  been  replaced  by  simple  linear  ones.  The  problem  (143) -(144)  is  the  dual 
problem  of  problem  (141)-(142) .  The  duality  theorem  states  that  as  long  as 
¥>(d)  >  for  some  positive  d^s  and  J(k)  <  00  for  some  k£S,  the  solution  to 
(143)-(144)  is  less  than  or  equal  to  the  solution  to  (141)-(142).  When  the 
solution  to  (  143)-(144)  is  strictly  less  than  the  solution  to  (141)-(142),  a 
duality  gap  exists  [61].  For  convex  functions  with  convex  constraints,  this 
difficulty  does  not  occur.  The  solution  to  (143)-(144)  is  also  the  solution 
for  (141)-(142).  However,  for  a  general  function,  J(k),  a  duality  gap  may 
exist  so  that  the  solution  to  (  143)  — (144)  is  a  lower  bound  on  the  solution  to 
(141)-(142),  rather  than  the  minimum  [61]. 

On  the  other  nand,  exterior  penalty  functions  solve  the  problem 


[151: 

min  J (k)  +  cH(k)  (146) 

k€si 

where  c  is  some  positive  constant,  S  is  the  region  in  R™  where  all  the 


constraints,  (142),  are  satisfied,  and  H(k)  is  a  functional  with  these 
properties : 


H(k)  >0  for  all  k€  Rm 

(147) 

H(k)  is  continuous 

(148) 

H(k)  =0  <=>  k£  . 

(149) 

As  long  as  kSS^  H(k)=0,  so  problems  (141)-(142)  and  (146)  are  identical. 
When  k  is  outside  S^,  the  function  J(k)+cH(k)  is  large.  As  c  becomes  large, 
the  minimum  of  J(k)  +  cH(k)  approaches  S^.  The  most  common  penalty  function 
is 

No 

H(k)  =  £  max[0,g. (k) ]  .  (150) 

i=l  1 

cor  this  function,  the  value  of  H(k)  is  the  sum  of  the  squares  of  the 
distances  by  which  each  constraint  is  violated,  so  the  penalty  term  increases 
rapidly  when  the  distance  k  is  outside  S^.  The  advantage  of  this  method  is 
that  problem  (146)  is  an  unconstrained  minimization  problem  which  is  often 
easier  to  solve  than  problem  (141)-(142).  The  disadvantage  of  this  approach 
is  that  to  obtain  a  good  approximation  to  problem  (141)-(142),  c  must  become 
large.  However,  as  c  approaches  infinity,  the  matrix  of  second  partial 
derivatives  of  J(k)+cH(k)  (the  Hessian)  becomes  increasingly  ill 
conditioned.  Many  algorithms  for  unconstrained  minimization  depend  on 
either  the  Hessian  or  an  approximation  of  the  Hessian  to  find  the  minimum. 

If  the  Hessian  is  ill-conditioned ,  these  algorithms  will  converge  very  slowly 

[61]. 

By  combining  penalty  function  methods  with  Lagrange  multiplier 
methods,  cne  can  eliminate  the  duality  gap  and  use  a  smaller  value  of  c,  thus 


142 


improving  the  conditioning  of  the  Hessian  at  the  solution  [62].  Using  the 
penalty  function  H(k)  ,  in  (150),  consider  the  problem: 


subject  to 


min  J(k)  +  cH(k) 
kSS 


gi(k)  <0  i =  1, . . . ,N. 


(151) 


(152) 


Given  the  properties  of  H(k)  ,  problem  (151)-(152)  is  equivalent  to  (141)-(142) 
The  dual  of  problem  (151)-(152)  is 


subject  to 


max  gc(d_) 


d.  >  0  i = 1, . . . ,N 

l  ’  * 


(153) 


(154) 


g  (d)  =  inf{J(k)  +  cH(k)  -f-  E  d.g.(k)}, 
c  1=1  1  1 


(155) 


Theorem  1:  There  exists  a  cs  R  with  0  <  c.  <  =°  such  that  the  solution  to  problem 

(153) -(154)  is  also  the  solution  to  problem  (151)-(152). 

Proof:  See  Bertsekas  [63]. 

Theorem  1  implies  there  is  no  duality  gap  for  problems  (151)-(152)  and  (153)- 

(154) .  Moreover,  since  the  value  of  c  needed  to  solve  this  problem  exactly 
is  finite,  the  structure  of  the  Hessian  is  more  favorable  for  solving  the 


oroblem. 


Bertsekas  [64]  discusses  a  variation  of  problems  (151)-(152)  and 


(153)  —  (159)  and  suggests  a  very  straightforward  way  to  solve  the  maxinizatii 
over  d_.  He  suggests  solving  the  problem 


max  g  (d) 
c  — 


(156) 


143 


subject  to 


d.  >  0  i  =  1,. .  .  ,N 

l 


(157) 


where 


g  (d,k)  =  infUQO+rV  Z  {maxfO.d. + c ' g. (k) ] }2-d2] } .  (158) 

kes  1=1  11  1 


For  the  case  where  the  ith  constraint  is  violated,  the  corresponding  term  in 
the  summation  is 

2 


d.g.(k)  +  c'/2  gj(k)  (159) 

which  is  identical  to  the  corresponding  term  in  (155)  for  c  =  c’/2  and  H(k)  as 
defined  in  (150).  For  the  case  where  the  ith  constraint  is  satisfied,  but 

d.  +  c ' g. (k)  >  0  (160) 

l  &i  — 

equation  (159)  also  applies;  and,  when 


di  +  c '  g„.  (k)  <  0 

he  corresponding  term  in  the  summation  is 


(161) 


1_  H2 
T  d- 


(162) 


2c’  i 

3e rtsekas  has  shown  that  the  solution  to  (156)-(157)  is  equivalent  to  that  of 
problem  (141)— (147)  for  all  values  of  c'  greater  than  some  lower  bound  c  (c 
exists  and  is  finite). 

One  can  solve  the  problem  (136)-(157)  iteratively,  viewing  the 
iteration  over  _d  as  a  fixed  stepsize  gradient  problem  [64].  The  partial 
of  g  (d,k)  with  resoect  to  d.  is 

C  ~~  1 

jg„  (d,jc) 


=  max  [-d .  /c  '  ,  g.(k)]  i=l,...,N. 

5a,  1  1  — 


(163) 


Hence  the  gradient  of  g  (d,k)  with  respect  to  _d  is  the  vector  of  these 

c  ■ 

partials.  The  appropriate  update  of  dj_1  is 


1 


iJ+1  =  dJ+c’Vgc(d,k). 


(164) 


For  J(k)  and  g^(k)  convex,  Bertsekas  has  shown  that  his  method  has  demon¬ 
strated  global  convergence  for  a  wide  range  of  step  sizes.  The  main  advan¬ 
tages  of  using  this  method  is  that  it  combines  the  advantages  of  both  penalty 
function  and  Lagrange  multiplier  methods  and  that  the  iterative  method  to 
solve  the  maximization  over  ji  is  very  simple. 

There  remains  the  problem  of  solving  the  minimization  over  k  for  a 
fixed  d.  This  problem  can  be  solved  using  a  variable  metric  algorithm  [65]. 
each  iteration  of  the  routine  the  user  must  supply  the  value  of  the  function 
co  bj  minimized  and  its  gradient.  From  this  information,  the  routine  builds 
up  an  approximation  to  the  inverse  Hessian  which  improves  as  the  routine 
gathers  information  from  more  points. 


To  solve  for  gc(d,k)  in  (158)  at  each  iteration,  one  must  solve 

for  J(k)  and  g .  (k)  ,  i  =  l,...,N.  For  problem  (134)-(140)  from  Section  2,  J(k) 

can  be  solved  using  (134)  and  (140).  The  constraints  g^(Mk))  are  chosen  by 

the  designer  and  thus  are  also  readily  available.  The  gradients  are  also 

needed  at  each  iteration.  Taking  the  partials  of  g  (d,k)  in  (158)  with 

c - 

respect  to  k, 


3g  (d,k)  3J(k) 

c -  - — 


•J(k)  3g.(k) 

—  +  max[0,  d  +c/g  (k)]  - 

i  '  i 


(165) 


Thus  to  solve  problem  (155),  the  gradient  of  the  constraints  and  of  the  cost 
with  respect  to  k  must  be  provided.  Section  6.4  discusses  the  computations 
required  to  solve  for  7 .  J(k)  and  7  g(k). 
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6.4.  Cost  and  Constraint  Gradient  Calculations 

As  described  in  Section  3  our  approach  to  the  solution  of  problem  (134) 
(140)  is  to  solve  the  equivalent  problem  (156)-(158)  iteratively.  To  solve 
the  minimization  over  k,  for  a  fixed  d,  the  gradient  of  J(k)  with  respect  to 
k  and  the  gradient  of  g^(A(k))  with  respect  to  jc  must  be  computed.  Using 
linear  operator  theory,  one  can  derive  a  fairly  simple  expression  for  the 
gradient  of  J(k).  Using  eigenvalue  sensitivity  theory,  one  can  derive  an 
expression  for  the  gradients  of  the  constraints  with  respect  to  k  [66]. 

The  gradient  of  J  (k)  with  respect  to  K  is  as  follows: 

Theorem  2:  V„J(k)  -  (RKCP-3TZP) CT 

K. 

where  P  and  E  are  solutions  of 


ATE  +  ZA  +  M(k)  =  0 

AP  +  PAT  +  X  *  0 
o 

and  A  and  M(k)  are  defined  in  (138)  and  (137)  respectively. 


(166) 


Proof :  Follows  directly  from  [60]. 

~,^J (k)  is  found  by  rearranging  7^ J (k )  (k  is  a  vector  of  the  elements  of  the 


matrix  K) . 


To  compute  the  gradient  of  the  constraints  with  respect  to  k, 
one  must  compute  the  gradient  of  the  eigenvalues  with  respect  to  k.  After 
finding  the  latter,  one  can  use  the  rules  of  implicit  differentiation  to 
find  the  former.  Consider  equation  (167). 


w  *  ( A-  I )  v  =  0 


(167) 


vr.ere  v*  is  the  left  eigenvector  of  A  and  v  is  the  right  eigenvector 

T* 

corresponding  to  eigenvalue  '■  .  Since  w1  and  v  are  the  eigenvectors  ci  A, 
(167)  is  true.  Vnen  the  Frechet  differential  of  a  function  F(k)  exists,  it 


<5  F(k,Ak) 


^F(k+£Ak)|£=0. 


(168) 


The  Frechet  differential  (168)  can  also  be  written  in  terms  of  an  inner 
product  as 


where 


5  F(k,Ak)  a  <7  F(k)  ,Ak> 


<7F(k),Ak>  -  tr{7F  (k)Akh 


(169) 


(170) 


Take  the  Frechet  differencial  of  both  sides  of  (167) 


dw^(A-AI)v  +  w^(A-AI)<Sv  +  w^(<5A-<5AI) v  *  0. 


(171) 


Again  since  w1  and  v  are  eigenvectors,  the  first  two  terms  are  identically 
zero,  thus 


T  T  - 

w  (JAv  *  w  5Av. 


Since  oA  is  a  scalar. 


T  T  ' 

6  Aw  v  =  w  <5Av. 


(172) 


(173) 


Since  w  v  is  also  scalar 

w^SAv 


(174) 


T  T 

Notice,  however,  that  if  w  v=  0,  equation  (174)  will  not  hold.  If  w  \;  =  0, 

T 

then  the  left  eigenvector,  w  ,  is  perpendicular  to  the  right  eigenvector,  v. 
This  only  happens  when  A  has  a  Jordan  block  of  dimension  greater  than  one. 


Continuing,  from  definition  (138)  and  (169) 


5A(K,AK)  =  A(K+  eAK)  j  , 


=  —  [A-B(K+  SAK)C]I  _n 

ct  t—'J 


-3AKC. 


(175) 


Using  definition  (169) 


<5A(K,AK)  -<  VX(k),  AK> 


tr{VAi(K)AK} . 


(176) 


Substituting  (175)  and  (176)  in  (174): 


T  vT(-BAKC)v 

tr{VA  (k)AK>  = 


T 

w  v 


-  trt 


/_ 


wT(BAK)cv 


T 

W  V 


-} 


T 

Cv  w  B 

-  tr{ AK  }  , 
w  v 


(177) 


Since  (177)  must  hold  for  an  arbitrary  AK, 


V  A  (K )  =  - 


CvwTB 


L  W  V  j 


(178) 


or 


VA(k) 


BTwvV 


T 

w  v 


(179) 


”,  A(k)  is  found  by  rearranging  7  A(K)  (k  is  a  vector  comprised  of  the  element 


of  K)  . 


Define  A.  in  terms  of  two  real  variables,  a.  and  a.. 

l  li 


A .  =  o .  +  j  w  . 
l  i  i 


(180) 


Then 


?o_.  5  A . 

liT  *  R*al(lki) 

1  1 


(131) 


u ,  3  A  . 

~  =  Imaginary  (— )  . 


3k . 

l 


(182) 


The  constraint  functions  g^.(A)  from  equation  (135)  are  considered  functions 
cf  the  two  real  variables,  j  and  u.  For  the  purposes  of  problem  (1 34)-(liO) 
each  constraint  will  be  a  function  of  onlv  one  eigenvalue.  If  all  the 
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eigenvalues  must  lie  inside  a  particular  boundary,  then  n  of  the  constraints 
will  be  the  equation  for  the  boundary  (one  for  each  eigenvalue).  Given  this 
situation. 


3k 


3gi(<?i  .UK) 
3o, 


^  3k  )  +  l  3uj 


(183) 


Since  the  regions  in  the  complex  plane  are  chosen  by  the  designer,  it  will 
be  assumed  that  the  regions  are  chosen  so  that  the  partials  with  respect  to 
o  and  u>  exist.  For  the  same  reason,  the  functions  g  (a,t»j)  are  known  explicitly, 
and  thus,  so  are  the  partials.  Thus  using  equations  (179)-(183),  the  gradient  l 
of  the  constraints  with  respect  to  the  feedback  gains,  k,  can  be  calculated. 

In  summary,  both  the  gradients  of  the  constraints  with  respect  to 
k  and  the  gradient  of  the  cost  J(k),  as  well  as  the  values  of  the  constraints 
and  the  cost,  can  be  calculated  given  a  point  k.  Using  this  information,  one 
can  find  the  solution  to  (134)-(140)  by  solving  the  equivalent  problem  (15 3)  —  (155 ) 
as  described  in  Section  6.3.  The  next  section  will  discuss  some  of  the  specific 
details  and  problems  involved  in  implementing  this  method  to  solve  the 
reformulated  problem  (134)-(140). 


6.5.  Second  Order  Example 

The  purpose  of  the  numerical  example  of  this  section  is  to  demon- 
strata  how  well  the  method  developed  in  Section  6.2  to  solve  problem  (129)- (134) 
works  on  a  simple  second  order  example.  The  problem  is  as  follows: 

min  J  *  */2 /  u  (t)u(t)dt  (134) 

k£S  O 


subject  to 


x(t)  =  Ax(t)  +  Bjj(t) 


(185) 


u(t)  =  -[^  k2]x(t) 

S  =■  (k/(A-Bk)  is  asymptotically  stable} 


■-1 

5  ' 

'o' 

B  * 

_-5 

-1_ 

1 

gi(a,tj)  <0  i  =  1,4 


(186) 

(187) 

(188) 

(189) 


where  a  and  w  are  the  real  and  imaginary  parts  of  the  eigenvalues  of  the 
closed  loop  system.  The  constraints  are  (see  Figure  6.1): 


g1(o.w) 

=  a)  -  2 . 6a 

(190) 

g2(o»w) 

-  w  +  2 . 6  a 

(191) 

g3(a,m) 

=  4.0804-  a2  -w2 

(192) 

g4(a,m) 

*  c2  +  u2  -  53.1441. 

(193) 

Each  of  these  four  equations  must  be  satisfied  for  both  eigenvalues  so  there 
are  actually  eight  constraints. 

Examining  a  second  order  system  with  a  single  input  is  particular! 
convenient  for  demonstrating  the  behavior  of  this  algorithm.  First,  for  a 
second  order  system  one  can  derive  explicit  equations  relating  the  feedback 
gains  to  the  eigenvalues  of  the  closed  loop  system.  Second,  equations 
mapping  the  boundaries  in  the  complex  plane  to  boundaries  in  k-space  (the 
space  of  feedback  gains)  can  be  obtained  using  the  mapping  method  described 
Chapter  4.  For  this  second  order  system  and  reasonable  boundaries,  like  thos 
given  in  (190)- (193),  the  boundaries  in  k-space  are  not  too  complex  (Figure  6. 
Since  the  minimization  is  actually  over  k  in  the  k-plane.  Figure  6.2  shows 
exactly  what  the  constraints  are  in  this  space. 


Figure  6.1.  Second  order  example:  constraint  region  in  the  complex  plane 


The  closed  loop  characteristic  equation  for  this  system  is 


X  +  (2  +  k2)X  +  (26  +  5^  +  ^) 


0. 


(194) 


Applying  the  quadratic  formula  to  (194)  yields 

2  l,t 

X  ,X  *  -<1+-5k2)  ±  -5(k2  -  20kx  -  100)  . 

Taking  the  partials  of  X^  and  X2with  respect  to  k^  and  k2» 


3X 


-  +  10  (k2  -  20kx  -  100)  "/2 


(195) 


(196) 


3X 


3k 


^  «  -.5  +  .5k2(k2-  20^-100) 


-1/, 


(197) 


From  (196)  and  (197),  one  can  see  that  these  partials  have  discontinuities 
precisely  at  the  boundary  where  the  closed  loop  system  poles  change  from  a 
complex  pair  to  two  real  poles  or  vice  versa.  Not  only  are  these  partials 
discontinuous  at  this  boundary,  their  magnitude  approaches  infinity  as  the 
poles  approach  this  boundary.  The  equation  of  this  boundary  in  k-space  is 

2 


k“  -  20kx  -  100  =*  0 


(198) 


which  corresponds  to  boundary  9  in  Figure  6.2. 

The  other  boundaries  in  Figure  6.2  correspond  to  the  boundaries  in 
the  complex  plane  (Figure  6.1) as  follows:  the  large  circle  (Figure  6.1)  maps 
into  the  triangle  formed  by  1,  2,  and  3  (Figure  6.2),  the  small  circle  maps 
into  the  triangle  formed  by  4,  5,  and  6, 'the  two  lines  into  boundaries  7  and 
8.  The  region  enclosed  by  the  solid  line  in  Figure  6.2  is  the  region  in  the 
k-plane  where  all  the  cons traints (190) - (19 3)  are  satisfied.  The  reason  for 
choosing  these  boundaries  in  the  complex  plane  is  that  such  boundaries  do 
occur  in  real  problems  (e.g.,  the  aircraft  example  in  the  next  section)  as 
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constraints  on  the  locations  of  closed  loop  system  poles.  A  nice  feature  of 
using  circles  for  boundaries  in  the  complex  plane  is  that  for  second  order 
systems  circles  map  into  triangles  in  the  k-plane.  Hence  some  of  the 
boundaries  in  the  k-plane  are  straight  lines  (Chapter  4). 

The  problem  (184)- (189) is  to  find  the  minimum  energy  control  subject 
to  the  indicated  constraints.  Since  system  (185)  is  stable,  the  minimum 
energy  control  without  constraints  corresponds  to  zero  gain.  With  the  feed¬ 
back  gains  set  to  zero,  the  poles  of  (185)  are 

VX2  *  ~L- j5‘  (199) 

With  the  given  constraints (Figure  6 . 1)  ,  the  minimum  energy  feedback  gains  were 
found  to  be 

k  =  (-.273  +1.68]  (200) 

which  places  the  closed  loop  poles  at 

,\  ,X2  =  -1.34  +  j4 . 788 .  (201) 

This  answer  makes  sense.  From  Figure  6.2,  one  can  see  that  this 
point  k  (points  b)  is  approximately  the  point  inside  the  constraint  region 
closest  to  the  origin.  From  Figure  6.1,  one  can  see  that  the  closed  loop 
eigenvalues  (points  a)  are  about  as  close  to  the  open  loop  eigenvalues 
(points  b)  as  possible  given  the  constraints.  The  algorithm  converged  to 
the  minimum  k  (200)  for  a  wide  range  of  initial  guesses  for  k.  Initial 
guesses  for  k  which  placed  the  closed  loop  poles  outside  the  large  circle 
and  to  the  left  of  boundary  1  were  the  only  ones  for  which  the  algorithm  did 
not  converge  to  the  value  of  k  given  in  (200). 

The  region  in  the  complex  plane  for  which  the  algorithm  did  not 
converge  corresponds  to  the  area  in  the  k-plane  (Figure  6.2)  to  t:,e  richt  of 


line  1  and  just  below  curve  9.  The  fact  that  the  algorithm  could  not  converge 
from  these  points  can  be  explained  by  the  discontinuities  in  the  partials  of 
the  eigenvalues  with  respect  to  k  mentioned  in  (196)- (198).  Consider  point  c 
(Figure  6.2)  as  a  typical  point  in  this  region.  It  corresponds  to  a  complex 
pole  pair  outside  the  large  circle  in  Figure  6.1.  Line  1  in  the  k-plane 
represents  the  boundary  for  a  complex  pole  pair  crossing  this  circle;  thus, 
the  negative  of  the  gradient  in  the  k-plane  for  these  points  points  towards 
line  1  and  nearly  perpendicular  to  it.  Moving  in  this  direction  should 
reduce  the  cost  function  (158) .  However  notice  that  from  point  c,  for 
example,  movement  in  this  direction  will  lead  to  guesses  for  k  which  fall 
above  or  on  curve  9.  A  point  on  curve  9  corresponds  to  a  double  real  root 
for  which  the  partials  of  the  eigenvalues  with  respect  to  k  are  infinite. 

This  will  obviously  cause  problems.  Notice  also  that  for  two  real  poles, 
the  direction  ot  decreasing  cost  is  determined  by  boundary  2;  whereas,  for  a 
complex  pair,  the  direction  of  decreasing  cost  is  determined  by  line  1.  For  all 
these  points  -points  for  which  the  routine  would  not  converge  to  (200)  -  the 
minimization  routine  found  points  which  approached  boundary  9.  However,  the 
algorithm  was  not  able  to  move  across  or  along  the  boundary .  In  summary', 
the  derivatives  of  the  constraints  with  respect  to  the  feedback  gains  are 
not  continuous.  This  fact  can  lead  to  convergence  problems.  However,  one 
can  avoid  these  problems  by  choosing  a  better  initial  guess  for  k  (and 
lower  values  for  c) . 

This  example  was  also  used  to  study  the  behavior  of  the  algorithm 
with  respect  to  changes  in  the  constant  c'  in  (158).  As  discussed  in  the 
section  on  nonlinear  programming,  for  a  very  large  value  of  c' ,  the  minimi- 
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zation  over  k.  converges  Co  a  solution  which  is  close  to  the  solution  to  the 

actual  problem,  in  this  case  (184)-(189).  For  a  smaller  c',  each  iteration 

over  k  stops  farther  from  the  real  solution  than  with  a  larger  c',  but  the 

iterations  over  d_  lead  more  quickly  to  the  true  solution  of  the  problem.  For 

this  problem,  c'  equal  to  2000  seemed  to  work  best.  The  minimization  over  k 

led  to  a  solution  which  was  very  close  to  the  final  solution  of  the  algorithm. 

The  iterations  over  _d  merely  served  to  bring  the  point  a  bit  closer  to  the 

-7  -3 

boundaries  (within  10  ,  instead  of  10  ) .  This  was  true  even  for  smaller 

values  of  c' ,  20  and  200.  The  smaller  values  of  c?  led  to  more  iterations 
over  _d,  but  fewer  over  k  at  each  substep.  For  c’  equal  to  2,  Che  first 
iteration  over  k  did  converge  to  a  solution  which  was  different  from  the 
solution  with  c'  equal  to  2000.  However  the  iterations  over  _d  led  to  the 
same  final  solution  as  with  c’  equal  to  2000.  In  terms  of  total  function  ■ 
evaluations,  c'  equal  to  2000  was  the  most  efficient;  moreover,  the  solutions 
for  smaller  values  of  c'  were  not  significantly  different  from  those  with  c' 
equal  to  2000. 

To  summarize  the  results  for  this  example,  this  algorithm  works 
provided  a  good  initial  guess  for  k  and  a  reasonable  value  of  c'  are  used. 
Provided  these  two  conditions  are  satisfied,  the  minimum  energy  feedback  gains 

for  problem  ( 184 )  —  ( 1S9 )  are 

k  =  [-.273  ml. 68].  (2C2) 

This  gain  places  the  closed  loop  system  poles  at 

•  =  -1.84  +  j  4 . 738 .  (20  3) 

The  next  two  sections  will  present  a  more  complex  example  and  discuss  some 
results  for  that  examrle. 
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6.6.  F-4E  Example 

The  problem  studied  in  this  example  is  the  same  as  that  of  Chapter  5. 
The  system  is  a  linearized  model  of  the  longitudinal  motion  of  the  F-4E 
Phantom.  The  model  equations  and  design  criteria  are  summarized  here  for 
clarity . 

The  system  equations  are: 


_d_ 

dt 


N 

z 

all 

a12 

a!3 

N 

z 

“ 

bl 

q 
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a21 

CN 

<N 

a23 

q 

+ 

0 

6 
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0 

0 

-a 

6 

e 

a 

(204) 


10  0 

Nz 

y(t)  = 

0  10 

q 

0  0  1 

6 

e 

(205) 


u  ( t)  =  - ky  ( t ) . 


(206) 


The  actuator  bandwidth,  a  was  assumed  to  be  14  rad/s.  The  values  a  and  b 

ij  1 

are  different  for  each  flight  condition  and  are  given  in  Appendix  I.  The  matrix 
k  is  to  be  determined. 

One  design  problem  for  this  airplane  was  to  choose  k  such  that  the 
closed  loop  eigenvalues  for  each  flight  condition  are  in  certain  regions  in 
the  complex  plane  (Figure  6.3).  Ideally  one  would  like  to  find  one  set  of 
gains  which  worked  for  all  four  flight  conditions.  The  constraints  on  the 
short  period  eigenvalues  are  given  by  restrictions  on  the  damping  and  the 
natural  frequency  of  the  short  period  mode.  The  characteristic  equation  for 


these  eigenvalues  is 


X  +  2£  w  X  +  a) 

sp  sp  s p 


(207) 


where  5^  is  the  damping  and  ui^  is  the  natural  frequency.  Under  normal 
operating  conditions  (i.e.  ,  no  sensor  failures),  £gp  and  to  are  required 
to  satisfy 


.35  <  C  <  1.3  (208) 

sp 

“a  <  usp  <  “b  (209) 

where  and  depend  on  the  flight  condition  (see  Table  6.1).  For  the  case 

when  one  of  the  sensors  fail,  C  and  u>  must  satisfy 

’  sp  sp  } 

.15  <  5  (210) 

sp 

cj  <  oj  (211) 

c  sp 

where  u>  depends  on  the  flight  condition  (Table  6.1).  For  the  emergency 
situation,  the  actuator  pole  is  also  required  to  satisfy  (210)-(211).  For  the 
nonemergency  situation,  the  actuator  pole  (A  )  is  required  to  satisfy 

cl 

0),  <  A  <  70  rad/s.  (212) 

b  a 

These  regions  in  the  complex  plane  are  shown  in  Figures  6.3  and  6.4. 

For  this  chapter,  an  additional  design  criterion  was  added.  The 
feedback  gains  were  to  be  chosen  such  that  the  minimum  total  control  energy 
is  used  given  the  constraints  on  the  locations  of  the  closed-loop  poles 
described  above.  The  appropriate  cost  function  to  minimize  is 

50  T 

min  J  *  /  u  (t)u(t)dt  (213) 

k£S  o 


where  S  is  the  set  of  feedback  gains  for  which  the  closed  loop  system  (204)- 
(206)  is  asymptotically  stable.  Taken  together  with  appropriate  equations 
for  the  constraints  in  Figures  6.3  and  64,  equations  (204)— (206)  and  (213)  repre 
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sent  a  problem  of  the  form  (129 ) — (134 )  in  Section  6.2.  Thus  the  method  of  soluti 
developed  in  the  previous  sections  of  this  report  can  be  applied.  The 
results  for  this  example  are  discussed  in  the  next  section. 


Table  6.1.  Frequency  limits  which  determine  boundaries 
in  Figures  4  and  5 


u  (rad/s) 
a 

ub (rad/s) 

oj  (rad/s) 
c 


Flight 

Condition  #1 

Flight 

Condition  #2 

Flight 

Condition  # 3 

Flight 

Condition  #4 

2.02 

3.50 

2.19 

3.29 

7.23 

12.6 

7.86 

11.8 

1.53 

2.65 

1.65 

2.49 

6.7.  Results  and  Discussion  for  Airplane  Example 


The  design  problem  for  the  F4-E  airplane  is  to  find 
one  set  of  constant  feedback  gains  for  which  the  closed 
loop  system  poles  are  in  the  appropriate  region  in  the  complex  plane  for  each 
one  of  the  four  flight  conditions  under  normal  operating  conditions.  After 
finding  such  a  solution,  the  next  problem  is  to  look  for  a  set  of  gains  that 
satisfies  the  above  criteria  and  also  is  robust  with  respect  to  sensor 
failures . 

First,  each  flight  condition  was  studied  separately  to  see  if  a 
fixed  gain  controller  could  be  found  to  satisfy  the  constraints  under  normal 
operating  conditions.  As  mentioned  in  Section  6.5,  a  good  initial  guess  for 


the  feedback  gains  is  important  for  the  algorithm  to  converge  properly. 

Since  the  problem  is  to  find  the  constrained  minimum  for  J  in  (213), 
one  sensible  starting  point  is  the  set  of  gains  which  correspond  to  the 
unconstrained  minimum  of  J.  One  can  find  these  gains  simply  enough  by 
solving  a  Riccati  equation  for  each  set  of  system  matrices  [68].  These  gains 
were  used,  and  gains  for  a  fixed  gain  controller  were  found  for  each  flight 
condition.  The  next  step  was  to  find  one  set  of  gains  which  would  work  for 
all  four  flight  conditions. 

Such  a  set  of  gains  was  found  both  for  the  case  when  all  three 
states  (204)  were  available  and  for  the  case  when  only  the  first  two  were 
available.  Chapters  considered  the  latter  case.  Rather  than  looking 

for  a  particular  set  of  gains,  the  procedure  used  was  to  map  the 
constraints  from  the  complex  plane  (Figure  6.3)  into  the  space  of  feedback 
gains.  The  entire  region  of  possible  gains  which  satisfy  the  constraints 
(Figure  6.5)  was  found.  Using  the  method  described  in  this  thesis,  the 
minimum  energy  k  using  only  two  gains  was  found  to  be  (see  Table  6.2) 


k  =  [-2.8281124 x 10”2  -2.0652172  * 10_1] . 


(214) 


This  point  is  marked  in  Figure  6.5  and  is  near  the  boundary  of  the  enclosed 
region  found  by  Franklin,  at  the  point  approximately  nearest  the  origin. 
Thus  the  results  presented  in  this  thesis  are  consistent  with  Franklin's. 


For  the  case  with  three  feedback  gains,  the  minimum  energy  Jc  was  found  to  be 
(see  Table  6.5) 

k  =  (-3.3498269  x  IQ-2  -2  .  75  74095  <  IC-1-  3 . 3295187  *  IQ-1]  (215) 


An  important  measure  of  svstem  oerforman.ee  is  the  C*  resoor.se  discussed 

h 


in  Charter  5. 


The  C*  is  a  line>r  combination  of  the  normal  acceleration 


Table  6.2.  Minimum  energy  feedback  gains:  normal  operating 
condition  solution  (2  gains) 


k-  [-2.8281124xl0-2  -2.0652172X10”1] 

Flight 
Condition  It 

Closed  Loop  Eigenvalues 

1 

-2.0483019,  -2.0200008,  -14.537826 

2 

-3 . 1624007  +  j  5 . 2143045 ,  -18 . 493321 

3 

-1.589176  +  J1. 7992842,  -14.547498 

4 

-2.1992634  + j 5. 8851873,  -16.308839 

Table  6.3.  Minimum  energy  feedback  gains:  normal  operating 
condition  solution  (3  gains) 


k-  [-3. 8498269x10” 2  -2 . 7574095x10” 1  3.3295187x10  1] 


Flight 
Condition  It 

— 

Closed  Loop  Eigenvalues 

1 

-2.0194637+ j 3. 3183791xl0'9,  -20.227559 

2 

-3.3792113  +  j 5. 0261318,  -25.502133 

3 

-1.59  7  7817  +  j 1.7690526,  -20.060981 

4 

-2.2968881  +  j 5. 7580029,  -22.569046 
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and  the  pitch  rate  of  the  plane,  given  by 

C*  -  (N  +  X2.43q)/k*  (216) 

N  z  n  c 

where  k  is  the  stationary  value  of  C*  and  is  used  for  normalization.  The 
c  N 

C*  response  to  a  step  input  should  fall  in  the  region  shown  in  Figure  6.6. 

N 

Figure  6.7  shows  this  response  for  each  of  the  four  flight  conditions.  The 
first  column  consists  of  the  responses  for  the  design  presented  in  this 
section  with  k  as  given  in  equation  (215) .  Comparing  these  responses  with 
Figure  6.6,  one  can  see  that  they  do  lie  within  the  required  region.  The 
second  column  contains  the  C*  responses  for  the  following  gain  matrix 


k  «  [-.115  -.8] 


(217) 


which  is  the  design  used  in  Chapter  5,  These  C^  time  responses  appear  to  be 
faster  and  to  satisfy  the  requirement  given  by  Figure  6-6  better  than  those 
presented  in  this  section.  This  is  not  surprising;  the  gains  (217)  were 

chosen  on  the  basis  of  the  C*  criterion.  The  criterion  used  to  choose  the 
gains  for  this  thesis  was  the  minimization  of  the  control  energy.  Thus, 
slower  C*  responses  should  be  expected. 

More  specifically,  the  design  criterion  used  to  choose  the  gains  in 
equation  (215)  was  the  minimization  of  the  control  energy  required  to  bring  the 
system  back  to  equilibrium  from  a  disturbance.  Figure  6.8  shows  u(t)  for  each 
of  the  four  flight  conditions  (the  first  column  contains  the  ones  for  this 
section; the  second  for  those  of  Chapter  4).  From  these  figures,  one  can  see  that 
the  controls  for  this  section  are  considerably  smaller  than  the  controls  which 
result  from  Chapters.  However,  the  system  is  stabilized  faster  (but  at  the 
expense  of  actuator  control)  for  the  gain  (217). 


I 


i 


I 


0.5 


1.0 

Time  (seconds) 


Figure  6.6.  C*  response  envelop 


Chapter  6 


Chapter  5 


Lgure  6.8.  Control  incuts  for  tvo  designs. 


The  feedback  gains  of  Chapter  5  were  chosen  by  looking  at  the  time  ^ 

responses  for  several  points  and  picking  the  best  one.  For  this  section,  the 
gains  (215)  were  chosen  by  using  a  cost  function  which  represented  the 
minimum  energy  control.  By  changing  the  cost  function  or  the  constraint  • 

boundaries  in  the  complex  plane,  one  could  easily  incorporate  the  C*  response 
criterion  into  the  design.  One  could  also  choose  the  cost  and  the  boundaries 
so  that  the  solution  would  be  a  compromise  between  the  minimum  energy  control  • 
and  a  fast  response.  Moreover,  the  results  from  this  design  indicate  where 
trade-offs  can  be  made  and  how  to  make  them.  In  short,  using  a  cost  function 
to  choose  a  set  of  feedback  gains  may  provide  more  insight  than  trial  and  * 

error  alone. 

Chapter  5r  also  mapped  the  emergency  regions  in  the  complex 

plane  (Figure  6.4)  into  the  gain  space  for  the  case  with  only  two  feedback  * 

gains.  Unfortunately, the  region  for  which  all  the  constraints  for  both  the 
normal  and  the  emergency  situations  are  satisfied  does  not  intersect  either 
axis  in  the  gain  space.  For  the  problem  with  only  the  first  two  states 
available  for  feedback,  this  means  no  set  of  gains  satisfying  the'constraints 
is  robust  with  respect  to  the  failure  of  either  sensor. 

In  light  of  the  results  of  Chapter  5,  the  problem  of  a  robust  control! 

was  considered  for  the  case  with  three  feedback  g3ins .  For  flight  condition  4, 

a  fixed  gain  using  only  the  first  two  states  was  found  which  satisfied  all 
the  constraints  and  was  robust  with  respect  to  either  the  first,  the  second, 
or  both  sensors  failing.  Adding  the  third  gain  set  equal  to  zero,  yields  a 
set  of  three  gains  which  is  completely  robust.  However,  this  result  is  not 
surprising  since  the  open  loop  poles  for  flight  condition  4 
already  satisfy  the  emergency  specifications.  For 


169 

the  other  three  flight  conditions  considered  individually,  no  gains  could  be 
found  which  were  robust  with  respect  to  the  first  sensor  failing.  Thus  none 
could  be  found  for  all  the  flight  conditions  taken  together.  When  considered 
individually,  a  set  of  gains  which  is  robust  with  respect  to  the  failure  of 
the  second  sensor  was  found  for  each  flight  condition.  A  set  of  gains  which 
is  robust  with  respect  to  the  second  sensor  failing  was  also  found  when  the 
first  three  flight  conditions  were  considered  together  (see  Table  6,4). 
Unfortunately,  when  all  four  flight  conditions  were  considered,  no  common 
solution  which  was  robust  with  respect  to  the  second  sensor  could  be  found. 

A  solution  w'hich  is  robust  with  respect  to  the  third  sensor  failing  is  just 
the  solution  given  in  Table  6.2  with  a  third  gain  equal  to  zero  added.  These 
results  seem  to  indicate  that  a  fixed  gain  controller  is  not  adequate  to 
satisfy  the  robustness  requirements  for  this  example. 

While  studying  this  example,  some  of  the  problems  in  implementing 
the  algorithm  for  the  second  order  example  in  Section  6.5  were  also  problems 
for  this  example.  First,  the  gradient  of  the  cost  with  respect  to  the 
feedback  gains  is  discontinuous  at  a  double  real  pole  (see  Section  6.4).  For 
the  minimum  energy  gains  in  Table  6.2,  the  eigenvalues  of  the  first  flight 
condition  seem  to  be  converging  to  a  double  real  pole  on  the  boundary  of 
the  constraints.  Since  the  gradient  is  discontinuous  at  this  point,  it  was 
necessary  to  try  initial  guesses  for  k  close  to  the  apparent  solution  but  on 
both  sides  of  the  discontinuity  to  be  sure  the  algorithm  was  not  hanging  up 
there.  The  algorithm  converged  back  to  the  double  pole  from  both  directions 
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Table  6.4.  ftinimum  energy  feedback  gains  which  are  robust  with 
respect  to  the  failure  of  the  second  sensor  (for  the 
first  3  flight  conditions  only) 


NF:  Se 
F:  Se 


k  =  [-5 .01384 77 xlQ-2  -4.0115944X10'1  5.0676513x10 


'Sensor  #2 

Closed  Loop  Eigenvalues 

NF 

-2.7084034,  -2.0196952,  -23.109953 

F 

-.72136579  + jl. 3337152,  -76.395321 

1  1  1  r 

NF 

-4.0968046+ j5. 1326,  -29.668796 

F 

- . 77106637  +  j 5.0911232  ,  -36 . 320273 

NF 

-1.8964697+ j 1.7951444,  -22.887455 

F 

-.49872359+  j2. 0622586,  -25.682948 

2  has  not  failed 
2  has  failed 


double  root  boundaries  away  from  the  solution.  In  that  case,  a  pole  placement 
which  places  the  poles  on  the  opposite  side  of  the  double  root  boundary  could 
be  done.  Using  those  gains  as  a  new  starting  point  might  allow  the  algorithm 
to  converge  to  the  real  solution. 

Another  observation  was  that  the  value  chosen  for  c'  in  (158)  affected 
the  final  solution  returned  by  the  algorithm.  For  large  values  of  c' 

(200-2000),  the  algorithm  converged  quickly  to  the  boundary  of  the  constraint 
region,  but  had  trouble  moving  along  the  boundary  to  the  minimum  with 
iterations  over  ji  (156)-(157).  As  explained  in  Section  6.3,  this  is  due  to  the 
ill  conditioning  of  the  Hessian  for  large  values  of  c'.  For  smaller  values 
of  c',  the  minimization  over  k  converged  to  a  solution  outside  the  boundary 
of  the  constraint  region.  The  maximization  over  _d  forced  the  solution  to 
the  boundary.  For  this  particular  example,  choosing  c'  equal  to  a  snail 
number  for  the  initial  iteration  over  k_  and  then  increasing  it  gradually  for 
subsequent  iterations  to  enforce  the  constraints  more  quickly  seemed  to  work 
well.  The  results  in  this  section  are  for  c'  equal  to  1  for  the  initial 
iteration  and  doubled  thereafter. 

Summarizing  the  results  of  this  section,  the  algorithm  developed  in 
this  chapter  was  applied  to  the  problem  of  designing  a  controller  for  the 
Fh-E  aircraft.  A  fixed  gain  controller  was  found  which  satisfied  the 
design  specifications  under  normal  operating  conditions.  However,  a  fixed 
gain  controller  which  was  robust  to  either  ser.sor  failure  was  not  possible. 

This  example  also  served  to  re-emphasire  some  of  the  inherent  problems  with 
this  design  technique:  the  discontinuities  of  the  gradient,  the  initial 


guess  for  k,  and  the  choice  of  c' 


6.8.  Conclusion  and  Summary 

This  chapter  has  dealt  with  one  method  of  solving  the  problem  of 
designing  a  fixed  gain  controller  for  a  linear  time  invariant  system  when 
some  of  the  design  criteria  are  represented  as  constraints  on  the  location  of 
the  closed  loop  system  eigenvalues  in  the  complex  plane  and  others  are  repre¬ 
sented  by  a  quadratic  cost  function  which  is  to  be  minimized.  First,  the 
original  problem  (129)-(133)  is  reformulated  to  yield  (134)-(140).  In  the 
form,  the  cost  function  is  easier  to  calculate.  Second,  problem  (134 ) — (140) 
solved  by  an  augmented  Lagrangian  method.  The  problem  is  a  max-min  problem. 

A  variable  metric  method  is  used  to  solve  the  minimization  over  ik.  A  fixed 
step  size  method  is  used  to  solve  the  maximization  over  d  [64],  Third, 
expressions  for  the  gradients  needed  to  solve  (153) -(155)  are  derived. 

Two  examples  are  studied:  a  simple  second  order  numerical  example 
and  a  model  of  the  longitudinal  motion  of  a  F4-E  plane.  Both  examples  serve 
to  point  out  several  problems  with  implementing  the  solution  of  (153)-(155). 
First,  there  are  discontinuities  in  the  gradients  for  the  case  when  a  complex 
pair  of  sy:  tern  poles  change  to  a  real  pair  or  vice  versa.  If  the  algorithm 
gets  hung  up  at  such  a  point,  a  new  starting  point  on  the  other  side  of  the 
boundary  may  help.  Second,  the  value  of  c'  in  (155)  mist  be  chosen  appro¬ 
priately  for  a  given  problem.  Too  large  a  value  of  c'  causes  slow  convergence 
of  the  algorithm.  Too  small  a  value  yields  a  solution  outside  the  boundary. 
The  best  approach  seems  to  be  to  choose  c'  small  for  the  first  iteration  and 
increase  it  thereafter.  Finally,  a  good  initial  guess  for  k  is  important  in 


order  for  the  algorithm  to  converge  properly. 


3 


There  are  advantages  to  this  design  method.  First,  it  has  the 
ability  to  incorporate  diverse  design  criteria  such  as  minimum  energy 
control,  rate  of  change  of  input,  constraints  on  the  location  of  poles  in 
the  complex  plane,  etc.  Second,  multi-input,  multi-output  systems  can  be 
considered.  Third,  it  has  the  ability  to  handle  larger  systems  than  some 
of  the  other  methods  used  to  solve  this  type  of  problem.  Finally,  it  provides 
insight  into  the  effects  of  the  various  design  constraints.  It  can  be  used 
to  determine  which  of  the  design  specifications  can  be  satisfied  and  which 


ones  may  be  too  stringent. 
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CHAPTER  7 


SUMMARY  AND  FUTURE  RESEARCH 


As  Chapter  2  illustrates,  the  robust  control  problem  statement 
encompasses  a  great  variety  of  problems.  However,  the  work  can  be  classified 
into  two  general  areas.  The  first  models  the  plant  perturbations  as  being 
largely  unstructured.  The  design  and  analysis  proceeds  for  the  worst  case 
situation  and  attempts  to  develop  aposteriori  bounds  for  the  perturbation. 

As  a  result,  the  designs  and  bounds  are  often  conservative.  The  second  area 
assumes  that  the  perturbations  can  be  modeled  apriori.  The  control  system 
is  then  developed  for  the  class  of  plants  described  by  the  perturbation  model. 
This  provides  hope  of  designing  less  conservative  (and  likely  more  efficient) 
control  systems.  However,  this  problem  appears  to  be  more  difficult  and 
much  less  research  has  been  conducted  in  this  area. 

The  research  presented  in  Chapters  4  through  6  of  this  report  is 
directed  at  the  second  area  described  in  the  preceding  paragraph.  The 
parameter  space  design  method  of  Chapter  4  provides  great  insight  into  the 
problem  for  the  designer.  Chapter  5  demonstrated  that  this  method  works  well 
for  single  input,  low  order  problems.  The  extension  of  this  technique  to 
higher  order  systems  appears  straightforward,  with  only  technical  problems 
to  be  overcome.  The  extension  to  multiple  input  systems  will  require  more 
thought.  One  possible  generalization  is  to  use  the  multivariable  pole 
placement  equation  [33]  which  exploits  the  characteristic  polynomial  matrix. 

The  optimization  procedure  of  Chapter  6  also  provides  considerable 
insight  for  the  designer.  It  also  allows  high  order,  multi-input  systems 


to  be  considered.  However,  the  computational  procedure  suffers  from  several 


technical  problems.  The  most  serious  of  these  problems  are  the  extreme 
sensitivity  of  closed  loop  eigenvalues  at  double  real  roots  and  the  require 
ment  that  an  initial  stabilizing  guess  be  supp’ led.  It  is  anticipated  that 
additional  computational  experience  with  the  algorithm  and  the  designer's 
insight  into  the  problem  can  alleviate  these  difficulties. 

In  conclusion,  it  should  be  emphasized  that  it  is  unlikely  that  a 
single  technique  can  be  developed  which  is  able  to  handle  all  problems.  It 
is  necessary  to  have  several  design  tools  available  and  to  use  the  method 
best  suited  to  the  problem  at  hand.  A  major  purpose  of  this  report  is  to 
present  methods  that  will  extend  the  domain  of  problems  which  can  be 
handled  by  control  system  engineers. 


APPENDIX  I 
AERODYNAMIC  DATA 

Table  1.1  contains  the  values  for  a^,  c..  defined  in  equations  (116). 
These  data  were  obtained  by  transforming  the  aero  data  in  [57]  to  be  compatible 
with  the  state  space  representation  of  equations  (116) . 

The  data  used  in  equations  (117)  are  shown  in  Table  1.2 

Table  1.1.  Aero  Data  for  Equations  (3.2) 


M  =  .5 

M-  .85 

M  =  .9 

M=  1.5 

500 O’ 

5000' 

35000' 

35000' 

>■ 

ai 

-  .8532 

-  1.514 

-  .6314 

-  .8527 

a2 

.9931 

.9940 

.9974 

.9982 

n 

a3 

-  .08756 

-  .1315 

-  .04332 

-  .04669 

•rf  ■ 

•S-1 

. 

a4 

0 

0 

0 

-  .02274 

a5 

4.641 

11.25 

1.488 

-18.50 

a 

a6 

-  .9876 

-  1.606 

-  .6680 

-  .8881 

U.  1 

- 

a7 

-10.25 

-26.15 

-  3.104 

-15.53 

a8 

4.246 

14.46 

4.590 

8.860 

■ 

- 

C1 

17.53 

51.11 

18.14 

26.83 

***■  ' 

C2 

-  .5152 

-  .8560 

-  .3576 

-  .4879 

C3 

-  5.078 

-12.95 

-  4.018 

-  7.842 

• 

c. 

2.723 

9.273 

2.944 

6.714 

- 

4 

177 


Table  1.2. 

Aero  Data 

for  Equations 

(3.3) 

M  =  .5 

M=  .85 

M  =  .9 

M  =  1.5 

5000' 

5000' 

35000' 

35000' 

all 

-  .9896 

-  1.702 

-  .6607 

.5162 

a12 

17.41 

50.72 

18.11 

26.96 

a13 

96.15 

263.5 

84.34 

178.9 

a21 

.2648 

.2201 

.08201 

.6896 

a22 

-  .8512 

-  1.418 

-  .6587 

-  1.225 

a23 

-11.39 

-  31.99 

-10.81 

-  30.38 

bl 

-97.78 

-272.2 

-85.09 

-175.6 
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APPENDIX  II 


A  TYPICAL  CONSTRAINT  MAPPING 


F~r  flight  condition  2  and  the  system  of  equations  (127),  the 
eigenvalue  constraint  regions  shown  in  Figure  II. 1  will  be  mapped  into  K-space. 


f  l  , ,  9 


I  \  III  \II 


a  =  - 3 . 5 


g=- 1 2 . 6 
f=-  50 


Figure  II. 1.  Eigenvalue  constraint  region  (not  drawn  to  scale) 


After  k^,  k^  were  arbitrarily  fixed  at  (-.03,0),  these  constraints  were 
mapped  into  the  regions  shown  in  Figures  II. 2.  Points  a'-g'  of  Figures  II. 2 
are  the  images  of  points  a-g  of  Figure  II. 1,  with  points  a',f',g'  corre¬ 
sponding  to  double  eigenvalues  at  a,f,g. 

The  boundaries  partition  K-space  into  regions  which  correspond  to 
the  number  of  eigenvalues  in  each  region  of  Figure  II. 1,  labeled  I,  II,  III. 
The  desired  combination  is: 

one  eigenvalue  in  region  I, 
two  eigenvalues  in  region  II, 
three  eigenvalues  in  region  III. 

Table  II. 1  lists  how  the  eigenvalues  are  distributed  among  regions  I,  II,  III 
for  gains  in  the  regions  labeled  A-M  in  Figure  II. 3. 


Table  II. 1.  Distribution  of  Eigenvalues 


Region  Number  in  I  Number  in  II  Number  in  III 


183 


REFERENCES 


H.  S.  Black,  "Stabilized  Feedback  Amplifiers,"  Bell  Systems  Technical 
Journal ,  1934,  pp.  1-18. 

J.  B.  Cruz,  Jr,  (Ed.),  "Effect  of  Feedback  on  Signal  Distortion  in  Nonlinear 
Systems,"  Chapter  3  of  Feedback  Systems,  McGraw-Hill,  New  York,  1972. 

C.  A.  Desoer  and  Y.  T.  Wang,  "Foundations  of  Feedback  Theory  for  Nonlinear 
Dynamical  Systems,"  IEEE  Trans,  on  Circuits  and  Systems,  1979,  to 
appear . 

H.  W.  Bode,  Network  Analysis  and  Feedback  Amplifier  Design,  D.  Van 
Nostrand  Company,  Inc.,  Princeton,  N.J.,  1945. 

J.  B.  Cruz,  Jr.  and  W.  R.  Perkins,  "A  New  Approach  to  the  Sensitivity 
Problem  in  Multivariable  Feedback  System  Design,"  IEEE  Trans,  on  Automatic 
Control ,  1964,  pp .  216-223. 

C.  H.  Hsu  and  C.  T.  Chen,  "A  Proof  of  the  Stability  of  Multivariable 
Feedback  Systems,"  Proc.  IEEE,  1968,  pp.  2061-2062. 

A.  G.  J.  MacFarlane,  "Return  Difference  and  Return-Ratio  Matrices  and  Their 
Use  in  Analysis  and  Design  of  Multivariable  Feedback  Control  Systems," 

Proc.  IFF ,  1970,  pp.  2037-2049. 

H.  H.  Rosenbrock,  Computer-Aided  Control  System  Design,  Academic  Press, 
London,  1974. 

J.  C.  Doyle,  "Robustness  of  Multiloop  Linear  Feedback  Systems,"  Proc.  1978 
IEEE  Conf-  on  Decision  and  Control,  San  Diego,  pp.  12-18. 

A.  J.  Laub ,  "Computational  Aspects  of  the  Singular  Value  Decomposition 
and  Some  Applications,"  Proc.  16th  Allerton  Conf.  on  Communication.  Control. 
and  Computing,  Univ.  of  Illinois,  October  1978,  pp.  432-442. 

N.  Sandell,  "Robust  Stability  of  Multivariable  Feedback  Systems,"  I^r qc_. 

16th  Allerton  Conf.  on  Communicat  ion ,  Control,  and  Comput  in  a; ,  Un  i  v  .  o  f 
Illinois,  October  197S,  pp.  471-479. 

G.  Stein  and  J.  C.  Doyle,  "Singular  Values  and  Feedback:  Design 
Examples,"  Proc.  16th  Allerton  Conf.  on  Communication ,  ContnyK  _and 
Computing,  October  1979.  pp.  461-470. 

M.  Safonov,  "Tight  Bounds  on  the  Response  of  Multivariable  Systems  with. 
Component  Uncertainty,  "  Proc.  16th  Allerton  Conf.  on  Communication,  Cent  r  ' '  . 
and  Computing,  October  1978,  pp.  431-460. 

R.  E.  Kalman,  '"..hen  is  a  Linear  Control  System.  Optimal  ?"  Tran-  . 
yj .  Basic  Eng.),  1964 ,  pp .  51-60. 


i 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS- 1963  A 


B.  D.  0.  Anderson  and  J.  B.  Moore,  Linear  Optimal  Control,  Prentice-Hall 
Englewood  Cliffs,  N.J.,  1971. 


M.  Safonov  and  M.  Athans,  "Gain  and  Phase  Margin  for  Multiloop  LQG 
Regulators,"  IEEE  Trans,  on  Automatic  Control,  1977,  pp.  173-197. 

0.  J.  M.  Smith,  Feedback  Control  Systems,  McGraw-Hill,  New  York,  1958. 

J.  B.  Doyle,  "Guaranteed  Margins  for  LQG  Regulators,"  IEEE  Trans,  on 
Automatic  Control,  1978,  pp.  755-756. 

J.  C.  Doyle  and  A.  Stein,  "Robustness  with  Observers,"  Proc.  1978  IEEE 
Conf.  on  Decision  and  Control,  San  Diego,  pp.  1-6. 

H.  Kwakernaak,  "Optimal  Low  Sensitivity  Linear  Feedback  Systems," 
Automatica,  1969,  p.  279. 

Jacques  Willems  and  H.  van  der  Voorde,  "The  Return  Difference  for  Discrete- 
Time  Optimal  Feedback  Systems,"  Automatica,  1978,  pp.  511-513. 

C.  A.  Harvey  and  R.  E.  Pope,  "Study  of  Synthesis  Techniques  for  Insen¬ 
sitive  Aircraft  Control  Systems,"  NASA  Contractor  Report  CR-2803,  April 
1977. 

C.  A.  Harvey  and  R.  E.  Pope,  "Insensitive  Control  Technology  Development," 
NASA  Contractor  Report  2947,  February  1978. 

A.  Vinkler  and  L.  Wood,  "A  Comparison  of  Several  Techniques  for  Designing 
Controllers  of  Uncertain  Dynamic  Systems,"  Proc,  1978  IEEE  Conf,  on 
Decision  and  Control,  San  Diego,  pp.  31-38. 

D.  M.  Salmon,  "Minimax  Controller  Design,"  IEEE  Trans,  on  Automatic 
Control,  1968,  pp.  369-373. 

U.  Ly  and  R.  H.  Cannon,  "A  Direct  Method  for  Designing  Robust  Optimal 
Control  Systems,"  Proc.  AIAA  Guidance  and  Control  Conf.,  Palo  Alto,  Calif. 
August  1978,  pp.  440-448. 

A.  Vinkler  and  L.  J.  Wood,  "Guaranteed  Cost  Control  of  Linear  Systems 
with  Uncertain  Parameters — Application  to  Remotely  Piloted  Vehicle 
Flight  Control  Systems,"  Proc.  AIAA  Guidance  and  Control  Conf.,  Palo 
Alto,  Calif.,  August  1978,  pp.  226-239. 

D.  L.  Kleinman  and  P.  K.  Rao,  "An  Intormation  Matrix  Approach  for  Aircraft 
Parameter  Insensitive  Control,"  Proc.  1977  IEEE  Conf.  on  Decision  and 
Control ,  New  Orleans,  La.,  pp.  316-325. 


K.  K.  D.  Young,  P.  Kokotovic,  and  V.  Utkin,  "A  Singular  Perturbation 
Analysis  of  High-Gain  Feedback  Systems,"  IEEE  Trans,  on  Automatic 
Control,  1977,  pp.  931-938. 


185 


30.  V.  Utkin,  "Variable  Structure  Systems  with  Sliding  Modes,"  IEEE  Trans,  on 
Automatic  Control,  1977,  pp.  212-222. 

31.  K.  K.  D.  Young,  "Design  of  Variable  Structure  Model  Following  Control 
Systems,"  IEEE  Trans,  on  Automatic  Control,  1978,  pp.  1079-1085. 

32.  E.  J.  Davison,  "The  Robust  Decentralized  Control  of  a  General  Servo¬ 
mechanism  Problem,"  IEEE  Trans,  on  Automatic  Control,  1976,  pp.  14-24. 

33.  J.  Ackermann,  Abtastregelung,  Springer,  Berlin,  1972. 

34.  J.  Ackermann,  "Parameter  Space  Design  of  Robust  Control  Systems," 
preliminary  version,  Proc.  of  the  JACC,  Denver,  Colo.,  June  1979. 

35.  N.  Franklin  and  J.  Ackermann,  "Robust  Flight  Control — A  Design  Example," 
in  preparation. 

36.  A.  T.  Fam  and  J.  S.  Meditch,  "A  Canonical  Parameter  Space  for  Linear 
Systems  Design,"  IEEE  Trans,  on  Automatic  Control,  1978,  pp.  454-458. 

37.  A.  Schy,  "Nonlinear  Programming  in  Design  of  Control  Systems  with 
Specified  Handling  Qualities,"  Proc.  1972  IEEE  Conf.  on  Decision  and 
Control,  New  Orleans,  La. 

38.  A.  Schy,  W.  M.  Adams,  and  K.  A.  Johnson,  "Computer  Aided  Design  of  Control 
Systems  to  Meet  Many  Requirements,"  AGARD  Conf.  Proc.  on  Advances  in 
Control  Systems,  Nr.  137,  Gailo,  Norway,  1973,  pp.  6. 1-6. 7. 

39.  F.  D.  Hauser,  "A  Nonlinear  Programming  Algorithm  for  Automated  Design 
and  Optimization  of  Flexible  Space  Vehicle  Autopilots,"  AIAA  Guidance 
and  Control  Conf,,  Key  Biscayne,  Fla.,  1973. 

40.  J.  S.  Karmarkar,  "A  Regulator  Design  by  Mathematical  Programming  Methods," 
Proc.  of  the  JACC.  1973,  pp.  699-710. 

41.  A.  Kanarachos,  "Computer  Aided  Design  of  Control  Loops  by  Parameter 
Optimization  Methods, "(in  German),  Regelungs technik ,  1978,  pp.  220-226. 

42.  A.  Kreisselmeier  and  R.  Steinhauser,  "Insensitive  Control  for  Large 
Parameter  Variations  for  a  Stabilizer  for  the  F-4C  Aircraft,"  to  appear 
as  DFVLR-Forschungsbericht . 

43.  A.  Kreisselmeier  and  R.  Steinhauser,  "Systematic  Control  Design  by 
Optimizing  a  Vector  Performance  Index,"  Proc.  IFAC  Symp.  on  Computer- 
Aided  Design  of  Control  Systems,  Zurich,  Switzerland,  August  1979. 

44.  J.  J.  Belletrutti  and  A.  G.  J.  MacFarlane,  "Characteristic  Loci  Techniques 
in  Multivariable  Control  System  Design,"  Proc.  IEE,  1971,  pp.  1291-1297. 

45.  D.  H.  Ovens ,  "Integrity  of  Multivariable  First-Order-Type  Systems,"  Int ■ 

J.  of  Control,  1976,  pp.  827-835. 


186 


46. 


47. 


48. 


49. 


50. 


51. 


53. 


54. 


55. 


56. 


57. 


58, 


59. 


60. 


O.  A.  Solheim,  "Some  Integrity  Problems  in  Optimal  Control  Systems,"  AGARD 
Conf.  Proc.  on  Advances  In  Control  Systems,  Nr.  137,  Gello,  Norway,  1973, 
pp.  4.4-4.10. 

P.  K.  Wong,  A.  Stein,  and  M.  Athans,  "Structural  Reliability  and 
Robustness  Properties  of  Optimal  Linear-Quadratic  Multivariable 
Regulators,"  Proc.  7th  IFAC  Congress,  Helsinki,  Finland,  1978,  pp.  1797- 
1805. 


A.  Kreisselmeier ,  "Considerations  on  the  Robustness  of  Control  Systems," 
European  Space  Agency,  Technical  translation  ESA-TT-453,  May  1978, 
(translation  of  DFVLR- IB-552-77/19) . 

D.  D.  Siljak,  Nonlinear  Systems,  The  Parameter  Analysis  and  Design, 

J.  Wiley,  New  York,  1969,  Chapts.  1  and  2. 

J.  Rissanen,  "Control  System  Synthesis  by  Analogue  and  Computer  Based  on 
the  Generalized  Linear  Feedback  Concept,"  Inti.  Seminar  of  Analog 
Computation  Applied  to  the  Study  of  Chemical  Processes,  Brussels, 

Presses  Academlques  Europeennes,  Nov.  1960. 

E.  J.  Davison,  "On  Pole  Assignment  in  Linear  Systems  with  Incomplete 
State  Feedback,"  IEEE  Trans,  on  Automatic  Control,  Vol.  AC-15,  1970, 
pp.  348-351. 

M.  Marden,  Geometry  of  Polynomials,  American  Mathematical  Society, 
Providence,  R.I.,  1966. 


H.  R.  Sirisena  and  S.  S.  Choi,  "Pole  Placement  in  Prescribed  Regions  of 
the  Complex  Plane  Using  Output  Feedback,"  IEEE  Trans,  on  Automatic  Control, 
Vol.  AC-20,  1975,  pp.  810-812. 


S.  N.  Franklin,  "Design  of  a  Robust  Flight  Control  System,"  M.S.  thesis 
University  of  Illinois,  Urbana,  1979. 

J.  Ackermann,  "On  the  Synthesis  of  Linear  Control  Systems  with  Specified 
Characteristics,"  Automatics,  Vol.  13,  1977,  pp.  89-94. 


V.  M.  Popov,  "Invariant  Description  of  Linear  Time-Invariant  Controllable 
Systems,"  SIAM  J.  Control,  Vol.  10,  1972,  pp.  252-264. 

R.  L.  Berger,  J.  R.  Hess,  and  D.  C.  Anderson,  "Compatibility  of  Maneuver 
Load  Control  and  Relaxed  Static  Stability  Applied  to  Military  Aircraft," 
AFFDL-TR-73-33,  April  1973. 

"Flying  Qualities  of  Piloted  Airplanes,"  MIL-F-8785B  (ASG) ,  August  7,  1969. 

D.  McRuer,  I.  Ashkenas,  and  D.  Graham,  Aircraft  Dynamics  and  Automatic 
Control,  Princeton  University  Press,  Princeton,  N.J.,  1973. 


D.  P.  Looze  and  N.  R.  Sandell,  Jr.,  "Gradient  Calculations  for  Linear 
Quadratic  Fixed  Structure  Control  Problems,"  IEEE  Trans,  on  Automatic 
Control,  Vol.  AC-25,  No.  2,  April  o980.  9 


187 


1 

k 

y- 

l 


61.  D.  G.  Luenberger,  Introduction  to  Linear  and  Nonlinear  Programming, 
Addison-Wesley ,  Reading,  Mass.,  1973. 

62.  M.  R.  Hestenes,  "Multiplier  and  Gradient  Methods,"  J.  Optimization  Theory 
and  Applications,  Vol.  4,  No.  5,  1969. 

63.  D.  P.  Bertsekas,  "Combined  Primal-Dual  and  Penalty  Methods  for 
Constrained  Minimization,"  SIAM  J.  Control,  Vol.  13,  No.  3,  May  1975. 

64.  D.  P.  Bertsekas,  "On  the  Method  of  Multipliers  for  Convex  Programming," 
IEEE  Trans,  on  Automatic  Control,  Vol.  AC-20,  No.  3,  June  1975. 

65.  J.  E.  Dennis  and  H.  H.  W.  Mei,  "An  Unconstrained  Optimization  Algorithm 
which  Uses  Function  and  Gradient  Values,"  Technical  Report  TR  75-246, 
Dept,  of  Computer  Science,  Cornell  University,  Ithaca,  N.Y.,  June  1975. 

66.  J.  B.  Mankin,  A  Study  of  Control  System  Sensitivity,  Management 
Information  Services,  Detroit,  Mich.,  1972. 

67.  D.  G.  Luenberger,  Optimization  by  Vector  Space  Methods ,  Wiley,  New  York, 
1969. 

68.  H.  Kwakernaak  and  R.  Sivan,  Linear  Optimal  Control  Systems,  Wiley- 
Interscience,  New  York,  1972. 


